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1 .  Introduction 

My  objective  in  this  paper  is  to  give  some  of  the  basic  results  in  the  theory  of  bifurca¬ 
tion  in  differential  equations.  It  is  difficult  to  trace  the  historical  development  of  any  im¬ 
portant  concept  and  bifurcation  theory  is  no  exception.  However,  a  careful  study  of  litera¬ 
ture  shows  that  Poincare  [1],  [2]  and  Liapunov  [1],  [2]  are  responsible  for  our  present 
basic  philosophy  as  well  as  several  of  the  fundamental  ideas  of  the  methods  that  we  pres¬ 
ently  employ.  These  two  persons  can  be  directly  linked  with  the  importance  of  exchanges 
of  stability,  the  occurrence  of  complicated  motions  in  dynamical  systems,  the  principle  of 
reduction  to  lower  dimensional  problems,  the  philosophy  of  genericity  and  the  transforma¬ 
tion  theory  of  differential  equations  that  is  so  important  in  obtaining  approximations  of  the 
center  manifold  and  the  flow  on  the  center  manifold,  in  many  respects,  we  are  still  exploit¬ 
ing  the  ideas  of  these  two  giants. 

After  the  initial  impetus  of  Poincard  and  Liapunov,  it  is  somewhat  surprising  that 
modern  bifurcation  theory  did  not  appear  at  an  earlier  date.  It  is  perhaps  true  that  the 
ideas  of  Liapunov  connected  with  bifurcation  theory  were  being  developed  more  extensively 
than  the  corresponding  ones  of  Poincard.  There  was  a  very  active  group  in  the  U.S.S.R. 
(consisting  of  Andronov,  Vitt,  Khaikin,  Bogoliubov,  Krylov,  Leontovich,  Malkin  and  others) 
working  on  critical  cases  in  stability  theory,  nonlinear  oscillations  and  the  general  theory  of 
integral  manifolds.  The  techniques  developed  from  the  study  of  these  areas  are  fundamental 
ingredients  in  dynamic  bifurcation  theory  (see  Hale  [I]  for  some  references). 

A  fundamental  step  towards  modern  bifurcation  theory  in  differential  equations  oc¬ 
curred  with  the  definition  of  structural  stability  of  Andronov  and  Pontrjagin  [1]  in  1937 
and  the  classification  of  structurally  stable  systems  in  the  plane.  With  these  concepts.  An¬ 
dronov  and  Leontovich  [1]  were  able  to  make  precise  definitions  of  types  of  bifurcation 
points  which  had  the  possibility  of  being  classified.  These  results  were  applied  extensively 
to  the  theory  of  nonlinear  oscillations  by  Andronov,  Vitt  and  Khaikin  in  1937  (the  second 
edition  of  this  book  is  Andronov,  Vitt  and  Khaikin  [I  j ).  As  Minorsky  [1 J  said  in  1962: 
"Having  established  the  initial  advance  in  the  field  of  nonlinear  oscillations  (up  to  1940),  the 
Russian  scientists  maintain  their  leadership  and  initiative  characterized  by  a  remarkable  co¬ 
ordination  of  efforts  between  the  mathematical  and  experimental  parts  of  these  fundamental 
researches.”  There  were  several  important  developments  in  this  intervening  period  by  Levin¬ 
son  [1),  [2],  Cartwright  and  Littlewood  (see  Cartwright  [I])  on  the  forced  van  der  Pol  and 
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Lienard  equations.  However,  in  western  Europe  and  the  United  States,  the  interest  in  this 
aspect  of  differential  equations  had  never  been  very  extensive.  In  addition,  there  was  little 
awareness  of  the  developments  that  had  been  made  in  the  U.S.S.R.,  and,  as  a  consequence, 
some  duplication  of  effort  occurred. 

Since  1960,  there  have  been  extensive  developments  in  the  abstract  theory  of  dynami¬ 
cal  systems.  At  the  same  time,  some  of  this  theory  has  been  applied  to  very  interesting 
problems  in  the  biological  and  physical  sciences.  In  an  attempt  to  explain  phenomena  that 
occur  in  nature,  it  has  been  necessary  for  researchers  to  discuss  the  dynamic  bifurcation  of 
specific  types  of  equations  in  great  detail.  This  has  led  to  an  exciting  interaction  between 
analytical  and  theoretical  methods. 

In  this  paper,  we  present  some  of  the  concepts  and  results  that  play  an  important  role 
in  these  areas.  When  the  dimension  of  the  system  is  one  or  two,  one  can  obtain  a  rather 
complete  theory  at  least  for  general  one  parameter  families  of  vector  fields.  For  either 
several  parameter  problems  or  for  the  dimension  of  the  system  greater  than  two,  only  par¬ 
tial  results  are  known.  On  the  other  hand,  the  results  in  low  dimension  are  applicable  to 
higher  dimensional  problems  (even  infinite  dimensional  ones)  when  the  discussion  is  restricted 
to  a  neighborhood  of  an  equilibrium  point  for  which  the  theory  of  center  manifolds  can  be 
employed . 

The  table  of  contents  expresses  in  general  terms  the  substance  of  this  paper.  The  first 
eight  sections  deal  with  structural  stability  and  bifurcation  in  the  low  dimensional  problems 
mentioned  above.  §9  is  devoted  to  the  formulation  of  some  of  the  basic  problems  in  the 
qualitative  theory  for  a  special  class  of  dynamical  systems  in  infinite  dimensions.  This  class 
is  general  enough  to  include  many  functional  differential  equations  and  partial  differential 
equations.  §10  is  concerned  primarily  with  some  types  of  bifurcation  that  occur  because 
the  base  space  is  infinite  dimensional.  Due  to  space,  very  few  proofs  are  given.  Also,  there 
are  several  important  omissions  of  topics  from  differential  equations  that  are  systematically 
used  throughout,  but  which  are  not  as  well  known  as  they  should  be.  Most  notable  among 
these  are  the  general  theory  of  integral  manifolds  (for  general  references,  see  Hale  [1]).  the 
center  manifold  theorem  (see  Kelley  [1]),  the  theory  of  transformation  to  normal  forms 
(for  references,  see  Bibikov  [1] ,  Br 'juno  [1] ,  Henrard  [1])  and  the  general  method  of 
averaging  (for  references,  see  Hale  [1]). 

The  author  is  endebted  to  many  colleagues  and  students  whose  ideas  have  been  incor¬ 
porated  into  these  notes-too  many  to  mention  by  name.  He  also  acknowledges  the  initia¬ 
tive  of  Professor  Laksmikantham  in  proposing  the  CBMS  Regional  Conference.  Finally, 
Sandra  Spinacci  has  exhibited  her  usual  patience  and  understanding  in  the  preparation  of  the 
final  manuscript. 


k 


2.  On  the  definition  of  bifurcation 

Suppose  X,  Z  are  topological  spaces,  U  C  X  is  open,  A  is  an  open  set  in  a  topological 
space  and  /:  U  x  A  — *  Z  is  a  given  continuous  function.  Let 

S  =  {(*.  X)  G  U  x  A:  /(*,  X)  =  0} 
be  the  set  of  solutions  of  the  equation  fix,  X)  =  0.  For  a  fixed  X,  let 

Sx  =  {*:  (x,  X)  e  S } 

be  the  “cross-section”  of  the  solution  set  at  X. 

The  basic  problem  is  to  discuss  the  dependence  of  the  set  Sx  on  X.  In  a  specific  prob¬ 
lem,  one  has  a  prescription  which  compares  Sx  with  for  different  X  and  p.  This  com¬ 
parison  is  usually  made  by  means  of  an  equivalence  relation  which  divides  the  sets  {$x.  A  £  A} 
into  equivalence  classes.  Given  the  function /and  an  equivalence  relation  ~,  we  say  X0  is  a 
bifurcation  point  for  (/,  ~)  if,  for  any  neighborhood  V  of  X0,  there  are  X, ,  X2  £  V  such 
that  SX|  -f  V  This  definition  is  less  general  than  the  one  in  Marsden  [1] . 

A  special  case  is  when  the  equivalence  relation  specifies  that  Sx  ~  if  the  sets  Sx 
and  SM  are  homeomorphic.  This  is  a  very  convenient  choice  when  studying  the  change  in 
the  structure  of  the  set  of  equilibrium  points  in  a  differential  equation  as  parameters  are 
varied.  In  this  case,  the  function  /  represents  the  vector  field  in  a  differential  equation 
dx/dt  =  fix,  X).  It  is  also  appropriate  in  differential  equations  for  the  study  of  the  set  of 
solutions  of  some  prescribed  type;  for  example,  periodic  solutions,  invariant  tori,  etc.  In 
this  latter  case,  the  topological  spaces  are  defined  so  that  they  include  only  those  functions 
which  exhibit  this  prescribed  behavior  and  the  function  /  could  be  the  differential  operator. 
fix,  X)  =  dx/dt  -  gix,  X). 

To  study  more  general  bifurcations  in  differential  equations,  the  equivalence  defined 
by  homeomorphism  is  not  sufficient.  Consider  a  differential  equation  du/dt  -  giu,  X)  =  0 
where  (u,  X)  £  SI  x  A  and  £2  is  an  open  set  in  some  Banach  space  E.  For  X,  Z  Banach  spaces 
of  functions  from  [0,  °»)  to  E,  let  U  C  X  be  defined  by  U  =  {u  £  X:  u(t )  £  S2,  r  £  [0.  °°)} . 
The  above  equation  can  be  written  formally  as  f(u,  X)  =  0  where  /:  U  x  A  — *■  Z,  f{u.  XHO 
=  du(t)jdt  -  g{u(t),  X).  Assuming  everything  can  be  made  rigorous  and  that  all  solutions 
are  obtained  in  this  way,  a  comparison  of  the  corresponding  sets  Sx  and  by  homeomor¬ 
phism  will  not  be  very  interesting.  Thus  an  alternative  approach  must  be  taken. 

Suppose  the  differential  equation  generates  a  strongly  continuous  semigroup  Tx(t), 
t  >  0,  on  n.  A  frequently  used  concept  of  equivalence  in  differential  equations  is  to  say 
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that  £(• ,  X)  ~  g(- ,  p)  if  there  is  a  homeomorphism  h:  fi  — ►  n  such  that  h  maps  orbits  of 
Tx(t)  onto  orbits  of  TM(t)  preserving  the  sense  of  direction  in  time.  A  vector  field  £(■.  X0) 
is  structurally  stable  if  there  is  a  neighborhood  V  of  X0  such  that  g(- ,  X)  ~  g{- ,  X0)  for  all 
X  €  V.  Thus,  X0  is  a  bifurcation  point  if  X0  is  not  structurally  stable. 

A  different  but  equivalent  formulation  of  the  above  concept  of  equivalence  in  differen¬ 
tial  equations  was  introduced  by  Andronov  and  Pontrjagin  [1]  in  1937  for  differential  equa¬ 
tions  in  the  plane.  They  gave  a  characterization  of  the  structurally  stable  vector  fields  which 
will  be  discussed  later.  Peixoto  [1]  generalized  these  results  to  arbitrary  compact  two  di¬ 
mensional  manifolds  and  proved  the  set  of  structurally  stable  vector  fields  is  open  and  dense. 
For  some  time,  it  was  the  feeling  that  this  same  property  should  hold  true  for  arbitrary  sys¬ 
tems.  Unfortunately,  it  was  shown  by  Smale  [1]  that  structurally  stable  vector  fields  are 
not  dense  in  dimension  >  4.  Williams  [1]  proved  the  same  result  for  n  >  3.  Since  many 
vector  fields  cannot  be  compared  by  this  equivalence  relation,  it  becomes  necessary  to 
weaken  the  concept  of  equivalence.  Each  new  definition  of  equivalence  leads  to  a  new  type 
of  stable  vector  field  (ones  which  are  equivalent  to  everyone  in  a  neighborhood  of  it)  and 
thus  a  new  type  of  bifurcation.  The  ultimate  goal  is  to  have  the  definition  restrictive  enough 
to  permit  classification  of  the  stable  ones,  but,  at  the  same  time,  to  have  the  stable  vector 
fields  generic;  that  is,  the  intersection  of  a  countable  sequence  of  open  dense  sets.  Much  of 
the  research  in  finite  dimensional  abstract  dynamical  systems  in  the  last  twenty  years  has 
been  devoted  to  this  general  problem.  Relevant  references  are  Smale  [2] ,  Peixoto  [2] ,  [3] . 
Palis  and  Melo  [l],Newhouse  [1  ] ,  Nitecki  f  1 J ,  Shub  [IJ.  Guckenheimer  [IJ.Arnofd  [1]. 
In  the  next  section,  we  give  more  specific  details. 

When  the  evolutionary  equation  is  infinite  dimensional,  several  new  problems  arise. 

This  case  will  be  discussed  in  a  later  section. 


3.  Structural  stability  and  generic  properties  in  R" 

Suppose  SI  is  an  open  set  in  R”  with  dSl  =  T,  f2  =  SI  U  T. 

The  space  C'Xft,  R")  is  the  Banach  space  of  functions  bounded  and  continuous  to¬ 
gether  with  all  derivatives  up  through  order  r  >  0  with  the  norm  of  /in  Cr(S2,  R")  being 
given  by  the  maximum  of  the  supremum  over  SI  of  the  norm  of  /and  its  derivatives  up 
through  order  r.  Let  X'  =  Xj,(S2)  be  the  set  of  elements  of  C(Sl ,  R")  which  are  transver¬ 
sal  to  T.  For  any  /G  X' ,  r  >  1 ,  the  differential  equation 

(3-D  *=/(*) 

defines  a  family  of  transformations  Tf(t)  on  SI  satisfying  the  semigroup  property  with 
Tf(t)x0  =  xft .  x0),  where  x(r,  x0)  is  the  solution  of  (3.1)  with  x(0.  x0)  =  x0.  Further¬ 
more,  for  each  x0  G  SI,  there  are  an  aXQ  <  0,  (JXQ  >  0,  such  that  the  maximal  interval  of 
definition  of  Tf(t)x0  is  [aJtQ,  PXQ).  The  number  is  either  that  value  aXQ  where 
T/(ax 0)x0  G  an  =  T  or  and  in  this  case  the  interval  [arX(),  PXQ)  is  (~°°,  PXQ)-  The 
number  PXq  is  defined  in  a  similar  way  in  the  positive  direction.  The  operator  Tf(t)  on  SI 
satisfies  7^(0)  =  I,  the  identity,  Tf(t  4-  s)x  =  Tf(t)Tf(s)x  for  those  t,  s  for  which  it  is 
meaningful  and  Tfft)x  has  continuous  derivatives  up  through  order  r  in  t,  x. 

The  orbit  y^(x)  of  f  through  x  is 

ff(x)  =  U { Tf{t)x,  te[ax,px]}. 

The  co-limit  set  uy(x)  and  a-limit  set  a^(x)  of  the  orbit  tf(x)  are  defined  by 
cof(x)  =  PI  cl  U  Tf(t)x,  af(x)  =  PI  cl  U  Tf(t)x. 

r>0  t>r  tCO  r<r 

An  equilibrium  point  or  critical  point  of  /  is  a  zero  of  /.  A  periodic  orbit  of  /  is  an  orbit 
which  is  a  closed  curve.  A  set  M  C  n  is  invariant  if,  for  each  x  G  M,  T^(t)x  is  defined  for 
t  G  (-oo,  oo)  and  belongs  to  M  for  t  G  (-<»,  «).  This  implies  Tf(t)M  =  M  for  t  G  (-  °°.  »). 

The  vector  fields  in  (3.1)  are  chosen  from  X J, ;  that  is,  are  transversal  to  T,  in  order 
to  eliminate  technical  difficulties  with  points  of  contact  on  T.  We  are  discussing  the  vector 
fields  in  R" ,  but  many  of  the  remarks  hold  for  vector  fields  on  compact  manifolds  M. 

Definition  3.1 .  Two  vector  fields  /  g  in  X^ ,  r  >  1 ,  are  equivalent,  f  —  g,  if  there 
is  a  homeomorphism  h:  Ti  —*  SI  such  that  h  maps  the  orbits  defined  by  /  homeomorphic- 
ally  onto  the  orbits  defined  by  g  with  the  sense  of  direction  in  time  preserved.  An  / G  Kr„ 
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is  said  to  be  structurally  stable  if  there  is  a  neighborhood  U  of  /such  that  f  ~  g  for  every 
g  £  U.  An  /€  X„  is  a  bifurcation  point  if  / is  not  structurally  stable. 

Two  important  remarks  need  to  be  made  about  this  definition.  Definition  3.1  would 
not  be  meaningful  without  the  condition  r>\.  In  fact,  for  r  =  0,  given  any  vector  field  / 
that  has  an  isolated  zero  at  x0  and  any  e  >  0,  there  are  a  6  >  0  and  a  function  g  such  that 
\f-g\  <  e  and  gfx)  =  0  for  |x  -  x0|  <  6.  Therefore,  no  f  with  an  isolated  zero  could  be 
structurally  stable. 

In  Definition  3.! ,  it  is  tempting  to  require  that  the  mapping  h  be  a  diffeomorphism. 
However,  if  /( 0)  =  0,  g(0)  =  0,  3/(0)/3x  =  A.  dg(0)/dx  =  B,  and  /  ~  g  in  a  neighborhood 
of  zero,  then  one  can  show  (see  Peixoto  [2j ,  (4J)  that  the  eigenvalues  of  A  and  B  must  be 
proportional.  Since  one  can  always  make  a  small  perturbation  that  will  change  one  eigen¬ 
value  of  A  and  not  the  other,  it  follows  that  no  vector  field  with  a  zero  could  be  structurally 
stable.  Thus,  the  Definition  3.1  would  have  little  meaning.  If  x0  is  a  critical  point  of /and 
A  =  3/(x0)/3x,  then  x0  is  said  to  be  hyperbolic  if  the  real  parts  of  the  eigenvalues  of  A 
have  nonzero  real  parts.  The  point  x0  is  a  saddle  point  of  order  k,  if  it  is  hyperbolic  and 
there  are  k  eigenvalues  of  A  with  positive  real  parts.  The  term  saddle  point  without  the 
designation  of  the  order  will  refer  to  any  saddle  point  of  order  k  with  k  &  0  or  n. 

If  n  =  2,  a  saddle  point  of  order  1  corresponds  to  the  usual  definition  of  saddle  point. 
For  n  =  2,  a  saddle  point  of  order  0  or  2  corresponds  to  a  node  or  focus  depending  upon 
whether  the  eigenvalues  of  A  are  real  or  complex. 

If  7  is  a  periodic  orbit  of  /,  then  one  can  define  a  Poincare  map  near  7  in  the  following 
way.  For  any  arc  C  transversal  to  7  at  p0  and  any  p  £  C  sufficiently  near  p0 ,  there  is  a 
unique  r(p)  >  0  such  that  T^(t{ p))p  e  C,  Tf(t)p  €  C  for  0  <  t  <  T(p).  The  map  p  <-» 
7/(r(p))p  is  called  the  Poincare  map  7r(p).  This  map  in  C  and  7r(p0)  =  p0.  The  periodic 
orbit  7  is  hyperbolic  if  no  eigenvalue  of  37r(p0)/3p  has  modulus  one. 

If  n  =  2,  the  periodic  orbit  7  is  hyperbolic  if  dir(p0)/dp  =£  1 .  It  is  instructive  to  give 
an  equivalent  definition  in  terms  of  the  vector  field  itself.  If  7  =  {0(f).  f  £  R]  where  0(f) 
is  periodic  of  least  period  to  and  0(f)  =  /(0(f))  then  the  linear  variational  equation  for  0  is 

(3.2)  y  =  A(i)y.  A{t)  =  3/(0(f))/dx. 

One  characteristic  multiplier  of  this  co-periodic  system  is  1  since  0  satisfies  (3.2).  If  X{t)  is 
a  principal  matrix  solution  of  (3.2),  then  the  product  of  the  multipliers  is  equal  to  det  X(u>). 
Thus,  if  py  =  exp  a >oy,  oy  real,  is  the  other  multiplier,  then 

(3.3)  a7  =  ifoUA(s)ds- 

One  can  then  easily  show  that  7  is  hyperbolic  if  and  only  if  oy  0,  unstable  if  oy  >  0  and 
asymptotically  orbitally  stable  if  ay  <0. 

For  two  dimensional  systems,  the  following  result  of  Andronov  and  Pontrjagin  [1] . 
Peixoto  [1],  completely  solves  the  problem  of  structural  stability  in 
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Thkorkm  3.2.  //  C  X  j,  r  >  I ,  is  the  set  of  structurally  stable  vector  fields  in  X2. 
then  f  €  Z,r2  if  arul  only  if  the  following  conditions  are  satisfied: 

(i)  The  critical  points  of  fare  hyperbolic. 

(ii)  The  periodic  orbits  of  f  are  hyperbolic. 

(iii)  There  is  no  orbit  of  f  with  both  the  a -  and  urlimit  sets  being  saddle  points. 

Furthermore.  is  open  and  dense  in  Xr2. 

The  fact  that  a  structurally  stable  vector  field  must  satisfy  (i)-(iii)  is  very  easy  to 
prove.  However,  the  converse  is  more  difficult  and  relies  heavily  upon  the  following  result 
of  Hartman  [1J,  (3]  Grobman  [lj,  and  its  extension  to  diffeomorphisms  which  is  valid  in 
the  space  of  n -dimensional  vector  fields  X  „ . 

Thkorkm  3.3  (Hartman-Grobman).  Iff  S  X'n.  r>  l,  f(x0)  =  0,  and  the  eigen¬ 
values  of  A  =  df(x0)/dx  have  nonzero  real  parts  then,  in  a  neighborhood  of  x0,  f  is  equiva¬ 
lent  to  the  linear  equation  x  -  Ax. 

In  Theorem  3.2,  the  fact  that  Xj  is  open  follows  from  the  definition  and  the  fact  that 
it  is  dense  follows  from  an  argument  in  transversality  theory.  See  Peixoto  [1]  for  a  com¬ 
plete  proof. 

Condition  (i),  the  Implicit  Function  Theorem  and  the  compactness  of  n  imply  that 
/€  Xj  has  only  a  finite  number  of  critical  points.  Using  (ii),  (iii)  and  similar  arguments,  one 
shows  there  is  only  a  finite  number  of  periodic  orbits. 

The  simplicity  of  the  description  of  the  structurally  stable  systems  in  two  dimensions 
given  by  Theorem  3.2  permits  a  complete  classification  in  terms  of  certain  distinguished 
graphs  (see  Peixoto  [4| ). 

To  what  extent  does  Theorem  3.2  hold  in  dimension  n  >  3?  As  remarked  earlier,  the 
structurally  stable  systems  are  not  dense  in  X'  for  n  >  3.  This  was  proved  by  Smale  |1 ) 
for  n  >  4  and  by  Williams  [1]  for  n  3*  3.  However,  there  are  structurally  stable  systems  in 
every  dimension  and  on  every  type  of  ^-dimensional  manifold. 

Even  though  is  not  dense,  it  is  very  important  to  classify  structurally  stable  vector 
fields  and  to  find  “simple”  classes  of  vector  fields  which  are  generic.  Let  us  turn  first  to 
the  problem  of  genericity. 

The  concepts  (i),  (ii)  in  Theorem  2  have  meaning  in  R".  Also,  (iii)  can  be  ex¬ 
tended  in  the  following  way.  For  any  hyperbolic  critical  point  or  periodic  orbit  of  a 
vector  field  /  6  Xrn,  one  can  define  the  global  stable  and  unstable  manifolds  in  the  fol¬ 
lowing  way.  The  stable  (unstable)  manifold  of  a  hyperbolic  critical  point  x0  is  the  set 
of  x  €  n  such  that  Tj(t)x  — *  jc0  as  I  — ►  +°°  (-<»).  Similar  definitions  are  given  for  a 
periodic  orbit. 

In  R2,  condition  (iii)  is  then  equivalent  to  the  statement  that  the  stable  and  unstable 
manifolds  of  all  critical  points  and  periodic  orbits  intersect  transversaily.  One  can  then  ask 
if  the  vector  fields  in  Xr„  which  satisfy  these  properties  are  generic  in  X'.  The  answer  is 
yes  and  is  the  famous  theorem  of  Kupka  (1J  and  Smale  [3j. 
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Thuokkm  3.4  (Kui’ka-Smalk).  The  set  of  vector  fields  in  Xr„  for  which  the  critical 
points  and  periodic  orbits  are  hyperbolic  with  stable  and  unstable  manifolds  intersecting 
transversally  is  generic. 

Any  vector  field  satisfying  the  conditions  of  Theorem  3.4  will  be  called  a  Kupka-Smale 
(KS)  vector  field.  They  can  have  only  a  finite  number  of  critical  points  with  the  proof  being 
the  same  as  in  two  dimensions.  However,  in  contrast  to  two  dimensions,  there  can  be  an 
infinite  number  of  periodic  orbits  if  the  dimension  is  >  3  (for  an  example,  see  Nitccki  [1) , 
Palis  and  deMelo  [1]). 

The  KS  vector  fields  are  dense,  but  all  KS  vector  fields  cannot  be  structurally  stable 
since  the  structurally  stable  systems  are  not  dense  in  dimension  >  3.  To  find  a  subset  of 
the  KS  vector  fields  which  are  structurally  stable,  one  must  put  some  further  restrictions  on 
the  behavior  of  the  a-  and  co-limit  sets  of  orbits. 

For  fe  X'.let 

La(f  )  =  {/>:  p  £  ct(q)  lor  some  <7} .  )  =  ip'  P  e  o>(d)  for  some  q) . 

Dki  inition  3.5.  Suppose  /€  X' .  A  point  p  6  ft  is  a  wandering  point  of  /  if  there 
are  a  neighborhood  V  of  p  and  t0  >  0  such  that  if  |/|  >  t0.  then  Tj-(t)  V  O  V  =  0.  In  the 
contrary  case,  p  is  a  nonwandering  point  off  The  set  of  nonwandering  points  off  is  de¬ 
noted  by  fi(/). 

In  Definition  3.5,  the  notation  |/|  >  t0  means  for  all  /  >  t0  and  all  r  <  — 10  as  long  as 
the  orbit  is  defined. 

We  remark  that  Sl(f)  D  La(f)  U  Lu(f).  but  it  is  easy  to  construct  examples  for 
which  the  inclusion  is  proper  (see,  for  example.  Palis  and  deMelo  [1]). 

Dki-inition  3.6.  A  vector  /E  XJ,  is  Morse-Smale  (MS)  if  it  is  KS  with  a  finite  num¬ 
ber  of  critical  points  and  periodic  orbits  with  fi(  /  )  equal  to  the  set  of  critical  points  and 
periodic  orbits. 

Some  of  the  basic  results  on  Morse-Smale  systems  are  due  to  Smale  (4 j ,  Palis  (1 )  and 
Palis  and  Smale  [1  j .  They  are  summarized  in  the  following  theorem  which  is  also  valid  for 
vector  fields  on  any  compact  manifold. 

Thiokkm  3.7.  (1 )  The  set  of  MS  systems  is  open  and  nonempty  in  X'  for  any  n. 

(2)  Any  f  &  MS  is  structurally  stable. 

(3)  The  set  of  gradient  vector  fields  which  are  MS  is  open  and  dense  in  the  set  of  all 
gradient  vector  fields. 

Since  the  MS  systems  are  structurally  stable,  they  cannot  be  dense  in  dimension  n  >  3. 
On  the  other  hand,  one  can  ask  if  there  are  any  other  structurally  stable  systems  which  are 
not  MS.  One  way  to  answer  this  question  is  to  construct  a  structurally  stable  system  with 
infinitely  many  periodic  orbits. 

To  see  how  such  a  situation  might  arise,  suppose  SI  C  R3  and  /E  Xr}(Sl)  has  a  hyper¬ 
bolic  periodic  orbit  y.  Let  Ws( y),  Wu(y )  be  the  stable  and  unstable  manifolds  for  y  and  let 
n  be  a  Poincare  map  of  some  transversal  r  of  7  at  p  and  W,(y)  -  W*( 7)  n  r,  W"(y)  —  W"(y)  n  r: 
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that  is,  that  part  of  the  stable  and  unstable  manifolds  in  the  transversal  r.  Then  ^(7), 
IV“(7)  are  *he  local  stable  and  unstable  manifolds  of  the  point  p  as  a  fixed  point  of  the  dif- 
feomorphism  rr.  There  is  the  possibility  that  Wf(y)  £>  W“(y)  contains  points  other  than  the 
fixed  point  p  of  it.  Any  such  point  q  is  called  homoclinic  to  p.  A  point  q  is  called  trans¬ 
verse  homoclinic  to  p  if  W*r(y)  is  transversal  to  W“(y)  at  q.  If  q  is  transverse  homoclinic  to 
p.  then  the  behavior  of  the  stable  and  unstable  manifold  is  very  bad.  In  fact,  since  itWf(y) 
C  ^(7).  *^"(7)  c  and  q  £  Wf( 7)  ft  ^“(7),  q  ^  p,  we  must  have  nnq  £  W’( 7)  O 

W“( 7)  for  all  n  =  0,  ±  1 ,  ±  2,  ...  and  itnq  — ►  p  as  n  — *•  «>.  If,  in  addition,  q  is  transverse 
homoclinic  to  p,  continuity  of  the  map  it  implies  that  the  picture  near  p  must  be  something 
like  the  one  in  Figure  I .  The  arrows  do  not  represent  the  direction  of  a  flow  as  for  vector 
fields,  but  only  that  points  move  in  the  direction  indicated  under  iterates  of  it.  In  Figure  1 , 
we  have  only  indicated  some  of  the  complications  that  are  arising  from  looking  at  the  for¬ 
ward  evolution  of  the  unstable  manifold.  The  same  type  of  thing  must  occur  with  the 
stable  manifold.  Note  that  there  will  be  infinitely  many  transverse  intersections  in  any 
neighborhood  of  the  homoclinic  point  q. 


Figuki  1. 

This  phenomenon  was  observed  by  Poincare  j  1  j .  Birkhoff  1 1 )  proved  that  every 
transverse  homoclinic  point  is  the  limit  of  periodic  points  (that  is.  points  x  such  that  n"x  =  x 
for  some  integer  n)  and  indicated  some  of  the  random  behavior  that  occurs  near  these 
points.  Smale(l)  carried  the  analysis  even  further.  We  briefly  describe  the  results  following 
Moser  [2] .  If  A  is  a  finite  or  countably  infinite  sequence  of  symbols,  let  S  be  the  collec¬ 
tion  of  doubly  infinite  sequences  s  =  {sk ,  k  =  0.  ±  1 ,  . . .}  with  each  sk  £  A.  The  shift 
automorphism  o  on  S  is  defined  by  os  =  (7* ,  k  =  0.  ±  I ,  . . .} ,  -  sk  ,  for  all  k. 

Near  a  homocliuic  point  q,  wc  can  construct  a  small  quadrilateral  Q,  two  of  its  sides 
consisting  of  parts  of  W“(p),  W}(p)  and  the  others  parallel  to  the  tangents  of  these  sets  at 
q  (see  Figure  1 ).  For  any  point  a  £  Q,  let  k  =  k(a)  be  the  smallest  positive  integer  such  that 
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ir*(a)  €  Q,  if  it  exists.  Let  D(Tt )  be  the  set  of  a  £  Q  for  which  such  a  k  exists  and  define 
not  =  rr*(a)  for  a  €Z)(  rr).  The  map  rr  is  called  the  transversal  map  of  rr  for  the  quadrilateral  Q. 

Theorem  3.8.  If  it  is  a  C°° -diffeomorphism  of  the  plane  with  a  point  q  transverse 
homoclinic  to  a  hyperbolic  fixed  point  p,  then  in  a  neighborhood  of  q,  the  transversal  map 
n  of  a  quadrilateral  possesses  an  invariant  set  I  homeomorphic  to  the  sequence  space  S  with 
an  infinite  number  of  symbols  by  a  map  t:  S  — ►  /  such  that  rfr  =  to.  Also,  there  is  an  in¬ 
teger  k,  an  invariant  set  I  of  itk  and  a  homeomorphism  7:  S  — ►  /,  where  S  is  the  sequence 
space  of  a  finite  number  of  symbols,  such  that  itk7  =  7o. 

Note  the  difference  in  the  two  conclusions  in  the  theorem.  In  the  first  part,  the  set  I 
is  invariant  for  rr  and  rr  is  equivalent  on  /  to  the  shift  automorphism  on  an  infinite  number 
of  symbols.  In  the  second  part,  the  set  I  is  invariant  under  a  fixed  power  k  of  it  itself  and 
rr*  is  equivalent  on  /  to  the  shift  automorphism  on  a  finite  number  of  symbols. 

It  follows  immediately  from  Theorem  3.8  that  there  are  infinitely  many  periodic 
points  in  a  neighborhood  of  the  transverse  homoclinic  point  and  they  are  dense.  Also,  there 
is  a  random  behavior  to  the  orbits  on  the  invariant  set  /  (or  /  )  since  knowing  the  early 
terms  of  a  sequence  tells  nothing  about  the  later  terms  of  a  sequence. 

For  examples  of  transverse  homoclinic  points  in  celestial  mechanics,  see  Moser  [2] . 
Transverse  homoclinic  points  also  occur  in  structurally  stable  systems  as  we  shall  see  below. 
More  examples  in  second  order  nonautonomous  differential  equations  will  be  given  later 
when  we  are  studying  analytical  methods  in  bifurcation  theory.  Now,  we  prefer  to  continue 
the  general  survey. 

To  describe  further  aspects  of  the  theory,  it  is  convenient  to  work  with  Diffr(A/), 
r  >  1 ,  the  space  of  diffeomorphisms  with  derivatives  up  through  order  r  on  a  smooth  com¬ 
pact  manifold  M.  This  can  be  related  to  differential  equations  in  several  ways.  One  of  the 
most  important  is  through  a  Poincare  map  for  periodic  orbits  as  described  above.  More  gen¬ 
erally,  if  /£  X'(£l),  M  C  S2  is  compact  and,  for  each  rEM,  there  is  a  t(jc)  >  0  such  that 
Tf{r{x))x  e  M,  Tf{t)x  &  M,  0  <  t  <  t(x),  then  the  map  x  h*  T^(t(x))x  is  in  Diff  r(M)  if  it 
has  the  required  number  of  derivatives. 

If  g  e  Diff  r(M),  a  point  p  £  M  is  a  periodic  point  of  g  if  there  is  a  positive  integer 
n  =  n(p)  such  that  Tnp  =  p.  The  periodic  orbit  is  hyperbolic  if  no  eigenvalue  of  dg(  p)/dx 
has  modulus  one.  For  each  hyperbolic  periodic  point  p  of  g,  one  can  define  the  global 
stable  manifold  Ws(p)  and  unstable  manifold  Wu(p)  in  a  manner  similar  to  the  definitions 
for  vector  fields. 

We  now  give  an  example  due  to  Thom  which  was  an  inspiration  for  many  further  de¬ 
velopments  in  dynamical  systems.  In  R2  identify  the  points  (x,  y),  (x  +  m,  y  +  n)  for  all 
integers  m,  n.  Any  unit  square  with  integer  vertices  may  be  identified  with  the  torus  T2 
and  any  mapping  of  the  plane  into  itself  yields  a  mapping  of  T2  into  T2  in  the  obvious  way. 
Let  L  be  a  2  x  2  matrix  with  integer  coefficients,  determinant  1  and  real  eigenvalues.  The 
eigenvalues  are  then  X,  X-1  with  X  <  1  irrational.  This  implies  that  the  linear  subspaces 
E’,  Eu  generated  respectively  by  the  eigenvectors  for  X,  X-1  have  irrational  slope.  For  any 
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x  6  R2 ,  each  of  the  lines  x  +  E1,  x  +  Eu  is  invariant  under  L.  The  map  L  on  R2  generates 
a  natural  map  tt  on  T2  obtained  from  the  above  identification  of  T2  with  unit  squares  in 
R2  with  integer  coefficients.  If  p  =  nx  £  T2.  x  G  R2,  let  W*(p)  =  n{x  +  £*),  H'u(  p)  = 
rr(x  +  Eu).  If  p  is  a  periodic  point  of  tt,  then  W*(p),  Wu(p)  are  respectively  the  stable  and 
unstable  manifolds  of  p.  Since  the  slopes  of  the  linear  subspaces  E*,  Eu  are  irrational,  the 
sets  W*(p)  and  Wu(p)  are  dense  in  T2  for  every  p&T2 .  Furthermore,  it  is  not  difficult  to 
show  that  every  point  of  intersection  of  these  sets  is  a  point  of  transversal  intersection.  Also, 
the  points  of  intersection  are  dense  in  T2 .  In  particular,  there  is  a  dense  set  of  points  trans¬ 
verse  homoclinic  to  the  critical  point  p  =  ir(0).  Since  each  transverse  homoclinic  point  is 
the  limit  of  periodic  points,  it  follows  that  the  periodic  points  of  n  are  dense  in  T2 .  A  sim¬ 
ple  direct  proof  of  this  last  result  is  contained  in  Palis  and  deMelo  (1,  p.  171]. 

With  it:  T2  — *  T 2  defined  as  above,  Anosov  (I],  [2]  showed  that  it  is  structurally 
stable.  A  more  elementary  proof  was  given  by  Moser  [1],  [2].  Since  n  is  structurally  stable 
and  contains  infinitely  many  periodic  orbits,  this  necessarily  implies  that  the  Morse-Smale 
systems  are  not  dense  in  the  set  of  structurally  stable  systems. 

The  above  example  was  generalized  by  Anosov  [2]  in  the  following  wgy. 

Definition  3.9.  Let  Af  be  a  compact  manifold.  An  f  £  Diffr(Af),  r  >  1 .  is  an  Ano¬ 
sov  diffeomorphism  if  the  tangent  space  at  each  point  x  of  M  is  a  direct  sum  Ex  ©  £“  in¬ 
variant  under  the  derivative  Df;  that  is,  DfxEx  =  DfxEx  =  and  there  is  a  Rie- 

mannian  metric  on  M  and  a  constant  X  €  (0,  1)  such  that  \Dfxv\  <  X|u|,  \Df~'u\  <  XI  u|  for 
allxGM,  uGEj.u  6  E“. 

Anosov  [2]  has  shown  that  these  diffeomorphisms  are  structurally  stable.  A  simpler 
proof  was  given  by  Moser  [I] ,  [2] .  For  a  discussion  of  the  restrictions  that  are  imposed  on 
the  manifold  M  in  order  for  it  to  admit  an  Anosov  diffeomorphism,  see  Palis  and  deMelo  |1] . 

The  next  important  step  in  the  abstract  theory  of  dynamical  systems  was  taken  by 
Smale  [2]  by  defining  systems  which  satisfy  Axiom  A.  Suppose  /  G  Diffr(Af)  and  AC  M 
is  a  closed  invariant  set.  The  set  A  is  said  to  have  a  hyperbolic  structure  if  the  tangent  space 
at  each  point  x  G  A  is  the  direct  sum  £*  ©  E“  invariant  under  Df  and  there  are  a  Rieman- 
nian  metric  and  X  G  (0,  1)  such  that  \Dfxv\  <  X|u|,  \Df~xu\  <  X|«|  for  x  G  A,  u  G  Ex. 
u  G  £x.  If  A  is  hyperbolic,  it  is  possible  to  define  stable  and  unstable  manifolds  for  the 
set  A  by  looking  at  asymptotic  orbits. 

Definition  3.10.  / 6  Diff(Af)  satisfies  Axiom  A  if  the  set  of  nonwandering  points 
Q, (/)  is  hyperbolic  and  the  periodic  point,  of  f  are  dense  in  f 1(f). 

If  /  satisfies  Axiom  A,  Smale  [2]  has  shown  that  J2(/)  =  fl,  U  fl2  U  ■  •  •  U  £lk  where 
each  is  closed  invariant  and  transitive;  that  is,  has  a  dense  orbit .  Robinson  ( 1 1  has  shown 
that  any  /G  Diffr(Af ),  r>  1 ,  is  structurally  stable  if  it  satisfies  Axiom  A  and  all  stable  and 
unstable  manifolds  intersect  transversally.  An  /G  Diffr(Af)  is  said  to  be  absolutely  stable  if 
there  are  a  neighborhood  V(f)  C  Diffr(Af)  of  / and  a  constant  K  >  0  such  that,  for  every 
g  6  V(f),  there  is  a  homeomorphism  h  of  M  such  that  hf  =  gh  and  | Ai  -  /|0  <  K\f  -  g|0 
where  |-|0  designates  the  norm  in  C°.  Results  of  Franks  [1  ] .  Guckenheimer  [1]  and 
Mafld  (1)  show  that  /  is  absolutely  stable  if  and  only  if  it  satisfies  Axiom  A  and  all  stable 


12 


JACK  K.  HALK 


and  unstable  manifolds  intersect  transversally.  Marie  (unpublished)  has  also  recently  shown 
that  Axiom  A  is  implied  by  structural  stability  and  a  technical  condition  on  the  characteris¬ 
tic  exponents  of  Liapunov  on  (2(/). 

It  is  also  possible  to  study  structural  stability  restricted  to  the  set  of  nonwandering 
points.  More  specifically,  /€  Diffr(M)  is  said  to  be  (2- stable  if  there  is  a  neighborhood  V(f) 
of  /such  that,  for  every  g  €  V{f),  there  is  a  homeomorphism  h:  (2(/)  — *■  Mg)  such  that 
hf  -  gh  on  (2(/).  If  /satisfies  Axiom  A,  Mf)  —  S2,  U  •  •  •  u  Slk,  then  a  cycle  of  ft  is  a 
sequence  p,  e  (2k) ,  . . . ,  ps  €  (2^  =  (2^  such  that  Ws(.Pj)  n  Wu(pl4  , )  #  0,  1  <  i  <s  -  1 . 
Smale  [2]  showed  that  Axiom  A  and  no  cycles  imply  (2-stable.  Palis  (2)  has  shown  that 
any  /  satisfying  Axiom  A  is  not  (2-stable  if  it  has  a  cycle.  It  is  not  known  if  /  (2-stable  im¬ 
plies  that  it  must  satisfy  Axiom  A. 


4.  Stability  and  bifurcation  at  a  zero  eigenvalue 

Suppose  A  is  an  n  x  n  constant  matrix  whose  eigenvalues  have  negative  real  parts, 

.x  €  R,  v  £  R".  X0  €  C*( R  x  R" ,  R),  F0  6  Ck(R  x  R",  Rn),  k  >  1,  are  functions  of 
(x,  y)  £  R  x  R"  which  vanish  together  with  their  first  derivatives  at  x  =  0,  y  =  0.  One  of 
our  objectives  in  this  section  is  to  determine  how  the  stability  properties  of  the  solution 
x  =  0,  v  =  0,  of  the  equation 

(4.1)  x  =  X0(x,  j>),  y  =  Ay  +  r0(x,  v) 

depend  upon  the  nonlinear  functions  X0,  Y0. 

For  a  given  neighborhood  U  of  x  =  0,  y  =  0,  and  X  £  Ck(U,  R),  Y  £  Ck(U,  R), 
k>  1,  we  also  discuss  how  the  number  of  equilibrium  solutions  of  the  perturbed  equation 

(4.2)  x  =  X0(x.  y)  +  X(x,  y),  y  =  Ay  +  F0(x,  y)  +  Y(x,  y) 

depends  upon  X,  Y  in  a  neighborhood  of  X  -  0,  Y  —  0.  The  stability  properties  of  these 
solutions  will  be  considered.  A  complete  classification  is  given  for  some  cases.  Die  neigh¬ 
borhood  U  will  be  fixed  and  our  analysis  always  will  be  in  some  other  neighborhood  V  of 
x  =  0,  y  =  0  which  will  be  a  subset  of  U. 

Lkmma  4.1 .  There  are  a  neighborhood  V  of  (x,  y)  =  (0,  0)  and  a  neighborhood  W 
of  ( X ,  Y)  =  (0,  0)  such  that,  for  ( X ,  Y)  £  W  and  any  equilibrium  point  a  =  ot.{X,  Y)  of 

(4.2 )  in  V,  the  stable  manifold  has  dimension  >  n.  Furthermore,  there  is  an  n-dimensional 
manifold  S( a,  X,  Y)  in  V  containing  a  such  that  any  solution  of  (4.2)  with  initial  value  in 
S(a,  X,  Y )  approaches  a  as  t  — »  00  exponentially.  Furthermore,  S{ a,  X,  Y)  is  Ck  in  (X,  K) 
and  the  tangent  space  to  S{ a,  X,  V)  at  a  approaches  the  space  {(0,  y).  y  £  R" }  as  (.V,  )' ) 
-»  (0.  0). 

A  proot  will  not  be  given.  It  is  technical  but  may  be  supplied  by  using  the  ideas  in 
Coddington  and  Levinson  (1),  Hale  [1)  or  Hartman  [2]  for  the  construction  of  stable  man¬ 
ifolds. 

For  equation  (4.2),  a  set  SC  Rn+i  is  a  local  invariant  manifold  if,  for  any  (x0,  y0) 

£  S,  there  is  a  T  >  0  such  that  the  solution  of  (4.2)  through  (x0,  y0)  remains  in  S  for 
|/i  <  T.  If  T  =  °°,  we  say  S  is  an  invariant  manifold.  If  there  is  a  C1  function  ^(x,  X,  )') 
defined  in  a  neighborhood  of  (x,  X,  Y)  =  (0,  0,  0)  such  that  ii/(0,  0,  0)  =  0.  9i^(0,0.  0)/3x 
=  0  and  S  =  S{X,  K)  =  {(x,  y):  y  =  i l/(x,  X,  F)}  is  a  local  invariant  manifold  for  (4.2). 
then  we  say  S  is  a  local  center  manifold  for  (4.2). 
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Because  the  matrix  A  has  eigenvalues  with  negative  real  parts,  it  is  reasonable  to  ex¬ 
pect  that  there  is  a  local  center  manifold  for  (4.2).  This  idea  is  implicit  in  tire  papers  of 
Krylov  and  Bogoliubov  written  in  the  1930’s  (see  Hale  (I)  for  references)  and  was  used  ex¬ 
plicitly  by  V.  Pliss  [1]  to  study  stability.  Kelley  [1] ,  [2]  has  given  a  very  complete  and 
readable  account  of  this  idea.  More  precisely,  he  has  shown  the  existence  of  a  center  mani¬ 
fold  which  is  C*  in  (x,  X,  T)  in  a  sufficiently  small  neighborhood  V  (depending  on  k)  of 
(x,  y)  =  (0,  0).  A  center  manifold  is  also  exponentially  stable  in  the  sense  that  any  solution 
which  remains  in  the  neighborhood  V  of  (x,  y)  =  (0,  0)  approaches  the  center  manifold  ex¬ 
ponentially  as  t  — *  oo. 

The  flow  on  a  center  manifold  S(X,  K)  is  given  by  the  scalar  equation 
x  =  u(x,  X,  Y ). 

(4.3) 

t/(x.  X,  Y )  =  X0(x.  0(x.  X,  Y))  +  X(x,  0(x,  X,  Y)). 

If  a  is  an  equilibrium  point  of  equation  (4.2),  then  o  e  S(X,  T),  a  =  (x0,  0(xo,  X,  Y))  for 
some  x0.  The  converse  is  also  true.  The  stability  properties  of  x0  as  a  solution  of  (4.3)  are 
easily  determined  from  the  behavior  of  u(x,  X,  >')  near  x  =  x0.  From  Lemma  4.1 ,  the 
stability  properties  of  x0  are  reflected  in  an  obvious  way  in  the  stability  properties  of  a  as 
a  solution  of  (4.2).  In  particular,  x0  being  a  stable  (asymptotically  stable)  (unstable)  equilib¬ 
rium  point  of  (4.3)  is  equivalent  to  a  being  a  stable  (asymptotically  stable)  (unstable)  equi¬ 
librium  point  of  (4.2). 

From  the  above  discussion,  the  complete  structure  of  the  flow  for  (4.2)  is  deter¬ 
mined  by  knowing  the  flow  for  (4.3).  Thus,  if  we  can  determine  a  center  manifold,  the 
problem  will  be  solved.  Generally,  this  is  extremely  difficult  to  do.  We  give  an  alterna¬ 
tive  procedure  for  determining  the  qualitative  properties  of  the  solutions  of  (4.2)  which 
is  very  close  to  the  original  method  of  Liapunov  [1],  It  uses  only  the  bifurcation  func¬ 
tion  obtained  by  an  application  of  the  method  of  Liapunov-Schmidt  for  the  equilibrium 
solutions  of  equation  (4.2). 

This  method  has  other  advantages  over  the  use  of  the  center  manifold.  The  bifurca¬ 
tion  function  has  the  same  smoothness  properties  as  the  original  vector  field.  This  is  not 
true  of  the  local  center  manifold.  It  inherits  the  C*  property  for  any  finite  k,  but  may  not 
be  C"  (or  analytic  even)  if  the  original  vector  fields  are  C“  (or  analytic).  For  examples, 
see  Carr  ( 1 1 . 

There  is  a  neighborhood  of  (0,  0)£R  «  R"  and  a  neighborhood  of  Y  =  0  such  that 
the  equation 

(4.4)  A<t>  +  Y0(x,  0 )  +  K(x,  0)  =  0 

has  a  unique  solution  0(x,  Y)  with  0(0.  0)  =  0.  This  function  is  Ck  in  (x,  K).  Let 

(4.5)  G(x.  X,  Y)  =  *0(x.  0(x,  K))  +  X(x,  0(x,  Y)). 

We  can  now  prove  the  following  result. 
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Theorem  4.2.  There  are  neighborhoods  F,  of  x  =  0,  F2  of  y  =  0,  IF,  of  X  -  0, 

W2  of  Y  =  0  anr/  a  Ck -function  tp:  F,  x  IF2  — ►  F2  arc/t  Zhar  the  equilibrium  points  (x0.  yQ ) 
of  (4.2)  in  F,  x  F2  /or  (X,  K)  €  IF,  x  1F2  are  g-rVezr  by 

(4.6)  G(x0,  X,  F)  =  0,  >0  =  4Hx0.  Y) 

where  G  is  defined  in  (4.5).  Furthermore,  the  stability  properties  of  the  solution  x0  of  the 
scalar  equation 

(4.7)  i  =  G(x,  X,  Y) 

are  the  same  as  the  stability  properties  of  x0  on  a  center  manifold.  Finally,  any  solution  of 
(4.2)  which  remains  in  the  neighborhood  F,  x  F2  for  t  3=  0  approaches  an  equilibrium 
point  of  (4.2). 

The  function  G(x,  X,  T)  is  easier  to  compute  than  the  vector  field  u(x,  X,  Y)  on  a 
center  manifold.  The  importance  of  the  result  is  that  the  function  G(x,  X,  Y),  called  the 
bifurcation  function,  carries  all  of  the  information  necessary  to  determine  the  stability 
properties  of  the  equilibrium  points  on  the  center  manifold.  The  conclusion  on  stability  in 
Theorem  4.2  is  due  to  deOliveira  and  Hale  [1] . 

Proof  of  Theorem  4.2.  It  is  clear  that  the  equilibrium  points  of  (4.2)  are  given  by 
(4.6).  If  u(x,  X,  Y)  is  the  vector  field  on  a  center  manifold,  then  u(x,  X,  Y)  and  G(x,  X,  Y) 
have  the  same  set  of  zeros  in  a  neighborhood  of  zero.  If  x0  is  an  isolated  zero,  let  us  first 
prove  that  G,  u  have  the  same  sign  in  a  neighborhood  of  x0.  Without  loss  in  generality,  we 
can  suppose  x0  =  0  and  it  is  simple  because  we  can  make  a  small  perturbation  in  X,  Y  to 
attain  this  property.  So,  assume  G( 0,  X,  Y)  =  0  =  u(0,  X,  Y),  9G(0,  X,  Y)/dx  >  0, 

3u(0,  X,  Y)jhx  <  0.  Now,  consider  the  functions  G(x,  X  +  e,  T),  u(x.  X  +  e,  Y)  for  e  a 
small  real  number.  Then  one  can  show  that 

G{ x,  X  +  e,  F)  =  G(x,  X,  Y)  +  e, 

u{ x,  X  +  e,  Y)  =  u(x,  X,  Y)  +  6(x,  y)e  +  o(|el) 

as  |e|  — *  0  where  6(x,  X,  Y)  >  0  for  x  in  a  neighborhood  of  zero.  The  assertion  about  G 
is  obvious,  but  the  one  about  u  is  nontrivial  and  uses  the  properties  of  a  center  manifold. 
Since  the  derivatives  of  G  and  u  with  respect  to  x  at  (0,  0,  0)  have  opposite  sign,  this  shows 

that  a  small  perturbation  of  e  gives  rise  to  distinct  zeros  of  G,  u  in  a  neighborhood  of  x  =  0. 

However,  this  is  a  contradiction  since  G,  u  must  always  have  the  same  set  of  zeros.  Thus,  if 
x0  is  an  isolated  zero  of  G,  u,  then  the  functions  G,  u  have  the  same  sign  in  a  neighborhood 
of  zero.  Thus,  the  stability  properties  as  a  solution  of  (4.7)  and  (4.3)  are  the  same.  If  x0 
is  not  isolated,  then  one  can  easily  discuss  all  the  ways  in  which  this  can  occur  to  complete 
the  discussion  of  the  stability  of  x0. 

If  a  positive  orbit  stays  in  F,  x  F2,  then  its  w-limit  set  must  be  on  the  center  mani¬ 
fold  and  it  must  be  a  connected  invariant  set.  Also.  Lemma  4.1  implies  that  each  equilib¬ 
rium  point  has  a  stable  manifold  of  codimension  one  transverse  to  the  center  manifold.  If 
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the  cu-limit  set  contains  two  distinct  points,  then  it  must  contain  the  arc  between  these 
two  points  on  the  center  manifold.  If  there  are  no  equilibrium  points  between  these  two, 
then  the  positive  orbit  must  remain  between  the  transverse  stable  manifolds  above  and  go 
back  and  forth  along  the  arc  as  t  — *■  °°,  an  obvious  contradiction.  The  same  type  of 
argument  can  be  used  even  if  there  are  any  number  of  points  in  the  arc.  The  theorem 
is  proved. 

The  proof  of  the  last  statement  in  the  theorem  can  be  modified  following  Hale  and 
Massatt  [1]  to  obtain  the  following  result. 

Theorem  4.3.  The  co-limit  set  of  a  bounded  orbit  of  the  gradient  system  x  =  fix), 
fix)  =  grad  F\x)  is  a  single  point  if.  for  every  zero  x0  of  f  which  belongs  to  a  continuum 
of  zeros  of  f  the  matrix  df(x0)/dx  is  either  nonsingular  or  has  zero  as  a  simple  eigenvalue. 

Proof.  We  only  indicate  a  proof.  Details  may  be  found  in  Hale  and  Massatt  [1], 
Differentiating  F\x)  along  the  solutions  one  sees  that  the  oj-limit  set  of  every  solution  must 
belong  to  the  set  of  equilibrium  points.  Since  the  co-limit  set  is  connected,  it  must  contain 
either  a  single  equilibrium  point  or  a  continuum  of  such  points.  If  x0  belongs  to  a  contin¬ 
uum  of  the  co-limit  set,  then  the  hypothesis  on  3/(x0)/9x  implies  the  existence  of  a  smooth 
center  manifold  through  x0  of  dimension  one.  Thus,  the  co-limit  set  near  x0  is  a  smooth 
arc.  Since  it  is  compact,  it  must  be  an  arc  or  a  closed  curve.  But  then  one  can  use  a  gen¬ 
eralization  of  lemma  4.1  to  show  there  is  a  tubular  neighborhood  of  the  co-limit  set  such 
that  any  orbit  which  remains  in  this  neighborhood  for  all  r  >  0  lies  on  a  stable  manifold 
for  an  equilibrium  point.  This  proves  the  result. 

As  a  consequence  of  Theorem  4.2,  we  have 

Corollary  4.4.  An  equilibrium  solution  (x0,  y0)  e  K,  x  V2  of  (4.2)  for 
iX,  Y)  €  W,  x  Wj  is  asymptotically  stable  if  and  only  if  there  is  an  e  >  0  such  that 
G(x0  +  u.  X,  Y)u  <  0  for  0  <  |u|  <  e.  It  is  unstable  if  and  only  if  there  is  an  e  >  0  such 
that  G(x0  +  u.  X,  Y)u  >  0  for  either  0<u<e  or  ~e  <u  <  0 .  If  there  is  an  e>  0  such 
that  G(x0  +  u.  X,  Y)  =  0  for  0  <  |w|  <  e,  then  the  solution  is  stable  and  there  is  a  Ck  first 
integral  in  a  neighborhood  of  (x0,  y0).  If  the  functions  A^,  X,  Y0,  Y  are  analytic  (or  C°°) 
then  this  first  integral  is  analytic  (or  C°°). 

Proof.  Everything  is  obvious  except  the  existence  of  the  first  integral.  Lemma  4.1 
implies  the  existence  of  an  » -dimensional  Ck  -stable  manifold  above  each  equilibrium  point 
x0  +  u,  |u|  <  e.  These  manifolds  are  analytic  if  the  functions  in  (4.2)  are  analytic.  Using 
these  stable  manifolds,  one  can  make  a  C*  (or  analytic)  change  of  variables  (x,  y)  t-»  (v,  w) 
so  that  equation  (4.2)  in  a  neighborhood  of  (x0,  y0)  becomes  ii  =  0,  w  =  Aw  +  W(u,  w). 
This  equation  has  a  first  integral  V(v,  w)  =  u.  Thus,  the  original  equation  has  a  first  integral 
and  the  proof  is  complete. 

Corollary  4.5.  Suppose  G(x0,  X,  Y)  =  0  and  y0  =  0(xo,  Y).  <t>  satisfying  (4.4). 

If  there  is  an  integer  q  >  1  and  a  (3  +  0  such  that 

(4.8)  G(x0  +  u.  X,  y)  =  0uq  +  «(|u|*)  as  |u|  0, 
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then  the  solution  (x0,  y0)  of  (4.2)  is  asymptotically  stable  if  and  only  if  0  <0,q  odd. 
Otherwise,  it  is  unstable. 

We  have  also  the  following  results  of  Liapunov  [1]  (see,  also,  Bibikov  [1]). 

Corollary  4.6.  Suppose  G(x0,  X,  Y)  =  0 ,y0  =  <p(x0,  Y),  <t>  satisfying  (4.4)  and 
A0,  X,  Y0,  Y  are  analytic.  Then  either  there  are  an  integer  q  >  1  and  a  0  d*  0  such  that 

(4.8)  is  satisfied  with  the  stability  properties  as  stated  in  Corollary  4.5  or  (x0,  y0)  is  a  stable 
solution  of  (4.2)  and  there  is  an  analytic  first  integral 

Proof.  Since  G(x,  X,  Y)  is  analytic,  only  the  two  cases  mentioned  can  occur. 

Corollary  4.7.  If  there  is  a  scalar  function  H(x,  y,  z)  continuous  for  (x,  y,  z)  £ 

R  «  R"  >  R"  such  that  H(x,  y,  0)  =  0  and 

X0(x,  y )  +  X(x,  y)  =  H(x,  y.  Ay  +  Y0( x,  y)  +  Y( x,  y)) 

then  the  zero  solution  of  (4.2)  is  stable  and  there  is  a  first  integral. 

Proof.  The  hypothesis  implies  G(x,  X,  Y)  =  0  for  x  in  a  neighborhood  of  X  =  0. 

Let  us  now  turn  to  the  problem  of  bifurcation  of  the  equilibrium  point  (0,  0)  treating 
X,  Y  in  (4.2)  as  parameters.  An  equilibrium  point  (x0,  y0)  of  (4.2)  is  called  a  saddle-node 
if  xQ  as  a  solution  of  (4.7)  or  (4.3)  is  asymptotically  stable  from  one  side  and  unstable  from 
the  other.  The  following  result  is  contained  in  Andronov  et  al.  [1] ,  Sotomayor  [  1  ] . 

Theorem  4.8.  Suppose  k  >  2  and  there  is  a  (id*  0  such  that 

(4.9)  G(x,  0,  0)  =  fix2  +  o(|*i2)  as  Ixl  —  0. 

Then  there  is  a  neighborhood  V  of  (x,  y)  =  (0,  0),  a  neighborhood  W  of(X,  T)  =  0  and  a 
Ck -function  y:  W  — *■  R  such  that  the  following  conclusions  hold: 

(i)  y(X,  Y)  >  0  implies  there  are  no  equilibrium  points  of  (4.2)  in  V. 

(ii)  y(X,  Y)  =  0  implies  there  is  one  equilibrium  point  of  (4.2)  in  V  which  is  a  saddle- 

node. 

(iii)  y(X,  T)  <  0  implies  there  are  two  equilibrium  points  of  (4.2)  in  V,  one  is  a  hy¬ 
perbolic  saddle  and  the  other  is  hyperbolic  and  asymptotically  stable. 

Furthermore,  the  set  T  =  {(A,  Y):  y(X,  Y)  =  0}  is  a  submanifold  of  codimension  1 
in  W  (see  Figure  2). 

i  roof.  The  function  G(x,  A,  T)  satisfies  8G(0,  0,  0)/8x  =  0,  d2G(0,  0,  0)/dx2  = 

2/3  d  0.  Thus,  there  are  a  neighborhood  W  of  (A,  Y)  =  (0,  0)  and  a  neighborhood  V ,  of 
x  =  0  such  that  the  equation  dG(x,  A,  Y)/bx  =  0  has  a  unique  solution  x*(A,  T)  with 
x*(0,  0)  =  0.  If  y(X,  Y)  =  (sgn  0)G(x*(X,  T),  A,  Y),  then  the  conclusion  of  the  first  part 
of  the  theorem  follows  from  Theorem  4.2. 

To  prove  the  last  part  of  the  theorem,  consider  the  special  case  where  A(x.  v)  = 

Y(x,  y)  =  0  for  all  x,  y  where  X  is  a  real  scalar  and  show  that  Dxy(\,  0)  d  0  for  X  =  0. 
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G(x,  0,  0)  =  0x3  +  o(lxl3). 


Figure  2. 

G(x,  0,  0)  =  0xJ  +  oflx|2). 

To  discuss  a  higher  order  bifurcation,  suppose  that  k  >  3  and 

(4.10)  G(x.  0,  0)  =  fix3  +  o(|x|3)  as  |x|  —  0,  0  *  0. 

The  argument  below  follows  Chow,  Hale  and  Mallet-Paret  [1]  and  Vanderbauwhede  [1] 
Then  there  are  a  neighborhood  F,  of  x  =  0  and  a  neighborhood  W  of  ( X ,  Y )  =  (0,  0)  such 
that  the  equation  3 2G(x,  A,  Y)/dx2  =  0  has  a  unique  solution  x*(A,  Y)  with  x*(0.  0)  =  0 
Let 


(4-11) 


?,(*.  n  = 


3 G(x*(X,  Y),  X,  Y) 


dx 


If  we  apply  the  same  argument  as  in  the  proof  of  Theorem  4.8,  then  for  (X,  Y)  €  W.  the 
equation  3 G(x.  X,  Y)ldx  has  no  solution  in  F,  if  07, (A-,  T)  >  0,  one  solution  if  7, (A,  K) 
=  0  and  two  simple  solutions  if  07,  (X,  Y)  <  0. 

If  07,  (X.  Y)  >  0,  then  G(x,  X,  T)  is  strictly  increasing  and  has  exactly  one  solution 
in  V, .  If  7,  (AT,  Y)  =  0,  let  x*(X.  Y)  be  the  unique  solution  of  3C(x,  X,  Y)ldx  =  0  in  F, 
and  define 

(4.12)  72(A’,  T)  =  G(x*(X,  K),  A,  Y). 

If  7, (A,  T)  =  0,  7j(A,  T)  =  0.  then  the  unique  solution  of  G(x,  A,  Y)  =  0  in  F,  is  a  tri¬ 
ple  zero.  If  either  07, (A,  Y)  >  0  or  7,(A,  T)  =  0,  72(A,  T)  ¥=  0,  then  this  unique  solution 
is  simple. 

Now  suppose  07, (A,  T)  <  0.  Then  there  are  two  simple  solutions  of  3 G(x,  A,  Y)/dx 
=  0  in  F, .  Again  using  the  argument  of  the  proof  of  Theorem  4.8.  one  can  show  that  they 
have  the  form 

x?  =  x*{X,  Y)  +  o *(A,  Y.  7, (A,  T))(-07,(A.  Y))'12 

where  o*(A,  Y,  0)  =  ±3~112  (see  Vanderbauwhede  [IJ  for  details). 

If 

(4.13)  7^,  Y)  =  (?<*;( A.  Y).  A,  T)C(x.*(A,  Y),  A.  V) 
then,  using  Theorem  4.2,  we  can  state  the  following  theorem. 
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Thforem  4.9.  Suppose  G(x,  0,  0)  satisfies  (4.10),  the  regions  Vt,  W  and  functions 
7 1 , 7  on  W  are  defined  as  above.  Then  the  following  conclusions  hold: 

(i)  If  either  (}y  , [X,  Y)  >  0  or  py ,  (X,  Y)  <  0  and  y(X,  T)  >  0,  then  equation  (4.6) 
has  one  zero  in  V  and  equation  (4.2)  has  exactly  one  equilibrium  point  in  a  neighborhood 
of  zero  which  is  asymptotically  stable  if  l 3  <  0  and  unstable  if  $  >  0. 

(ii)  // Pyl(X,  Y)  <  0  and  y(X,  Y)  =  0,  then  equation  (4.6)  has  two  zeros  and  equa¬ 
tion  (4.2)  has  two  equilibrium  points  in  a  neighborhood  of  zero,  one  being  a  saddle-node 
and  the  other  hyperbolic.  The  hyperbolic  one  is  asymptotically  stable  if  p  <  0  and  a  saddle 
point  of  order  one  if  (3  >  0. 

(iii)  If  Pyt(X,  Y)  <  0  and  y(X,  Y)  <  0,  there  are  three  simple  solutions  of  equation 
(4.6)  and  equation  (4.2)  has  three  equilibrium  points  in  a  neighborhood  of  zero,  all  hyper¬ 
bolic,  two  saddles  of  order  one  if  i 3  >  0  and  one  saddle  of  order  one  if  (3  <  0 ,  the  others 
being  asymptotically  stable. 

In  the  above  theorem,  the  set  T  =  {(A",  K):  pyx(X,  Y)  <  0,  7(Af,  K)  =  0}  is  the  bi¬ 
furcation  set  in  that  is,  the  set  where  the  number  of  equilibrium  points  of  (4.2)  changes 
from  one  to  three  and  where  the  topological  structure  of  the  trajectories  of  (4.2)  changes  in 
a  neighborhood  of  zero. 

It  is  possible  to  show  that  there  are  positive  nonzero  functions  o+ ,  o_  on  W  such  that 
saying  that  (X,  f)6F  implies  either 

72(AT,  K)  =  ct+(AT,  r)(-(37,()(,  Y))3'3 
or 

7j(*.  Y)  =  -oJX,  YX-py^X,  y))3/2 

where  y2  is  defined  in  (4.12).  This  shows  that  the  set  T  is  like  a  cusp  surface  in  W  as  in 
Figure  3  (see  Vanderbauwhede  [1]  for  details). 

With  G  satisfying  (4.9),  Theorem  4.8  implies  the  complete  behavior  of  the  solutions 
of  (4.2)  in  a  neighborhood  of  zero  is  characterized  by  one  function  of  X,  Y.  If  G  satisfies 
(4.10),  then  two  functions  were  needed  to  obtain  this  characterization.  If  the  degeneracy  of 
G(x,  X,  Y)  at  ( X ,  Y)  =  (0,  0)  is  higher  order,  then  it  becomes  extremely  difficult  to  give  an 
explicit  description  of  how  the  structure  of  the  zeros  depends  on  X,  Y.  This  is  due  to  the 
fact  that  it  is  always  difficult  to  discuss  the  zeros  of  polynomials  of  degree  greater  than  three. 

Singularity  theory  can  be  used  to  say  that  the  problem  actually  is  reducible  to  the  dis¬ 
cussion  of  a  polynomial.  If  we  suppose  all  functions  in  (4.2)  are  C",  then  G(x,  X,  Y)  is 
C~  in  x.  *,  Y.  If 

(4.14)  G(x,  0,  0)  =  (3xq  +  o(\x\q  +  1 )  as  !x|  —  0,  j3  *  0. 

then  the  Banach  space  version  of  the  Malgrange  Preparation  Theorem  due  to  Michor  [1|  im¬ 
plies  there  are  a  polynomial  P(x,  X,  T),  P{x,  0,  0)  =  (3x?  of  degree  q  in  x,  and  a  positive 
function  F.(x,  X,  T),  E(x,  0,  0)  =  1,  such  that 

(4.15)  G(x,  X,  Y)  =  E(x,  X,  Y)l\x,  X.  Y) 
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and  these  functions  are  C"  in  x,  X,  Y.  Thus  the  zeros  of  G(x,  X,  Y)  as  well  as  the  signs 
between  zeros  coincide  with  those  of  P(x,  X,  Y).  Thus,  we  can  state  an  analogue  of  Theo¬ 
rem  4.2  in  the  following  way.  The  equilibrium  points  of  (4.2)  are  given  by 

(4.16)  I\x0.  X,  Y)  =  0,  y  =  0(xo,  Y) 

and  the  stability  properties  of  the  equilibrium  points  are  determined  by  the  stability  proper¬ 
ties  of  the  equilibrium  point  xg  of  the  scalar  equation 

(4.17)  x  =  /\x,X,Y). 

Such  a  result  could  not  be  obtained  by  reducing  the  discussion  of  (4.2)  to  the  center 
manifold  (4.3)  and  then  applying  the  Malgrange  theorem.  The  reason  is  that  the  center 
manifold  is  not  necessarily  C“  (see  Carr  [  1 } ). 

For  the  special  case  q  -  2,  one  can  take  P{x,  X,  Y)  =  fix2  +  y(X,  Y)  where  y  is  the 
function  in  Theorem  4.9.  For  q  —  3,  one  can  take  P(x,  X,  Y)  =  fix3  +  yt(X,  Y)x  + 
y3(X,  T)  where  y3  is  essentially  the  same  function  as  y2  in  (4.12). 

In  the  applications,  one  frequently  encounters  the  situation  where  the  perturbed  vec¬ 
tor  field  ( X ,  Y)  is  not  an  arbitrary  function  in  a  neighborhood  of  zero,  but  belongs  to  a 
prescribed  k  parameter  family  of  vector  fields  X(x,  y,  \),  Y{x,  y,  X),  X  £  R*,  which  vanish 
for  X  =  0.  The  bifurcation  function  is  then  a  function  of  X,  G  =  G(x,  X).  The  number  of 
equilibrium  points  and  their  stability  properties  are  determined  by  the  manner  in  wliich  the 
k  parameter  family  of  vector  fields  crosses  the  bifurcation  surfaces  obtained  from  taking 
arbitrary  perturbations  X,  Y.  Of  course,  it  is  clear  that  the  same  conclusions  will  prevail  if 
the  original  problem  is  phrased  in  terms  of  the  restricted  family  depending  on  X.  This  latter 
method  is  often  more  transparent. 

As  a  final  remark,  it  is  possible  to  extend  everything  above  to  the  case  where  the 
eigenvalues  of  A  in  equation  (4.1)  have  nonzero  real  parts. 
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5.  Stability  and  bifurcation  from  a  focus 

In  this  section,  we  discuss  bifurcation  from  an  equilibrium  point  in  the  spirit  of  the 
previous  section  except  under  the  hypothesis  that  the  linear  approximation  of  the  vector 
field  has  two  purely  imaginary  roots  and  the  remaining  ones  have  negative  real  parts. 

The  unperturbed  equation  is  given  as 

(5.1)  x  =  Bx  +  X0(x,  y ),  y  =  Ay  +  K0(x,  y) 

where 


A  is  an  n  x  n  matrix  whose  eigenvalues  have  negative  real  parts,  X0  G  C*(R2  x  Rn,  R2), 
y0  €  Ck( R2  x  R",  R"),  k  >  1,  are  functions  vanishing  together  with  their  first  derivatives 
at  x  -  0,  y  =  0. 

For  a  given  neighborhood  f/  of  x  =  0,  y  =  0  and  X  G  Ck{U,  R2),  Y  G  Ck(U,  R"), 
the  perturbed  equation  is 

(5.2)  X  =  Bx  +  2f0(x, y)  +  X(x,  y),  y  =  Ay  +  K„(x,  y)  +  Y(x,  y). 

For  (X,  Y)  in  a  sufficiently  small  neighborhood  of  (0,  0),  equation  (5.2)  has  a  unique 
equilibrium  point  x  =  a(X.  Y),  y  =  (i(X.  Y)  in  a  neighborhood  of  x  =  0,y  =  0,a(0,0)  =  0 
0(0,  0)  =  0.  The  functions  a,  0  are  C* -functions  of  X,  Y.  By  a  translation  of  variables,  we 
can  therefore  assume  that  this  equilibrium  point  is  zero.  Thus,  we  will  assume  without  loss 
of  generality  that 

(5.3)  X(0. 0)  =  0,  no,  0)  =  0. 


Our  objective  is  to  study  how  the  stability  properties  of  the  unperturbed  equation 
(5.1)  are  determined  by  the  nonlinearities  X0,  Y0.  Also,  we  determine  the  number  and 
stability  properties  of  the  periodic  orbits  that  can  bifurcate  from  zero  as  the  perturbation 
terms  X,  Y  are  varied  in  a  neighborhood  of  zero.  The  existence  and  the  number  of  periodic 
solutions  follow  from  the  bifurcation  function  obtained  by  an  application  of  the  method  of 
Liapunov-Schmidt.  We  prove  that  this  bifurcation  function  also  carries  the  information  on 
the  stability  of  periodic  orbits.  The  proof  of  this  latter  fact  uses  the  center  manifold 
theorem. 
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Li  mma  5.1.  There  are  a  neighborhood  V  of  (x,  y)  =  (0.  0)  and  a  neighborhood 
W  of  ( X ,  Y  )  =  (0,  0)  such  that,  for  ( X ,  Y )  £  W  and  any  periodic  orbit  7  =  y(  X,  Y  )  of 
(5.2)  in  V,  the  stable  manifold  has  dimension  >  n  +  1 .  Furthermore,  there  is  an  (n  +  !  )- 
dimensional  manifold  S{ 7,  X,  Y)  ( either  a  generalized  Mobius  band  or  the  cross-product 
of  a  circle  and  an  n-dimensional  ball)  containing  7  such  that  7  is  exponentially  asymptot¬ 
ically  orbitally  stable  with  asymptotic  phase  relative  to  initial  values  in  S( 7.  X,  Y ).  Fur¬ 
thermore,  S( 7.  X,  Y)  is  Ck  in  X,  Y  and  is  diffeomorphic  to  the  local  stable  manifold 
for  any  periodic  orbit  of  the  linear  equation  x  =  Bx,  y  =  Ay. 

The  proof  of  this  result  may  be  supplied  by  using  the  ideas  for  the  construction  of 
stable  manifolds  near  periodic  orbits  in  Coddington  and  Levinson  [1),  Hale  [1]  (see.  also. 
Fenichel  [  1 1 ). 

The  manifold  S(7,  X,  Y)  in  Lemma  5.1  has  codimension  one.  Thus,  the  complete  be¬ 
havior  of  the  solutions  near  a  periodic  orbit  is  determined  by  what  happens  in  one  other 
direction.  An  analysis  of  the  Poincare  map  on  a  center  manifold  will  take  care  of  this  direc¬ 
tion. 

There  is  a  center  manifold  for  equaiton  (5.2)  which  is  Ck  in  ( x ,  X,  Y)  in  a  sufficiently 
small  neighborhood  V  of  ( x ,  y)  =  (0,  0)  and  a  neighborhood  W  of  (A,  Y)  ~  (0,  0).  If  a 
center  manifold  is  given  by  S  =  {(x,  y):  y  =  <Jt(x,  X.  T)}  then  the  flow  on  the  center  mani¬ 
fold  S  is  given  by 


(5.4) 


x  =  Bx  +  X{x). 

X{x)  =  XQ(x.  Hx.  x,  n)  +  X(x,  *(*.  X,  K)). 


If  7  is  a  periodic  orbit  of  (5.2),  then  7  belongs  to  S  and  conversely.  Titus,  the  periodic- 
orbits  of  (5.2)  can  be  determined  by  discussing  the  periodic  orbits  of  (5.4).  Also,  from 
Lemma  5.1,  the  stability  properties  of  a  periodic  orbit  of  (5.2)  are  completely  determined 
by  the  stability  properties  of  the  corresponding  periodic  orbit  of  (5.4). 

By  the  manner  in  which  a  center  manifold  is  constructed,  one  can  also  obtain  a  priori 
bounds  on  \j/.  In  fact,  for  any  e  >  0,  5  >  0,  there  is  a  constant  K  =  A'(e,  6)  such  that,  if 
y  =  ip(x,  X,  Y )  is  a  center  manifold  of  (5.2)  for  lx|  <  e.  (A-.  Y)  in  a  5-neighborhood  of 
zero. then 


(5.5)  h/t(x,  X.  Y)l  <  Xlxl. 

Also,  K(t,  5)  — ♦  0  as  (e,  6)  — *•  (0,  0).  Thus,  in  a  neighborhood  of  (x,  y)  =  0,  (A',  Y)  =  0, 
every  periodic  orbit  of  (5.2)  must  have  the  y  coordinate  satisfying  (5.5).  For  x  =  (x,.  x2), 
this  justifies  the  transformation  of  variables 


(5.6)  x,  =  u  cos  0,  x2  -  ~u  sin  9,  y  =  uv 

in  (5.2).  The  new  equations  for  (6.  u,  v)  are 
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0=1-  (A,  sin  cos  d)lu  =  1  +  0(0,  u,  v.  A,  K), 

(5.7)  u  =  3*,  cos  0  -  A2 

6  =  Av  +  Y/u  -  (A,  cos  0  -  A2  sin  0)v/u 

where  A'  =  (A*,,  Aj),  X  =  X0  +  X,  Y  =  Y0  +  Y  and  all  functions  are  evaluated  at  (u  cos  0 . 
-u  sin  0,  mu). 

Since  the  function  0  satisfies  0(0, 0,  u,  0,  0)  =  0,  it  follows  that  0  >  'A  for  u  in  a 
fixed  compact  set  and  (u,  X,  T)  in  a  sufficiently  small  neighborhood  of  (0,  0,  0).  Thus,  we 
may  replace  t  by  0  to  obtain 

(5.8)  du/dd  =  /(0,  u,  u,  A,  K),  <ib/</0  =  /(u  +  g(0.  u,  u,  A,  T). 

The  functions  /,  g  are  2rr-periodic  in  0  and 

(5.9)  /(0,  0.  u.  A,  r)  =  0,  3/(0,  0,  u,  0,  0)/du  =  0,  g(6.  0,  u,  0,  0)  =  0. 

Any  2ir-periodic  solution  of  equation  (5.8)  corresponds  to  a  periodic  orbit  of  (5.2) 
through  the  transformation  (5.6)  and  conversely. 

For  equation  (5  .8),  there  is  the  standard  procedure  of  alternative  problems  or  Liapunov 
Schmidt  for  determining  the  2rr-periodic  solutions  for  ( u ,  A,  Y)  in  a  neighborhood  of  zero 
and  o  in  a  compact  set  (see  Cesari  (1|  for  a  general  discussion  of  the  alternative  method  as 
well  as  references).  To  describe  the  method,  let  Pln  =  {w:  R  — ♦  R"  +  1 ,  27r-periodic,  con¬ 
tinuous)  with  the  supremum  topology  and,  for  e  >  0,  let  W(e)  be  the  e-neighborhood  of 
(A.  Y)  =  0.  For  a  €  R  fixed,  consider  the  equation,  for  (u.  u)  6  P2n, 

i/o’ u{s)ds  =  * 

(5  10>  —  =  ne,  u,  v,  X.Y)-^  u(s),  u(s),  A.  Y)ds. 

dfQ  =  Bv  +  g(0,  m.  u.  A.  Y) 

An  application  of  the  Implicit  Function  Theorem  shows  there  are  an  e  >  0  and  a 
unique  C*-function  ( u(-,a ,  A,  Y ),  u(-,  a.  A,  T))  G  P2n  for  lo|  <  e,  (A,  Y)  G  K'fe),  vanishing 
for  (a.  A,  T)  =  (0.  0,  0)  and  satisfying  (5.10).  If  we  define 

(5.11)  G(a,  A,  Y)  =  ^~  (’"As,  u(s,  a,  X.  Y),  u(s,  o.  A,  Y),  A,  Y)ds 

then  the  above  2ff-periodic  function  will  satisfy  (5.8)  if  and  only  if 

(5.12)  G(a,  A.  Y)  =  0. 

It  is  also  easy  to  show  that  every  2n-periodic  solution  of  (5.8)  for  (n,  A,  )’)  in  a  neighbor¬ 
hood  of  zero  and  u  in  a  compact  set  can  be  obtained  through  this  process. 
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The  function  G(a,  X,  Y)  is  called  the  bifurcation  function  and  equation  (5.12)  is  called 
the  bifurcation  equation. 

We  can  now  prove  the  following  result.  The  conclusions  on  stability  are  due  to  deOli- 
veira  and  Hale  ( I  j . 

Theorem  5.2.  There  are  a  neighborhood  U  of  zero  in  P2  „ ,  an  e  >  0  and  Ck -functions 
u:  R  x  V(e)  x  H'(e)  — ►  R, 
v:  R  x  F(e)  x  W(e)  — *  Rn, 

G  V{e)  x  W(e)  —  R, 

Re)  =  (~e,  e),  W(f)  the  e-neighborhood  of  (X.  Y)  -  (0.  0).  (u(\  a,  X.  )'),  o(  ,  a.  X,  ))) 

£  P2,  such  that  equation  (5.8)  has  a  2it-periodic  solution  (u{0).  v(0))  in  U  if  and  only  if 
(u(8),  v(0))  =  (u(6,  a,  X,  Y).  u(0.  a,  X,  T))  and  a  satisfies  (5.12).  Furthermore,  the  stability 
properties  of  a  2  tt-periodic  solution  u(8.  aQ.  X,  >’),  i40,  a0.  X.  Y)  coincide  with  the  stability 
properties  of  a0  as  a  solution  of  the  scalar  equation 

(5.13)  a  =  G(a.  X,  Y). 

In  particular.  u(8.  a.  X,  Y),  v(8.  a,  X.  K)  is  stable  (asymptotically  stable)  ( unstable )  if  and 
only  if  a0  is  stable  (asymptotically  stable )  (unstable). 

PiftooE.  We  only  outline  the  proof  of  stability.  The  details  can  be  found  in  deOliveira 
and  Hale  { 1 }  or  Hale  [2] .  The  theorem  is  first  proved  for  the  variable  t>  absent  ;  that  is,  a 
scalar  equation 

(5.14)  du/dd  =  f(0.  u,  X.  Y). 

Let  71(0.  a.  X.  T),  G(a,  X.  T)  be  the  27r-periodic  function  and  bifurcation  function  construc¬ 
ted  for  this  equation  by  the  alternative  method.  Then  dt7(0,  a.  X,  Y)/da  =  1  at  a  =  0. 

(3f,  Y)  =  (0,  0)  and  one  can  make  the  transformation  of  variables  u  i~«  b.  «  =  u(8,  b.  X.  T) 
to  obtain 

Since  3 2736  >  0  if  b,  X  are  small,  one  obtains  the  result  for  u  absent.  When  u  is  present, 
there  is  a  center  manifold  with  the  flow  on  the  center  manifold  given  by  an  equation  of  the 
form  (5.14).  One  now  has  two  bifurcation  functions  G(a.  X,  T)  and  G(a,  X,  V')  correspond¬ 
ing  respectively  to  (5.14)  and  (5.8).  These  functions  must  have  the  same  /.eros  and  one  can 
use  an  argument  similar  to  the  one  in  the  proof  of  Theorem  4.2. 

Theorem  5.2  has  obvious  interpretations  about  the  existence  and  stability  of  periodic 
orbits  of  (5.2).  Using  the  same  type  of  argument  that  was  used  in  the  proof  of  Theorem  4.2, 
one  can  also  show  that  the  w-limit  set  of  an  orbit  in  a  neighborhood  of  zero  consists  of 
either  zero  or  a  single  periodic  orbit.  These  results  are  summarized  in 
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Theorkm  5.3.  Let  u(8,  a,  X,  Y),  iAfi,  a,  X,  Y)  be  the  ln-periodic  functions  satis¬ 
fying  (5.10)  and  let  G(a,  X,  Y)  be  defined  in  (5.11).  Then  there  are  a  neighborhood  V  of 
(x,  y)  =  (0,  0)  and  a  neighborhood  W  of  (X,  X)  =  (0,  0)  such  that  equation  (5.2)  for  (X,  X) 
in  W  has  a  periodic  solution  (jr(r),  y(r))  in  V  if  and  only  if 

jc(r)  =  (u(0(f),  a0,  X,  X)  cos  0(f), -u(0(t),  a0.  X,  X)sin  0(f)), 

At)  =  A0(t),  a0,  X,  X)u(0(f),  aQ,  X,  X) 

where  G(a0.  X,  X)  =  0,  0(f)  =  I  +  0(0,  a0,  X,  X),  u(0,  a0,  X,  X))  and  0(0),  aQ  are  uniquely 
determined  by  x(0),  _v(0).  The  stability  properties  of  the  periodic  orbit  coincide  with  the 
stability  properties  of  a0  as  a  solution  of  (5.13).  Finally,  the  co-limit  set  of  an  orbit  of 
(5.2)  in  V  is  a  single  periodic  orbit. 

One  can  show  (see,  for  example,  Chafee  [1))  that 

(5.15)  G(a,X,Y)=-G(-a.X.Y). 

Thus,  if  a  is  a  solution  of  G(a,  X,  X)  =  0,  then  so  is  -  a.  One  shows  that  these  two  solu¬ 
tions  will  correspond  to  the  same  periodic  orbit  of  (5.2).  Thus,  we  only  need  to  be  con¬ 
cerned  with  positive  roots  of  G(a,  X,  X)  =  0. 

Using  this  fact,  Theorem  5.2  and  the  ideas  of  the  proof  of  Corollary  4.4,  one  obtains 
the  following  interesting  consequence  for  the  stability  of  the  zero  solution  of  the  unper¬ 
turbed  equation  (5.1). 

Corollary  5.4.  The  zero  solution  of  (5.1 )  is  asymptotically  stable  ( unstable )  if  and 
only  if  there  is  an  e  >  0  such  that  aG(a,  0,  0)  <  0  (>  0)  for  0  <  |a|  <  e.  There  is  an  e  >  0 
such  that  G(a,  0,  0)  =  0  for  |a|  <  e  if  and  only  if  the  zero  solution  is  stable  and  there  is  a 
first  integral  in  a  neighborhood  of  zero. 

This  corollary  was  first  proved  by  Liapunov  [1}  (see  also  Bibikov  [1])  for  the  case 
when  X0,  X0  are  analytic.  In  particular.  Corollary  5.4  implies  that,  if 

(5.16)  G(a.  0.  0)  =  P0a2q  *  '  +  AM2"  +  ' )  as  |al  —  0,  J30  *  0. 

then  the  zero  solution  of  (5.1 )  is  asymptotically  stable  if  0O  <  0  and  unstable  if  /30  >  0.  By 
using  the  theory  of  normal  forms  for  equation  (5.1)  (see,  for  example,  Bibikov  [1]  or 
Takens  [ I J )  one  can  relate  the  sign  of  0O  to  the  stability  of  equation  (5.1)  under  certain 
types  of  high  order  perturbation  (for  details,  see  Negrini  and  Salvadori  [1  j ,  Bcrnfeld  and 
Salvadori  (l|). 

For  k  >  2q  +  I ,  we  now  discuss  equation  (5.2)  with  (7(o,  0.  0)  satisfying  (5.16).  There 
are  functions  a;:  W  — ►  R,  or^( 0,  0)  =  0,  0  </'  <  q  -  1 ,  0^(0,  0)  =  p0  such  that 

G(a,  X,  X)=  J  afX,  Y)a2**  1  +  o(|<z|2«  +  ’) 

/= o 


(5.17) 
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as  \a\  — *  0.  Let 

T°  =  {(X  Y).  a0(X,  Y)  =  0],  r4  =  {(X.  Y):  a0(AT,  Y)  >  0), 

(5.18) 

r  =  {(AT.  Y):  a0(X.  K)  <  0} . 

All  of  the  coefficients  a^X,  T)  can  be  computed  in  terms  of  the  derivatives  of  X,  Y  of  or¬ 
der  <  2/  +  l.  For  the  matrix  of  the  linear  approximation  of  (5.2)  at  zero,  the  function 
aQ(X,  Y)  represents  the  real  parts  of  the  eigenvalues  which  are  purely  imaginary  for  (  V.  Y  ) 

=  (0,  0).  Thus,  in  order  to  have  a  periodic  orbit  bifurcate  from  zero  by  a  variation  in 
(X.  Y)  near  some  point  (Af,  T),  the  point  (Af,  Y)  must  have  a0(Af,  T)  =  0,  that  is  (A",  Y) 
belongs  to  F°. 

For  some  special  types  of  one  parameter  families  of  vector  fields,  the  bifurcation  at 
(Af0,  F0)  is  very  simple  as  shown  in  the  following  result,  essentially  due  to  Negrini  and  Sal¬ 
vador  [  1  ] . 

Theorem  5.5.  Suppose  G(a,  0,  0)  satisfies  (5.16),  (Af(p).  )'(/ j )).  p  £  R,  is  a  C 1  -family 

clef 

of  vector  fields,  vanishing  at  p  =  0,  and  suppose  the  functions  aAXlp).  Y(p))  =  a .(  p )  arc 
defined  in  (5.17).  If  a^O)  0,  then  there  are  a  neighborhood  V  of  lx,  y)  =  (0.  0)  and  a 
p0  >  0  such  that,  if  0  <  |p|  <  p0,  then  equation  (4.2)  has  a  periodic  orbit  in  V  if  and  only 
if  Og(0)/30  <  0.  When  this  condition  is  satisfied  the  orbit  is  unique.  It  is  asymptotically 
stable  if  j}0  <0  and  unstable  if  1 30  >  0. 

Proof.  The  bifurcation  function  G{a,  X{p),  T(p))  =  G(a,  p)  satisfies 

G{a,  p)  =  £  a ,(p)a2/+  1  +  <d> \7q  +  1 )  as  |a|  —  0 
/— o 

where  a;(0)  =  0,  0  </  <  q,  0^(0)  =  (?0  ^  0.  It  is  not  difficult  to  show  that  there  is  a  neigh¬ 
borhood  U  of  (a,  p)  =  (0,  0)  such  that  any  solution  of  G(a,  p)  =  0  in  U  has  the  a  priori 
bound  la|  <  Ac|jat| 1  f .  If  we  let  a  =  |p| l,2qb,  then 

±-(C(\p\'/2«b,  p)  =  F{b,  p)  =f  £  y.(p)b2'  +  1  +  o(lpl) 

ImI  /to  ' 

as  Ipl  —*  0,  where  -y0(0)  =  a q(0),  y^O)  =  0,  0  <  /  <  q,  yql0)  -  pn  =*  0.  The  Implicit 
Function  Theorem  implies  the  result  stated  in  the  theorem  on  existence.  The  stability  fol¬ 
lows  from  Theorem  5.2. 

Theorem  5.5  says  that  nothing  very  complicated  can  occur  from  variations  in  a  one 
parameter  family  of  vector  fields  regardless  of  the  degeneracy  in  G(a.  0.  0)  provided  that  the 
eigenvalues  cross  the  imaginary  axis  with  a  definite  speed;  that  is,  a0(p)  =  Qq(0)p  +  n(|p|) 
as  p  — 1 ►  0,  <*0(0)  *=  0.  It  is  only  possible  to  obtain  one  periodic  orbit  and  the  amplitude  of 
the  orbit  (front  the  proof)  is  of  order  |pt| 1  ^ as  Ipl  — *  0. 

Theorem  5.5  does  not  imply  that  there  is  at  most  tine  periodic  orbit  of  (5.2)  for  every 
(X,  Y)  in  a  neighborhood  of  (0,  0).  lit  fact,  if  q  >  1 ,  we  will  see  later  that  there  will  always 
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be  some  (AT,  T)  in  any  neighborhood  of  zero  such  that  there  are  q  periodic  orbits.  Assum¬ 
ing  one  parameter  families  of  perturbations  in  (S.2),  Flockerzi  |1]  has  used  Newton’s  poly¬ 
gon  to  determine  the  number  of  solutions  of  G(a,  p)  =  0  and  their  stability  when  a0(p)  = 
OqU*  +  o(lul*)  as  lul  — ►  0.  We  obtain  results  for  arbitrary  perturbations  (AT,  Y)  which  are 
more  in  the  spirit  of  Chafee  (1] . 

The  first  result  for  q  =  1  is  referred  to  as  the  generic  Hopf  Bifurcation  Theorem. 

Theorem  5.6.  Suppose 

G{a,  0,  0)  =  0oa3  +  o(|a|3)  as  |a|  — ►  0,  fS0  *  0, 

and  a0(AT,  Y)  is  defined  in  (5.17).  Then  there  are  neighborhoods  W  of  (AT,  Y)  =  (0,  0),  V 
of  (x.  A1)  =  (0,  0)  such  that  for  (AT,  Y)  in  W,  equation  (5.2)  has  a  periodic  orbit  in  V  if  and 
only  if  a0(AT,  T)( 30  <  0.  When  this  condition  is  satisfied,  the  orbit  is  unique  and  is  asymp¬ 
totically  stable  ( unstable )  if  and  only  if  0O  <  0  (>  0). 

Proof.  In  a  neighborhood  of  zero,  G(x,  AT,  K)  =  0  either  has  one  positive  root  or 
no  positive  root.  The  condition  for  the  existence  of  one  is  a0(AT,  y)0o  <  0.  The  stability 
of  the  orbit  follows  from  Theorem  5.2. 

With  the  notation  as  in  (5.18),  if  (AT.  Y)  £  r  +  (AT,  Y),  then  a  periodic  orbit  can  exist 
for  (AT,  Y)  in  r+(AT,  Y)  if  and  only  if  0O  <  0  and  then  it  is  asymptotically  stable  (the  super¬ 
critical  case)  with  the  origin  being  unstable.  It  can  exist  for  (AT,  K)  £  T"  if  and  only  if 
(70  >  0  and  then  it  is  unstable  (the  subcritical  case)  with  the  origin  being  stable. 

Suppose  there  is  a  one  parameter  family  of  perturbations  Af(X),  X(X),  X  £  R,  which 
vanish  for  X  =  0  and  this  curve  of  perturbations  has  (AT(X),  T(X))  in  f  for  -e  <  X  <  0  and 
in  T+  for  0  <  X  <  e.  Then  the  origin  is  asymptotically  stable  for  X  <  0  and  unstable  for 
X  >  0.  There  are  an  unstable  periodic  orbit  and  subcritical  bifurcation  at  X  =  0  if  0O  >  0 
and  an  asymptotically  stable  periodic  orbit  and  supercritical  bifurcation  at  X  =  0  if  0O  <  0. 
Notice  that  it  is  not  required  that  the  curve  cross  r°  transversally;  that  is,  the  eigenvalues  of 
the  linear  part  of  (5.2)  cross  the  imaginary  axis  with  a  definite  speed.  Contrast  this  with 
Theorem  5.5. 

If  the  first,  second  and  third  derivatives  of  G(a,  0,  0)  vanish  at  a  =  0,  Theorem  5.6 
cannot  be  used.  The  discussion  of  the  case  of  higher  order  degeneracy  of  G(a,  0.  0)  is  much 
more  complicated  and  one  is  forced  to  use  the  Malgrange  Preparation  Theorem  on  G(a,  AT,  >') 
for  the  general  case  as  we  did  in  §4.  There  are  some  more  specific  results  presented  below 
that  can  be  obtained  without  this  theorem;  namely  the  case  q  —  2  in  (5.16). 

Suppose  q  =  2  in  (S.16)  and  define  T°,  r+,  T~  by  (5.18).  if  o0(AT,  Y)0O  <  0.  then 
a0(Af,  y)a3(Af,  T)  <  0  for  (AT,  Y)  in  a  neighborhood  of  zero  and  G(a,  X,  Y)  has  exactly 
one  positive  zero.  Thus,  equation  (5.2)  has  exactly  one  periodic  orbit.  If  a0(AT,  >’)0O  >  0. 
there  may  be  more  than  one  positive  zero  of  G(a,  AT,  Y).  To  determine  conditions  for  when 
there  are  more  than  one,  define  P(r,  AT,  Y),  r>  0,  by  the  relation 


(5.19) 


G(a,  AT,  Y)  =  af\a\  AT,  Y). 
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Extend  P(r,  X,  Y)  as  an  even  function  for  r  £  R.  Then  the  positive  zeros  of  G  are  deter¬ 
mined  from  the  positive  zeros  of  f\r,  X,  Y).  Since  p0  #  0,  the  Implicit  Function  Theorem 
implies  there  are  neighborhoods  V  of  r  =  0,  W  of  ( X ,  K)  =  0  and  a  unique  function  6: 

W  — ►  R+  such  that  5(0, 0)  =  0  and  dP(&(X,  Y),  X,  Y)(dr  =  0.  Let 

(5  20)  y(X,  Y)  =  />(5(X,  Y),  X,  Y). 

If  y(X,  Y)P0  <  0,  there  are  two  positive  simple  zeros  of  P(r,  X,  T)  and  thus  two  periodic 
orbits  of  equation  (5.2).  The  orbit  corresponding  to  the  smaller  zero  is  unstable  (asymptot¬ 
ically  stable)  if  p0  <  0  (>  0)  and  the  other  is  asymptotically  stable  (unstable).  If  y(X,  Y)P0 
=  0,  S(X,  Y)  >  0,  then  there  is  a  unique  periodic  orbit  which  on  a  center  manifold  is  stable 
on  one  side  and  unstable  on  the  other.  If  y(X,  Y)P0  >  0,  there  is  no  positive  zero  of 
P(r,  X,  Y)  and,  thus,  no  periodic  orbit  of  (5.2).  These  results  are  summarized  in 

Theorem  5.7.  Suppose 

G(a.  0,  0)  =  (30i7s  +  <?(|d|5)  as  la]  — *  0,  0O  s*  0, 

a0(Af,  K)  is  defined  in  (5.17)  and  y(X,  Y)  is  defined  by  (5.20)  when  a0(X,  Y)P()  >  0.  Then 
there  are  neighborhoods  W  of  (X,  Y)  =  (0,  0),  V  of  (jc,  v)  =  (0.0)  such  that,  for  ( X ,  K)  £  IV. 
equation  (5  .2)  satisfies  the  following  properties  in  V: 

(i)  a0(Af,  Y)P0  <  0  implies  a  unique  periodic  orbit  which  is  asymptotically  stable 
(unstable)  if  P0  <  0  (>  0). 

(ii)  a0(Ar,  Y)p0  >  0,  y(X,  K)/30  <  0  implies  there  are  two  periodic  orbits,  the  one 
corresponding  to  the  smaller  value  of  a  is  unstable  (asymptotically  stable)  with  the  other 
one  being  asymptotically  stable  (unstable)  if  p0  <  0  (>  0). 

(iii)  aQ(X,  y)/30  <  0,  y(X,  T)  =  0,  h(X,  T)  >  0  implies  there  is  a  unique  periodic 
orbit  which  on  a  center  manifold  is  stable  on  one  side  and  unstable  on  the  other. 

(iv)  a0(X,  Y)P0  >  0,  y(X,  y)0o  >  0  implies  no  periodic  orbit. 

If  we  let  2  =  {(AT,  Y):  6(Ar,  T)  =  0},  2  has  codimension  1  and  is  tangent  to  T0  at 
(0,  0).  Define  P°  as  in  (5.18).  Then  Theorem  5.7  states  that  the  bifurcation  set  for  equa¬ 
tion  (5.2)  for  q  =  2  in  (5.16)  is  given  by  l’°  U  2.  Tire  neighborhood  W  is  divided  into  the 
components  shown  in  Figure  4  with  the  number  of  periodic  orbits  indicated.  The  set  r+  is 
drawn  to  the  right  of  T0  and  the  intersection  of  F0.  2  contains  the  point  (X,  T)  =  (0,  0). 


Figure  4a.  P0  <0.  Figure  4b.  p0  >  0. 
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From  these  pictures,  one  can  see  clearly  the  meaning  of  Theorem  5.7  for  q  -  2.  Any 
curve  of  perturbations  (A'lp),  T(p))  crossing  P°  transversally,  that  is,  P'(0)  ^  0,  can  never 
enter  the  region  containing  two  orbits.  If  the  curve  is  not  transversal  to  P°.  it  is  possible 
for  it  to  enter  this  region.  As  remarked  earlier,  this  type  of  problem  with  one  parameter 
problems  was  discussed  by  Flocker/.i  [1 1 . 

There  remains  the  problem  of  computing  the  number  p0  in  (5.1b)  and  the  functions 
otj(X,  Y)  in  (5.17).  One  method  is  to  approximate  the  Taylor  series  in  a  for  the  functions 
u(a,  X,  Y),  u(a,  X,  >'),  G[a,  X,  Y)  directly  from  the  defining  equations  (5.10),  (5.1 1 ).  An 
alternative  procedure  is  to  transform  the  vector  field  (5.2)  by  a  change  of  coordinates  to  a 
convenient  normal  form.  The  method  of  Lie  transforms  developed  originally  for  applica¬ 
tions  to  celestial  mechanics  is  very  convenient  for  obtaining  this  normal  form.  The  best  tel- 
erence  is  Henrard  [1|.  This  method  also  is  presented  in  the  forthcoming  book  of  Chow  and 
Hale  [1 1 .  It  leads  naturally  to  implementation  on  the  computer  (see  Meyer  and  Schmidt  1 1  ]  ) 
and  could  lead  to  considerable  improvement  over  the  present  methods  used  for  the  computa¬ 
tion  of  periodic  orbits. 

Using  the  same  ideas,  the  appropriate  extensions  of  the  results  to  the  case  where  A 
has  eigenvalues  with  nonzero  real  parts  are  only  a  technicality. 

In  the  applications,  it  often  happens  that  the  perturbed  vector  field  depends  in  a  spe¬ 
cial  way  on  a  finite  number  of  parameters  varying  in  a  neighborhood  of  /ero  The  bifurca¬ 
tion  diagrams  in  the  parameter  space  are  then  obtained  by  studying  the  way  in  which  the 
family  of  vector  fields  crosses  the  general  bifurcation  curves  given  above.  For  example,  if 
there  is  a  single  parameter  p  and  the  lowest  order  terms  in  the  bitincation  function  are 
pa3  +  y\  with  p  +  0,  y  0,  then  this  is  the  case  of  generic  Hopf  bifurcation,  and  the  de¬ 
pendence  on  X  through  the  linear  term  7X  implies  there  is  only  one  crossing  of  the  bifurca¬ 
tion  curve  as  we  noted  before.  When  there  are  several  parameters  in  the  generic  Hopf  bifur¬ 
cation  problem,  then  the  family  of  vector  fields  can  cross  seveial  times  and  in  a  vety  com¬ 
plicated  manner.  The  case  where  the  degeneracy  in  a  is  ol  order  5  (</  -  2  in  (5.1b))  is  even 
more  difficult.  Also,  it  is  often  the  case  that  one  parameter,  say  p.  plays  a  distinguished 
role  (sometimes  called  a  modal  parameter)  and  one  is  interested  in  studying  the  bifurcations 
when  this  parameter  is  fixed  and  the  other  parameters  arc  varied.  This  question  has  been 
discussed  in  detail  for  third  order  and  fifth  order  degeneracy  by  (lolubitsky  and  Langford  1 1 1 . 

i 

1 
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6.  First  order  bifurcation  in  the  plane 

Suppose  ft  C  R2  is  open  and  let  X2  =  X^(^)  be  the  set  of  C'r(ft,  R2)  vector  fields. 
r  >  1,  which  are  transverse  to  the  boundary  of  ft.  The  set  of  structurally  stable  vector 
fields  in  Xj  is  characterized  in  Theorem  3.2  and  1.2  is  open  and  dense  in  X 2.  Any  /£  Xr, 
which  is  not  structurally  stable  is  a  bifurcation  point.  Our  objective  in  this  section  is  to  dis¬ 
cuss  some  of  the  most  elementary  bifurcation  points. 

Following  Andronov  and  Leontovich  [1]  (see  also  Andronov  et  al.  [1])  except  for 
terminology,  we  say  /£  X2  is  a  bifurcation  point  of  degree  0  if  it  is  structurally  stable.  It 
is  a  bifurcation  point  of  degree  1  if  it  is  not  of  degree  zero  and  every  £  in  a  neighborhood 
of  /  is  either  of  degree  zero  or  equivalent  to  /.  It  is  a  bifurcation  point  of  degree  2  if  it  is 
not  of  degree  0  or  1  and  every  £  in  a  neighborhood  of /is  either  of  degree  0  or  1  or  equiva¬ 
lent  to  f.  Similarly,  one  defines  a  bifurcation  point  of  degree  k. 

The  bifurcation  points  of  degree  one  are  the  simplest  types  that  can  occur.  They  cor¬ 
respond  to  those  vector  fields  which  satisfy  all  of  the  conditions  for  structural  stability  ex¬ 
cept  one  and  the  condition  that  is  violated  must  be  done  so  in  the  simplest  way.  To  be 
more  specific,  suppose  an  equilibrium  point  x0  =  0  of  /is  not  hyperbolic.  This  can  occur 
only  when  A  =  3/(0)/3x  has  an  eigenvalue  on  the  imaginary  axis.  Now  consider  all  £  near 
/which  have  an  equilibrium  point  near  zero  which  is  not  hyperbolic.  In  order  for  a!)  such 
£  to  be  equivalent  to  f  it  is  necessary  that  A  has  either  a  simple  zero  or  a  pair  of  simple 
purely  imaginary  eigenvalues.  If  A  has  a  simple  zero  and  /  has  degree  1  in  X2.  r  >  2.  then 
the  bifurcation  function  G(a,  /)  in  §4  for  equilibrium  points  near  x0  =  0  must  have  the 
form  G(a,  f)  =  0oa2  +  o(|a|)  as  |a|  — ►  0,  /30  =£  0,  and  Theorem  4.8  is  applicable.  There 
are  a  neighborhood  W  of  /  and  a  submanifold  T  of  codimension  one  such  that  (V\F  = 

(/,  U  U2  where  g  €  (7,  implies  no  equilibrium  point  near  x0,  £  6  U2  implies  a  saddle  and 
hyperbolic  node  near  x0  and  g  6  r  implies  a  saddle-node  near  x0,  which  we  call  an  elemen¬ 
tary  saddle-node. 

If  A  has  a  pair  of  purely  imaginary  roots  and  /  has  degree  1  in  X\,  r  >  3,  then  the 
bifurcation  function  G(a.  f)  in  §5  for  the  existence  of  periodic  orbits  near  x0  must  have 
the  form 

G(a,  f)  =  0oa3  +  oQal3)  as  |a|  —  0,  0O  *  0. 

This  corresponds  to  the  generic  Hopf  bifurcation  (Theorem  5.6).  There  are  a  neighborhood 
W  of  /  and  a  submanifold  T  of  codimension  one  such  that  K,\T  =  Uf  U  U2  and  ££(/, 


HtECKMWC  PACK  BLANK -NOT 


JACK  K.  HALE 


32 

implies  no  periodic  orbit  near  x0  and  g  €  U2  implies  a  unique  periodic  orbit  near  x0.  We 
refer  to  the  equilibrium  point  of  g  near  x0  as  an  elementary  focus  if  g  €  T. 

The  condition  f  r  f  to  be  a  bifurcation  point  also  can  occur  when  a  periodic  orbit  7 
becomes  not  hyperbolic.  Let  7  =  {0(r,  p).  (£R)  where  0(r,  p),  0(0,  p)  =  p,  is  a  periodic 
solution  of  x  =  f(x).  Let  L0  be  a  transversal  to  7  at  p  and  it:  L0  — *•  L0  be  the  Poincare 
map.  Then  ir(p)  =  p  and  not  hyperbolic  implies  rr'(p)  =  1 .  If  /  has  degree  one  in  X2. 
r  >  2,  then  it  must  necessarily  be  true  that  it"(p)  ^  0.  This  means  that  rr(p  +  u)  -  (p  +  u) 
in  a  neighborhood  of  u  =  0  behaves  like  a  quadratic  function  in  u.  Then  there  will  be  a 
neighborhood  W  of  f  and  a  submanifold  T  of  codimension  one  such  that  W\r  =  U{  U  U2 
and  g  €  Ut  implies  no  periodic  orbits  near  7  (no  fixed  points  of  it  near  p)  and  g  &  U2  im¬ 
plies  two  hyperbolic  periodic  orbits  near  7  (two  fixed  points  of  it  near  p),  one  unstable  and 
the  other  asymptotically  stable.  If  g  £  T,  then  there  is  a  unique  periodic  orbit  near  7,  stable 
on  one  side  and  unstable  on  the  other.  Having  g  £  T  is  equivalent  to  having  a  double  zero 
of  ir{p  +  u)  -  (p  +  u)  near  u  =  0.  This  result  on  periodic  orbits  can  also  be  obtained  from 
the  method  in  §4  by  introducing  a  local  coordinate  system  x  (0,  p),  x  =  0(0,  p)  +  pi>(0), 
where  u(0)  is  a  unit  vector  orthogonal  to  3 0(0,  p)/30.  If  t  is  replaced  by  0  in  the  new  equa¬ 
tions,  one  obtains  a  scalar  equation  for  p  as  a  periodic  function  of  0,  a  special  case  of  equa¬ 
tion  (5.8).  Theorem  5.2  is  then  applicable. 

If  the  equilibrium  points  and  periodic  orbits  are  hyperbolic  and  /  is  a  bifurcation  point, 
then  there  must  be  a  trajectory  connecting  saddle  points.  However,  if  /has  degree  one,  then 
one  can  show  that  there  must  be  an  orbit  whose  a  -  and  w-limit  sets  are  the  same  point-a 
homoclinic  orbit.  This  situation  is  more  complicated  to  understand  than  the  previous  ones 
because  it  is  a  global  problem  for  which  the  discussion  cannot  be  restricted  to  the  considera¬ 
tion  of  only  the  fixed  points  of  a  map. 

To  describe  the  behavior  near  a  homoclinic  orbit,  suppose  zero  is  a  saddle  point  of  /, 
/(0)  =  0,  3/(0)/3.x  has  eigenvalues  X  <  0  <  p.  If  Wf,  are  the  stable  and  unstable  mani¬ 
folds  of  0,  then  there  is  a  homoclinic  orbit  through  zero  if  and  only  if  (WJ.O  W“)\{0)  0. 

H  W“,  p0  =£  0,  then  the  orbit  O(p0)  through  p0  approaches  zero  exponentially 
as  t  — ►  ±°°.  If  7  =  O(p0)  U  {0} ,  then  the  invariant  set  7  can  have  either  of  the  configura¬ 
tions  shown  in  Figure  5  with  respect  to  the  position  of  the  stable  and  unstable  manifolds  of 
zero.  In  order  to  be  specific,  we  suppose  the  situation  in  Figure  5a  occurs.  The  other  case 
is  discussed  in  a  similar  way. 


Figure  5a.  Figure  5b. 
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Let  L0  be  a  transversal  to  y  at  p0.  Let  Lq  (Lq)  be  that  part  of  L0  interior  (exterior) 
to  y  in  Figure  6.  There  is  a  neighborhood  U  of  p0  such  that,  for  any  p£(/n/,J,  there 
is  a  r  =  r{p)  >  0  with  x(t(p),  p)  €  LjJ,  x (f,  p)  &  L£  ,0  <  t  <  t(p),  where  ^r(/.  p)  is  the 
solution  of  jr  =  fix)  through  p.  If  tr0(p)  =  x(r(p),  p),  then,  for  any  integer  k,  positive  or 
negative,  one  can  define  trj(p)  on  some  subset  of  depending  on  k.  The  map  n0  is  like 
a  Poincard  map  relative  to  the  interior  of  y.  It  is  not  defined  at  p0  and  r(  p)  — *■  «>  as  p  — *•  p0 . 


Figure  6. 


Definition  6.1 .  The  homocttnic  orbit  is  asymptotically  stable  ( unstable )  if  p) 

— ►  p0  as  k  — ►  <*>  (k  — ►  -<»).  It  is  exponentially  asymptotically  stable  ( exponentially  un¬ 
stable)  if  the  approach  to  p0  is  exponential. 

One  can  show  that  the  stability  (or  instability)  of  y  is  exponential  if  and  only  if 
o0  =  tr  3/(0)/ 3jc  =  X  +  p  <  0  (or  >  0)  (see  Andronov  et  al.  [  1  ] ,  Chow  and  Hale  (1 1  ). 

There  is  a  neighborhood  W  of  /  such  that  each  g  €  W  has  a  unique  zero  near  x  =  0. 

By  translation  of  variables,  we  can  assume  g(0)  =  0.  Let  W*.  W“  be  the  stable  and  unstable 
manifolds  of  zero.  Since  these  manifolds  are  smooth  in  g,  there  are  points  pu  g  of  first  con¬ 
tact  with  L0  (pJg  of  last  contact  with  L0)  with  respect  to  increasing  r  so  that  one  of  the 
pictures  in  Figure  7  prevails. 


Figure  7b. 
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If  case  (b)  prevails,  then  g  has  a  homoclinic  orbit  through  0.  If  T  =  {g  €  W:  pu  = 
ps  ,  then  one  can  show  that  T  has  codimension  one  for  a  sufficiently  small  neighborhood 
W  of  /  If  mr  =  Ut  U  U2.  suppose  g  €  Ul  coincides  with  Figure  7a  and  g  £  U2  coin¬ 
cides  with  Figure  7c.  If  g  £  £/, ,  then  there  is  a  neighborhood  V  of  7  containing  only  the 
critical  point  zero,  such  that  every  solution  with  initial  value  on  Ig,  the  segment  between 
pu  and  ps  ,  leaves  V  in  negative  time.  Thus,  the  behavior  of  solutions  near  7  is  deter¬ 
mined  by  what  happens  to  the  positive  orbit  0+(pu  )  of  pu  .  If  7  is  asymptotically  stable, 
then  0+(pu  )  remains  in  V  and  the  Poincare-Bendixson  Theorem  implies  there  is  a  periodic 
orbit  in  V.  If  it  is  further  assumed  that  y  is  exponentially  asymptotically  stable,  then  there 
is  a  unique  periodic  orbit  in  V  (see  Andronov  et  al.  [1]  or  Chow  and  Hale  [1]).  If  7  is  ex¬ 
ponentially  unstable,  then  0+(pu  g)  leaves  V  in  finite  time. 

If  g  £  U2,  that  is,  Figure  7c  prevails,  then  the  behavior  of  solutions  near  7  is  deter¬ 
mined  by  the  negative  orbit  0~(psg)  of  p  If  7  is  unstable,  then  there  is  a  periodic  orbit 
in  V  and  it  is  unique  if  7  is  exponentially  unstable.  If  7  is  exponentially  asymptotically 
stable,  then  0~(pu  )  leaves  V  in  finite  time. 

Finally,  one  can  show  that  /  having  hyperbolic  equilibrium  points  and  periodic  orbits 
is  a  bifurcation  point  of  degree  1  if  and  only  if  o0  =  tr  9/(0)/9x  0;  that  is,  the  stability 

or  instability  of  7  is  exponential  (see  Andronov  et  al.  [1]).  These  results  are  summarized  in 

Theorem  6.2.  A  vector  field  f  e  X.2.  r  >  3,  is  a  bifurcation  point  of  degree  1  if  and 
only  if  there  are  a  neighborhood  W  of  f  and  a  submanifold  T  of  codimension  one  in  W  such 
that  W\r  =  (/,  U  U2  where  each  g  €  U/  is  structurally  stable  but  g  f  h  if  g  €  Ut.  Ii  £  U2. 
Then  only  one  of  the  following  situations  prevails : 

(i)  g  €  T  has  an  elementary  saddle-node  at  x0,  there  are  no  equilibrium  points  of  g 
near  jc0  if  g  G  (/,  and  a  saddle  and  node  near  x0  if  g  £  U2 . 

(ii)  g  €  T  has  an  elementary  focus  at  x0  ,  there  is  no  periodic  orbit  of  g  near  xQ  if 
g  £  17,  and  a  periodic  orbit  near  x0  if  g  &  U2- the  generic  Hopf  bifurcation. 

(iii)  g  £  T  has  a  periodic  orbit  7  which  is  stable  from  one  side,  unstable  from  the 
other,  g  G  {/,  has  no  periodic  orbit  near  7  and  g  £  U2  has  two  hyperbolic  periodic  orbits 
near  7. 

(iv)  cr0  =  tr  df(0)/dx  0,  g  £  T  has  a  homoclinic  orbit  containing  a  saddle  point  x0, 
get/,  has  a  saddle  near  x0  and  no  periodic  orbit  near  7,  g  &  U2  has  a  saddle  point  and  a 
unique  hyperbolic  periodic  orbit  near  7  which  coalesce  as  g  — *■  F. 

(v)  There  is  a  connection  between  two  distinct  saddle  points. 

Each  of  the  cases  (i)-(v)  is  shown  in  Figure  8. 

There  is  another  interesting  phenomenon  that  can  occur  with  a  bifurcation  point  of 
degree  I .  In  case  (i),  the  unstable  manifold  of  the  saddle  could  be  connected  to  the  stable 
manifold  of  the  node  to  form  a  closed  curve  as  in  Figure  9a.  They  may  then  coalesce  and 
disappear  to  form  a  hyperbolic  periodic  orbit. 
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Using  Theorem  6.2,  Sotomayor  |1]  has  proved  the  following  interesting  result  Let 
H  =  {g:  [0,  1  ]  ►  X  j,  g  a  C' -function,  r  >  5}  with  the  topology  of  uniform  convergence 

of  the  function  and  its  derivatives.  Let  Ht  C  H  be  the  subset  of  functions  which  have  the 

property  that,  for  any  g  e  //, ,  there  is  a  finite  set  of  points  p, . pm  in  [0,  1 1  such  that 

g(p)  is  structurally  stable  for  p  ¥=  p^,  ginf)  is  a  bifurcation  point  of  degree  1  in  X  j  and  the 
curve  defined  by  g  is  transversal  to  the  set  1^  of  Theorem  6.2  corresponding  to  the  bifurca¬ 
tion  point  g{Uj).  Sotomayor  [  1  ]  has  proved  that  //,  is  residual  in  H,  that  is,  one  can  as¬ 
sume  that  the  bifurcation  points  on  a  given  curve  of  vector  fields  genetically  will  be  of  the 
type  stated  in  Theorem  6.2. 

Teixeira  [  I  ]  has  generalized  the  theory  of  first  order  bifurcations  in  the  plane  to  the 
case  where  the  vector  field  is  required  only  to  be  transversal  to  a  disk  except  at  a  finite  num 
ber  of  points.  Newhouse  and  Peixoto  [1J  have  shown  the  interesting  fact  that  one  can  al¬ 
ways  pass  from  one  Morse-Smale  system  to  another  by  a  one  parameter  family  of  vector 
fields  for  which  every  bifurcation  point  is  a  saddle-node.  This  result  is  valid  in  dimension  n. 
Generalizations  are  in  Newhouse  [2] . 

The  higher  order  bifurcation  points  are  difficult  to  classify  and  only  the  case  of  degree 
2  has  been  completely  resolved.  The  local  problem  near  a  generalized  saddle-node  has  been 
discussed  in  §4  and  near  a  generalized  focus  in  §5.  Leontovich  [1]  has  discussed  the  gen¬ 
eral  case  near  a  homoclinic  orbit.  For  further  results,  see  Andronov  et  a!.  [1] ,  Takens  [2] . 
An  example  will  be  given  later  of  a  higher  order  bifurcation  near  an  equilibrium  point. 
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7.  Two  dimensional  periodic  systems 

In  this  section,  our  objective  is  to  discuss  subharmonic  bifurcation,  bifurcation  of 
homoclinic  orbits  and  the  bifurcation  of  tori  for  periodic  two  dimensional  systems.  The 
equation  is  chosen  to  have  a  simple  form  in  order  to  minimize  the  notation.  Using  the 
same  methods,  generalizations  are  possible  (see.  for  example.  Chow  and  Hale  jl]). 

Suppose  /:  R  — *  R  is  a  Cr-function,  r  >  2.  /•':  R  — *  R  is  continuous  and  periodic  of 
period  one.  The  equation  to  be  considered  is 

(7.1 )  x  =  y,  y  =  ~J\x)  tty  +  kF(t) 

where  |  X.  p)  €  R~  varies  in  a  neighborhood  of  (0.  0).  The  unperturbed  equation  is 

(72)  x  =  y.  v  =  J\x). 

Suppose  equation  (7.2)  has  a  periodic  orbit  PA  of  least  period  A.  where  A  is  a  positive 
integer. 

Problem  SB.  Find  a  neighborhood  U  of  I*  and  a  neighborhood  ('  of  (  X,  p )  =  0  such 
that,  for  each  (X.  p)  €  V,  one  knows  exactly  the  number  of  k -peril >du  solutions  oj  (7 .1 ) 
which  belong  to  V  as  well  as  their  stability  properties 

A  periodic  solution  of  (7.1  )  which  is  k  times  the  period  of  the  forcing  function  /•(/)  is 
called  a  subharmonic  solution  of  order  k. 

A  solution  of  Problem  SB  requires  some  hypotheses.  If  1'^  ■(  pkU).  PkUU.  /UR 

where  (  pk{t).  pk(t ))  is  a  A-periodic  solution  of  (7.2).  then  the  derivative  of  tins  function  is 
a  A-periodic  solution  of  the  linear  variational  equation 

(7.3)  x  ~  y,  y~  (<’)/(  />^(/))A'l.v(  r. 

We  suppose  the  set  of  A  periodic  solutions  of  (7.3)  is  one  dimensional.  The  periodic  orbit 
\\  of  the  Hamiltonian  system  (7.2)  must  belong  to  a  smooth  family  of  periodic  orbits  l'w(a( 
of  period  oj(a)  where  a  is  a  real  parameter  in  (~e,  e)  with  co(0)  =  A.  One  can  show  the 
hypothesis  on  equation  (7.3)  is  equivalent  to  co'(0)  0.  This  hypothesis  is  generic  with  re¬ 

spect  to  the  class  of  vector  fields  f  6  C2( R,  R),  n  >  2,  with  the  Whitney  topology.  Also, 
if  is  sufficiently  close  to  a  homoclinic  orbit,  then  the  above  hypothesis  is  always  satisfied 
(See  Brunovsky  and  Chow  |1]). 

If  we  define  the  I -periodic  function  /i^(o)  by 
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(7.4) 


nk  =  /‘pJW*. 

V“)  =  1*  '  J/nW^  "  a)dt. 


then  we  can  prove  the  following  result  due  to  Hale  and  Taboas  [1]  (for  an  abstract  version, 
see  Hale  (9|). 

Thkorem  7.1 .  Let  am ,  aM  be  respectively  the  value  of  a  in  [0,1)  for  which  hk[a) 
is  a  minimum,  maximum  and  suppose  hk(am)  >  0,  h'k(aM)  <  0.  Then  there  are  neighbor¬ 
hoods  U  of  P*.  V  of  ( X,  p)  =  (0,  0)  and  two  curves  C£ ,  Cj^  i  V  containing  (0,  0),  tangent 
respectively  to  the  lines  X  =  hk(am)p,  X  =  hk(aM)p  at  (0,  0),  dividing  K\{0}  into  two  dis¬ 
joint  open  sets  S*,  5*  such  that  equation  (7.1)  has  no  subharmonic  solution  in  U  for 
p)  £  5*  and  at  least  2k  for  (X,  p)  in  S*. 

We  remark  that  the  hypothesis  on  h  is  generic  in  f  6  C(R,  R).  f\t  +  1)  =  F(t). 

In  the  proof  of  this  theorem,  it  is  shown  that  2k  subharmonic  solutions  appear  as  X/p 
decreases  through  a  maximum  (increases  through  a  minimum)  of  hkfa).  Thus,  if  the  func¬ 
tion  hk(a)  has  no  other  maxima  or  minima  in  [0,  I),  then  the  conclusion  in  Theorem  7.1 
can  be  strengthened  to  say  there  are  exactly  2k  subharmonic  solutions  for  (X.  p)  £  S*. 

Only  the  idea  of  the  proof  will  be  given.  There  is  a  neighborhood  f/of  T*.  8  >  0. 
such  that  the  mapping 

x  -  pia)  -  ap(a),  v  =  p(a)  +  ap(a) 


is  a  one-to-one  C’  transformation  from  U  onto  [0,  k)  x  {|a|  <  6} .  Thus,  it  is  sufficient  to 
consider  only  A: -periodic  solutions  of  the  form 

x(t)  =  p(t  +  a)  +  :(t  +  a),  y(t )  =  x(t). 


r(a)  =  -apfa),  z(a)  =  ap(a) 

with  |n|  <  5,  a  £  [0,  k).  If  this  change  of  variables  is  performed,  then  the  method  of 
Uapunov-Sclrmidt  can  be  applied  to  obtain  the  usual  bifurcation  function  B( a.  X.  p).  This 
function  has  the  form 

B(a,  X.  p)  =  —  Xr?fc  +  prikhk{a)  +  0(|r)|  +  |p|)J  as  |X|,  |p|  --*■  0. 

The  remaining  argument  is  a  careful  analysis  of  the  zeros  of  this  function  (see  Hale  and 
Taboas  [1 1 ). 

For  a  given  o,  X,  p  for  which  B(a,  X.  p)  =  0,  one  obiains  a  k -periodic  solution  which 
is  close  to  (p(/  +  a).  p(t  +  a)).  The  constant  a  corresponds  to  a  phase  shift  and  is  approxi 
matcly  determined  from  the  relation  hk(a)  =  X/p.  If  hk(a)  0  for  a  F  am,  aM.  then  one 
can  associate  a  unique  a  =  a*(X/p)  with  X/p.  As  X,  p  — »  0.  this  function  a*(X/p)  will  ap¬ 
proach  a  limit  if  and  only  if  X/p  approaches  a  constant.  This  implies  that  the  initial  values 
of  the  corresponding  X-periodic  solution  (which  are  approximately  p(a  *(  X/p),  /Ha  *(  X/p))) 
will  converge  to  a  point  on  P*  as  X,  p  — *  0  if  and  only  if  X/p  approaches  a  constant  as 
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X,  p  — *•  0.  This  points  out  one  of  the  disadvantages  of  treating  the  original  two  parameter 
problems  in  (X,  p)  as  a  one  parameter  problem  along  a  ray  (X.  mX),  m  fixed,  in  parameter 
space. 

Using  the  results  of  §4,  if  hk(ct)  has  a  finite  number  of  extremal  values  in  |0.  I ),  one 
can  show  that  there  are  at  least  2k  hyperbolic  subharmonics  of  order  k  with  k  being  saddles 
and  k  being  either  nodes  or  foci  (see,  for  example,  Chow  and  Hale  jl  | ). 

A  specific  example  of  equation  (7.1 )  is 

(7.5)  x  -  x  +  x1  ~  -\x  +  nF(t). 


Y 


,»j 


x 


Figurf  10. 

For  (X,  p)  =  (0,  0),  the  phase  portrait  is  shown  in  Figure  10  with  any  neighborhood  of  the 
homoclinic  orbit  containing  periodic  orbits  Fw((|j  with  period  iota),  where  (a,  0)  is  the 
initial  value  of  the  solution  at  t  ~  0.  One  can  show  that  to  (a)  — »•  °°  as  a  a„  monotom- 
catly.  Thus,  there  is  a  sequence  aK  —*  such  that  t o(ok)  =  k  for  k  >  k0  (=  7).  Tire  pre¬ 
vious  results  can  be  applied  to  each  such  Ffc  for  the  subharmonic  bifurcation  curvcsC^1 ,  C J,. 
This  suggests  that  the  behavior  of  solutions  of  (7.1)  near  r„  could  be  very  complicated.  We 
show  that  this  is  actually  the  case. 

Suppose  /(0)  =  0,  f'(0)  <  0  and  there  is  a  homoclinic  orbit  through  (0.  0). 

P„  =  {(p«,(f),  p  „(!)),  l  €  R}  U  {(0,  0)} ,  where  p„(f)  — 0  as  r  — ♦  ±  °°  is  a  solution  of 
(7.2).  Since  (0,  0)  is  a  saddle  point,  there  are  a  neighborhood  V  of  (X,  p)  -  (0,  0)  and  a 
I -periodic  solution  (0(/,  X,p),  0(/,  X.p))  of  (7.1)  for  (X,p)  e  V,  0(1.  0,  0)  =  0.  Further¬ 
more.  the  trajectory  yK  u  C  R  x  R2  of  this  solution  has  stable  manifold  Sk  and  unstable 
manifold  t/KjJ  with  1 -periodic  cross-sections 

=  «*'  3’)  e  VMf-  UxJ<)  ~  {<-».  >•>:  <*.  x.  .)■)  e  UK  (l  | . 

Let  n  be  the  period  one  map  of  equation  (7.1).  Then  pk  -  (0(0,  X.  p),  0(0.  X.  p))  is  a 
fixed  point  of  rr  and  M(0),  Uku(0)  arc  the  stable  and  unstable  manifolds  of  p. 

Problem  HB.  Find  a  neighborhood  U  of  and  a  neighborhood  V  of  ( X,  p)  -  (0.  0) 
such  that,  for  each  ( X.  p)  C  V,  one  knows  whether  or  not  there  is  a  point  q  6  V  which  is 
transverse  homoclinic  to  pk 


k 
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As  noted  in  §3,  the  existence  of  a  point  which  is  transverse  homoclinic  to  pK  im¬ 
plies  a  type  of  random  motion  occurs  in  a  neighborhood  of  Our  primary  concern  here 
is  to  show  how  the  transverse  homoclinic  point  occurs  through  bifurcation  and  the  connec¬ 
tion  with  the  previous  problem  of  subharmonic  bifurcation. 

Define  the  1 -periodic  function  /i_(a)  by 

(7.6)  pi(t)dt.  h„(a)  =  1  J p~(t)h\t  -  a)  Jr 

and  assume  that 

(7.7)  h"(am)>0.  h'\aM)  <  0 

where  am.  aM  are  the  values  of  a  at  which  h  assumes  respectively  its  minimum,  maximum. 
We  have  the  following  result  of  Chow,  Hale  and  Mallet -Paret  |2| 

Tm  oKi  M  7.2.  There  are  neighborhoods  V  of  r_ .  V  of  ( X.  p)  =  (0.  0)  and  two 
curves  C“.  £  V containing  (0,  0)  tangent  respectively  to  the  lines  X  =  h„(am  )p.  X  = 

X’„(a M)p  at  (0,  0).  dividing  F\{(0.  0)}  into  two  disjoint  open  sets  S“.  Sf  such  that  equa¬ 
tion  (7.1)  has  no  homoclinic  points  in  S“  and  has  transverse  homoclinic  points  in  Sf 

In  any  neighborhood  of  any  point  (X.  p)  in  the  bifurcation  set  C“  U  C^.  there  arc 
infinitely  many  subharmonic  bifurcations  of  the  type  in  Theorem  7.1 . 

Finally,  if  h'„(a)  has  only  a  finite  number  of  zeros  in  (0.  I ).  then  the  period  one  map 
rr  has  infinitely  many  hyperbolic  saddles  and  nodes  (or  foci)  in  U  for  (\.  p)  €  S~ 

Note  the  conditions  on  /i„  are  generic  in  F. 

The  proof  of  the  first  part  of  this  theorem  follows  along  the  same  lines  as  the  proof 
of  Theorem  7.1  for  subharmonic  solutions.  One  must  replace  the  Fredholm  alternative  for 
periodic  solutions  by  a  Fredholm  alternative  for  solutions  bounded  on  R.  This  gives  rise  to 
the  function  h^ja)  in  (7.6),  a  function  whose  importance  in  this  problem  was  emphasized 
by  Mel'nikov  (I).  To  prove  the  part  on  subharmonic  solutions,  one  shows  that  condition 
(7.7)  implies  the  corresponding  conditions  in  Theorem  7.1  are  satisfied  for  each  hk(cs)  for 
which  l'k  £  U.  See  Chow,  Hale  and  Mallet -Paret  [1]  for  details. 

Holmes  [1]  has  also  proved  the  first  part  of  Theorem  7.2  dealing  with  the  existence 
of  homoclinic  points.  Closely  related  results  for  Hamiltonian  systems  have  been  obtained  by 
Churchill,  Pccetli  and  Rod  [1] .  Holmes  and  Marsdcn  |l  |  have  extended  the  first  part  of 
this  theorem  to  certain  types  of  partial  differential  equations  and  have  applied  the  results  to 
the  equations  for  a  beam.  In  this  case,  technical  difficulties  primarily  arise  from  the  fact 
that  the  semigroup  generated  by  the  linear  part  of  the  unperturbed  system  has  spectrum 
which  can  include  the  complete  unit  circle.  A  new  phenomenon  arises  because  the  spectrum 
contains  one.  This  can  also  occur  in  the  finite  dimensional  case.  For  example,  consider  the 
equations  in  R3, 

x~x  +  x2=-Xx  +  uF(t)  +  H(x.  v),  y  =  -  \v  +  pC,(t)  +  /.(x,  .»•) 
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where  F,  G  have  period  1 ,  H,  L  vanish  together  with  their  first  derivatives  at  x  =  0,  y  =  0 
and  H(x,  0)  =  0(ix|3)  as  |x|  — ►  0,  IXx,  0)  =  0.  For  X  =  p  =  0,  this  equation  has  a  homo¬ 
clinic  orbit  in  the  (x,  x)-plane.  The  equilibrium  point  x  =  x  =  0,  y  -  0  is  not  a  saddle  since 
one  eigenvalue  is  zero.  To  discuss  homoclinic  orbits  for  the  perturbed  equation,  one  must 
find  first  a  periodic  solution  of  period  one.  The  difficulties  involved  are  the  same  as  the 
ones  encountered  in  finding  a  1-periodic  solution  which  approaches  zero  as  X,  p  — <•  0  of  an 
equation  of  the  form  y  =  -\y  +  pG(f)  +  M(y)  where  M(y)  =  Odyl2)  as  |>>|  — *  0.  Given 
the  first  nonvanishing  coefficient  of  Mi  y)  and  assuming  that  G(t)dt  #  0.  this  puts  re¬ 
strictions  on  the  parameters  X,  p.  For  example,  if  M(y )  =  0yk  +  (X\y\k  +  1 ),  0  =*  0,  then 
k  =  2  implies  p/X3  approaches  a  constant  as  X  — ►  0.  If  k  -  3,  then  p/X2  approaches  a 
constant  as  X  — *  0.  Once  the  l -periodic  solution  has  been  obtained,  it  is  not  too  difficult 
to  adapt  the  previous  ideas  to  obtain  a  transverse  homoclinic  point  in  some  region  of  the 
parameter  space. 

It  often  occurs  that  one  is  interested  in  studying  problems  where  the  orbit  is  not  a 
homoclimc  orbit,  but  is  the  closure  of  an  orbit  connecting  two  distinct  equilibrium  points. 
To  obtain  a  generalization  of  Theorem  7.2  for  this  case,  we  give  a  procedure"  motivating  the 
appropriate  formula  for  the  corresponding  function  A„( a).  Consider  the  equation  (7.5)  with 
X  =  0.  Suppose  there  is  a  homoclinic  orbit  for  X  small.  It  must  be  close  to  /?„.  If  it  has 
the  form  p_  +  pp  +  o(|pl)  then  q  satisfies  the  equation  q  +  /'(p_)q  =  F(t).  Multiply 
by  P«.  integrate  from  0  to  «>,  and  integrate  the  term  pmq  by  parts  twice  The  result  is 

-p.(0)<?(0)  =  f~  Flt)p„(r)dr 

Doing  the  same  thing  but  integrating  from  -  •»  to  zero,  one  has  p_(0)<j(0)  =  S° ^h\t)paa{t)(Jt. 
Since  p„(0)  #  0,  this  implies  [~mF\i)pmU)dt  =  0,  which  is  the  same  as  /iM(0)  =  0. 

Let  us  turn  our  attention  now  to  the  bifurcation  of  tori  Suppose  G:  R2  — ►  R2, 

F:  R  x  R2  x  R*  — ►  R2  are  C-functions,  r  >  2,  F\l,  i\  e)  periodic  in  t,  F(t.  v,  0)  =  0.  and 
consider  the  equation 

(7.8)  y  =  G(y)  +  F\t.  y,  e). 

For  e  =  0,  suppose  this  equation  has  an  io-periodic  orbit  T.  In  a  neighborhood  of  P  one 
can  introduce  a  moving  orthonormal  coordinate  system  (p,  0)  which  is  w-periodic  in  0.  The 
new  equations  for  (p,  0 )  have  the  form 

(7.9)  0  =  1  +  9(f,  0,  p,  e),  p  =  R(r,0,p,e) 

where  8,  R  are  periodic  in  i,  6,  8 (/,  6,  0,  0)  =  0,  R{t,  0,  p,  0)  =  ap  +  Ofipl2)  as  Ipl  0. 
If  a  #  0,  that  is,  the  orbit  T  of  (7.8)  for  e  =  0  is  hyperbolic,  then  the  classical  results  on  the 
theory  of  integral  manifolds  (see,  for  example.  Hale  [1])  imply  the  existence  of  a  function 
p*(f,  6,  e),  p*( t,  0,  0)  =  0  defined  for  all  t,  e  and  periodic  in  both  variables  such  that 
(0(f),  p*(f,  0(f),  e))  is  a  solution  of  (7.9)  for  every  solution  0(f)  of 

0  =  1  +  8(f,  0.  p*(f,  0,  e),  e). 
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For  equation  (7.8),  this  implies  the  existence  of  an  invariant  cylinder  in  R  x  R2  with  peri¬ 
odic  cross-section  near  the  cylinder  R  x  T.  Since  the  cross-section  is  periodic,  the  flow  on 
the  cylinder  is  equivalent  to  the  existence  of  a  flow  on  a  torus.  If  the  original  orbit  T  is  not 
hyperbolic,  one  must  determine  more  explicitly  how  the  nonlinear  terms  influence  the  flow. 
We  return  to  this  problem  later. 

Suppose  H:  R3  x  R*  — *  R3  and  consider  the  equation 

(7.10)  i  =  H(z.  e). 

Suppose  for  e  =  0  this  equation  has  a  smooth  invariant  torus  T2  on  which  the  flow  is  parallel; 
that  is,  cither  the  rotation  number  is  irrational  or  every  orbit  is  periodic.  It  is  possible  to  intro¬ 
duce  coordinates  (p,  0,  f )  in  a  neighborhood  of  T 2  such  that  equation  (7.10)  is  equivalent  to 

(7.11)  f  =  1  +  Z(f.  0,  p,  e).  6=  1  +0(f,  d.p.e),  p  =  tf(f,  0.  p,  e) 

where  all  functions  are  periodic  in  f,  8  and  vanish  for  (p,  e)  =  (0,  0).  Any  integral  mani¬ 
fold  of  equation  (7.1 1)  which  can  be  represented  by  a  function  p*(f,  6.  e),  periodic  in  f,  8 
with  the  same  period  as  the  vector  field,  will  correspond  to  an  invariant  torus  of  (7.10). 

If  the  original  torus  T 2  is  hyperbolic,  then  equation  (7.10)  will  have  a  unique  invari¬ 
ant  torus  in  a  neighborhood  of  T2.  We  will  discuss  some  aspects  of  the  problem  when  T2 
is  not  hyperbolic.  Note  that  equation  (7.9)  can  be  considered  as  a  special  case  of  equation 

(7.1 1)  by  putting  t  -  f. 

To  be  specific,  suppose  e  =  (X,  p)  6  R2  and 

(7.12)  K(f,  O.p.  A.p)  =  p2  +  p/(0,f)  +  Xgtd,r)  +  ^(f.  0.P- X.p) 

where  R  =  0(|p|3  +  (|p|  +  |A|)(lp|  +  |A|  +  |p|))  as  |A|.  |p|.  |p|  —  0.  For  (X,  p)  =  (0.  0). 
equation  (7.1 1)  has  the  invariant  torus  p  =  0  which  is  asymptotically  stable  from  one  side 
and  unstable  from  the  other.  One  would  suspect  that  there  are  no  invariant  tori  for  some 
(X,  p)  near  zero  and  two  invariant  tori  for  other  (X,  p);  that  is,  there  is  a  bifurcation  of  tori 
similar  to  the  saddle-node  type  bifurcation  for  equilibrium  points.  We  will  confirm  some  of 
these  suspicions  and  also  point  out  that  other  phenomena  can  occur. 

Let  us  first  consider  the  one  parameter  problem  with  X  =  0.  If 

,7  I3)  fo  =  j:  fQr/(8  +  t,  f  +  f)dt 

is  independent  of  ( 0 ,  f )  (this  is  a  nonresonance  condition),  then  the  scaling  p  — *  \J\p\p 
can  be  justified.  Using  the  standard  method  of  averaging  (see  Hale  [!|,  Diliberto  [ I J ).  there 
is  a  transformation  of  variables  which  takes  (7.1 1)  into  an  equivalent  equation 

f  =  1  +  f4lpl'/2),  0  =  1  +  0dpi1'2), 

^  p  =  lpi,/2(p2  +  /oSgn  p)  +  o(|p|,/;) 


as  |p|  — *  0.  Integral  manifold  theory  implies  there  are  no  invariant  tori  if  /0  sgn  p  >  0  and 
two  hyperbolic  invariant  tori  if  /0  sgn  p  <  0. 
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If g0  =  liirij. /org(0  +  /,{■+  t)dt  is  independent  of  0,  f  and  X  =  7/},  then  one 
can  apply  the  same  result  to  the  complete  equation  (7.1 1)  with  the  only  change  being  that 
f0  sgn  g  is  replaced  by  (/0  +  -y£0)sgn  /a.  If  this  quantity  is  >  0,  there  are  no  invariant  tori 
and  if  this  quantity  is  <  0,  there  are  two  invariant  tori. 

From  the  above  discussion,  we  have  obtained  the  following  result.  If  f0  +  0,  g0  +  0, 
then,  for  any  y  for  which  f0  +  gg0  4=-  0,  there  is  a  gQ  =  g0(y.  f0,  g0)  such  that  on  the 
line  segment  Ly  =  {(X,  g):  X  =  yg,  0  <  Igj  <  g0),  there  is  no  invariant  torus  of  (7.1 1 ) 
in  a  neighborhood  of  zero  if  (/0  +  y?0)sgn  g  >  0  and  there  are  two  hyperbolic  tori  if 
U'o  +  r^o)s8n  M  >  0. 

This  result  implies  a  complete  solution  of  the  problem  of  bifurcation  of  tori  except  in 
a  sector  SyQ  containing  the  line  segment  /.  ,  y0  =  ~f0/g0  (see  Figure  1 1).  This  sector 

must  be  tangent  to  Lyo  at  (0,  0).  If  one  assumes  the  vector  field  is  C"  with  the  coeffi¬ 
cients  in  the  Taylor  series  in  p  being  trigonometric  polynomials  and  the  ratio  of  the  periods 
in  6,  f  are  irrational,  then  it  is  possible  to  show  that  the  excluded  sector  Syo  is  tangent  to 
Lyo  at  (0,  0)  to  infinite  order. 


/ 


Figurf  II, 

What  happens  in  SyQ  is  very  complicated.  In  fact,  when  the  two  tori  come  too  close  to¬ 
gether,  then  genetically  the  rotation  number  is  rational  and  hyperbolic  periodic  orbits  ap¬ 
pear.  The  two  invariant  sets  then  can  meet  at  these  hyperbolic  periodic  orbits  and  then 
change  topological  structure.  An  illustrative  example  of  the  cross-sections  of  these  sets  is 
shown  in  Figure  12  as  the  generic  bifurcation  takes  place.  The  interior  set  is  exponentially 
stable  in  Figure  12a  and  the  exterior  one  is  unstable  exponentially.  The  points  correspond¬ 
ing  to  the  periodic  orbits  are  hyperbolic-all  being  saudle  points  except  lor  the  nodes  .4.  /?, 
C,  D.  In  Figure  1 2b,  the  saddles  remain  and  thepoints  A,  B  are  saddle-nodes.  In  Figure  12c. 
the  invariant  sets  that  were  similar  to  tori  have  disappeared  leaving  only  hyperbolic  orbits. 


8.  ifigher  order  bifurcation  near  equilibrium 

In  86  the  vector  Helds  in  the  plane  corresponding  to  bifurcation  points  of  order  one 
were  completely  characterized .  Higher  order  bifurcation  in  any  dimension  were  discussed  m 
SS4  5  for  the  case  of  a  generalized  saddle-node  or  focus.  In  this  section,  we  give  two  illus- 

plele  analysis  co«.p.s  mhicli  ».  global  in  naiine-m  *«<*■;  »**" 

l,  homoli.  orbits  and  invade  ioii.  On.  «ec.o,  r.aid  is  «o  d— al 

matdx  of  the  lino.,  app.osim.tion  n...  »ro  has  bolh  tools  z.ro  but  with  nonsimplt  «■ 

„  i  j  mi  The  second  illustration  is  a  three 

mentary  divisors  and  is  based  on  Howard  and  K  pe  l)'  having  two 

dimensional  vector  field  with  the  matrix  of  the  linear  approximation  near  zero  having 
“aginary  and  one  zero  root.  The  relevant  literature  for  this  latter  example  is  Ung- 

ford  [1  ] ,  Guckenheimer  [3] ,  Chow  and  Hale  ll]  ■ 

Consider  the  equation 

(8.1)  x  =  y,  y  -  +  ezy  +  ax2  +  fry 

where  «  <  0.0  >  0  are  fixed  constants  and  e  =  («,.€,)  «•  «al,  varying  in  a  neighborhood 
°f  “lowing  the  behavior  of  (8.0  will  also  determine  the  behavior  for  the  caseo^  0^ 


-r,  c. 


•  -e-> 


Ain  since  the  change  of  variables  t  .  ■ 

above  case  Also,  the  form  (8.1)  is  chosen  for  simplicity  and  the  qualitative  structure 
will  be  unchanged  if  perturbations  (X,  Y)  in  the  vector  field  are  made  which  sa  i 
y  -  Offlxl  +  l  vl)J  +  UKlxl  +  |y|)),  Y  =  0«\x\  +  W)  0  +  lel»- 

'  Tlbit « *. <—•« »« ■*  *» <8 *>  ”  *  n”8h- 

borhood  of  (x,  y)  =  (0, 0)  for  e  varying  in  a  neighborhood  of  zero. 

Tw°. w.  Lid.,  ,h.  case  . ,  >  0-  Th.  ihimduciiph  »f  •»  s»W  »nabi.s 


(8.2) 


e,  =  6J,  e2  =  kS1.  8  >0, 

tv-S-'f,  x  i-62|ar,x.  y  t-53la|"V 


2  +  u6y  +  Syxy 


| 

I 


i 


leads  to  the  new  equations 

(8.3)  *  “  y>  y  =  *~x 

where  y  =  dial  1  ■ 
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For  6=0,  equation  (8.3)  is  conservative  with  the  first  integral  V(x,  y)  =  y2/ 2  - 
x2I2  +  x3/3.  The  equilibrium  point  (0,  0)  is  a  saddle  point  with  a  homoclinic  orbit  r„ 
through  it  while  the  equilibrium  point  (1, 0)  is  a  center.  The  first  step  is  to  analyze  the 
periodic  orbits  of  equation  (8.3). 

Lemma  8.1 .  Every  periodic  orbit  of  equation  (8.3)  must  intersect  the  segment  (0,  1) 
x  {0}  in  the  (x,  yY plane.  There  are  a  continuous  positive  function  6 0(Z>),  b  €  (0,  1),  and 
a  continuously  differentiable  function  p  *(b,  6),  b  £  (0,  1),  |5|  <  &0(b)  such  that  equation 
(8.3)  has  a  periodic  orbit  if  and  only  if  p  =  p*(b,  8).  Furthermore,  dp.  *(b,  0  )/db  <  0, 
p\b,  0)  — *  -y  as  b  — *  1 ,  p \b,  0)  — *  ~6y/7  as  b  — ►  0.  Finally,  if  p  =  p*(b,  8)  fora 
fixed  b  e  (0,  1)  and  |5|  <  50(6),  then  the  periodic  orbit  through  ( b ,  0)  is  the  only  one  cor¬ 
responding  to  this  p,  8. 

Only  a  sketch  of  the  proof  is  given  (see  Chow  and  Hale  [1]  for  details  based  on  Carr 
(1)).  The  first  statement  is  obvious  since  the  index  of  a  periodic  orbit  is  one.  If  b  €  (0,  1) 
is  on  a  periodic  orbit  P,  then  /r  Vdt  =  0.  Using  the  fact  that  the  orbit  T  must  be  symmetric 
with  respect  to  the  x-axis,  this  implies  that 

P  y2dx  +  y  f*(b)  xy2  dx  =  0 

where  ( c(b ),  0)  €  T,  c(b)  >  1 .  This  defines  p*(b,  8).  The  other  assertions  require  a  number 
of  computations  which  will  not  be  given. 

Let  r„  =  {(p„(/),  p„(f)),  t  €  R}  U  {(0,  0)}  where  p„  is  a  solution  of  (8.3)  for 
e  =  0,  p„(r)  — *  0  as  t  —> ►  ±  °o.  Using  the  type  of  analysis  in  §7,  one  can  find  a  unique 
curve  in  the  (p,  6)-plane  for  5  small  along  which  equation  (8.3)  has  a  homoclinic  orbit.  This 
gives  a  curve  C„  in  the  e  =  (et ,  e2)-plane  defined  for  0  <  |e|  <  e0  parametrically  by  a  func¬ 
tion  e2  =  c(e{^2)e,  where  c(0)  =  -yv,  v  =  J~^q0,ql,lf“a.q2 ■  On  C„,  the  homoclinic  orbit 
is  asymptotically  stable  from  §6.  Using  the  formulas  which  define  the  homoclinic  orbit  from 
Chow,  Hale  and  Mallet-Paret  [2] ,  one  shows  that  there  is  a  unique  periodic  orbit  to  the  left 
of  C„.  To  the  right  of  C„.  every  solution  of  (8.3)  leaves  a  neighborhood  of  T„. 

To  the  left  of  C„  and  close  to  C„  where  the  periodic  orbit  exists,  this  orbit  passes 
through  ( b ,  0)  where  b  is  close  to  zero.  Lemma  2.1  implies  that  this  orbit  continues  to  exist 
and  is  the  unique  periodic  orbit  on  the  ray  e2  =  p*(b,  6)e, .  Thus,  as  long  as  b  remains  in 
an  interval  (0,  I  -  d] ,  d  E  (0,  1 ),  we  have  the  existence  on  a  uniform  6-interval  0  <  |6|  < 
80(d).  To  obtain  a  uniform  5-intcrva)  for  b  E  (0,  1).  we  analyze  carefully  the  neighborhood 
of  the  equilibrium  point  (1,0).  The  matrix  of  the  linear  variational  equation  near  the  equil¬ 
ibrium  point  (1,0)  has  complex  eigenvalues  with  real  parts  8(p  +  y)  and  are  ±i  at  6  =  0. 

The  bifurcation  function  G(a,  6,  p)  for  periodic  orbits  obtained  as  in  §5  is  an  analytic  func¬ 
tion  of  (a,  5,  p)  near  (a,  0,  -y).  Also,  G(a,  0,  p)  =  0  for  all  ( a ,  p)  since  equation  (8.3)  for 
6  =  0  has  a  first  integral  (Corollary  5.4).  Since  (7(0,  5,  p)  =  0  for  all  6,  p.  this  implies 
that  the  bifurcation  equation  G(a,  6,  p)  =  0  can  be  replaced  by  the  equivalent  equation 
H{a,  6,  p)  =  0,  where  H{a,  6,  jj)  =  G(a,  6,  p)/a8.  Also,  H(0,  0,  p)  =  p  +  6  +  (terms 
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vanishing  for  (a,  6)  =  (0, 0)).  The  Implicit  Function  Theorem  implies  there  is  a  unique  func¬ 
tion  p*(a,  6)  such  that  H(a,  6,  8))  =  0  for  la]  <  a0,  |6|  <  60.  This  gives  the  addi¬ 

tional  uniformity  in  5  for  b  near  1 . 

Combining  all  of  this  information,  the  half-plane  t,  >  0  can  be  divided  into  four  sec¬ 
tors  as  shown  in  Figure  13a  for  which  the  corresponding  flows  are  given  in  Figures  1 3b,  c.  d. 


y> 


l 

Figurk  13b.  Sector  1. 


Figure  13c.  Sector  2.  Figure  13d.  Sector  3. 


Arnol'd  [1,  p.  255]  has  considered  the  same  nonlinearities  as  in  (8.1)  but  with  the 
perturbation  term  e2  +  e,x. 

Carr  [1]  has  given  a  complete  discussion  of  the  bifurcations  in  the  more  complicated 
equation  x  -  y.  y  =  e,x  +  e2  y  +  ox3  +  0x2y. 

The  next  example  is  concerned  with  a  differential  equation  in  R3  for  which  the  linear 
variational  equation  near  the  equilibrium  point  zero  has  two  purely  imaginary  eigenvalues 
and  one  zero  eigenvalue.  To  simplify  the  situation,  it  will  be  assumed  also  that  a  certain 
type  of  symmetry  prevails.  More  specifically,  consider  the  equation 

(8.4)  x  =  A(\)x  +  /\x.  y),  y  -  0y  +  g(x,  y) 

where  X,  0  are  small  real  parameters,  x  €  R:,  y  €  K,f,  g  are  C4 -functions. 
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f(x,  y)  =  0(1x10x1  +  |>-|)),  g(x,  y)  -  0(1*1  +  |.y|)2 

as  |x|,  lyl  — *  0.  The  hypothesis  on  /  implies  that  the  symmetry  condition 

(8.6)  /(0,  jy)  =  0 

is  satisfied. 

Since  /( 0,  y)  =  0,  it  is  legitimate  to  introduce  polar  coordinates  for  x  =  (x,,  x2 )  as 
x,  =  p  cos  0,  x2  =  -p  sin  8.  If  this  is  done  and  t  is  replaced  by  8.  one  obtains  the  equa¬ 
tions 

(8.7)  p  =  Xp  +  R(0,  p,  y,  a,  0),  y  =  0y  +  Y(8,  p,  y.  a,  0) 

where  R  =  0(|p|(ip|  +  |jp|)),  Y  =  0((|p|  +  lj-1)2)  as  p.  —  0. 

The  problem  is  to  determine  the  behavior  of  the  solutions  of  (8.7)  in  a  neighborhood 
of  (p.  y)  =  (0,  0)  for  (X,  0)  in  a  neighborhood  of  (0,  0).  To  discuss  (8.7).  we  suppose  that 
an  application  of  the  theory  of  normal  forms  to  (8.7)  yields  the  system 

(8.8)  p  =  p(\  +  ay  +  bp2)  +  0(p(\p\  +  l.vl)J),  v  =  0y  +  cy 2  +  dp 2  +  0((lp |  +  |v|)4) 

as  p.  i'  — *  0.  The  generic  situation  is  to  have  a,  b,  c,  d  nonzero  constants. 

To  simplify  the  computations  and  also  to  consider  the  most  interesting  case,  we  sup¬ 
pose  a  =  2,  b  -  l,c  -  d  =  -  l,so  that  the  truncated  form  of  (8.8)  using  only  the  lowest 
order  terms  is 

(8.6)  P  =  p(X  +  2 y  +  p2),  v  =  0v  -  y2  -  p2. 

The  first  objective  is  to  discuss  the  behavior  of  the  solutions  of  (8.9)  as  a  function  of  o.  0. 

It  is  convenient  to  introduce  the  change  of  variables  y  0/2  +  y,  X  +  /?*-•  a,  to  obtain  the 
more  symmetric  form 

(8.10)  p  =  p(a  +  2.1’  +  pl).  y  =  /J2/4  -  y2  -  p2. 

If  we  perform  the  scalings 

(8.11)  p  hf  ep,  y  *-*  ev,  0  •-»  e.  a  / 1— ►  e _  1  r 

the  new  equations  become 

(8.12)  p  =  2pv  +  ap  +  ep3.  y  =  1/4  -  y2  -  p2. 

Equations  (8.12)  must  be  discussed  for  all  p  >  0,  y  €  R,  a  e  R  and  e  in  a  neighborhood  of 

c  =  0. 

for  q  =  0,  e  =  0,  the  function 
(8  13)  V(p.  y)  =  pl*  ~  py1  ~  P3H 

is  a  first  integral  of  (8.12).  Differentiating  this  function  along  the  solutions  of  (8.12).  it  is 
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not  difficult  to  see  that  no  periodic  orbits  can  exist  unless  a  =  0  when  e  =  0.  This  remark 
makes  much  of  the  discussion  of  (8.12)  very  simple  since  it  is  only  necessary  to  determine 
the  structure  of  the  equilibrium  points  away  from  the  critical  value  a  =  0,  e  =  0. 

The  equilibrium  points  of  (8.12)  are  p  =  0,  y  =  ±  1/2,  for  all  a,  e  and  the  point  y  = 
-all,  p2  =  (1  -  a2)/4  for  e  =  0,  a2  <  1.  For  a2  =*  1  the  points  p  =  0,  y  =  ±  1/2  are  hy¬ 
perbolic  saddles  or  nodes.  Further  analysis  of  the  other  equilibrium  point  for  a2  <  1  shows 
that  the  original  (X,  0)-plane,  \  =  0a  -  0, 0  =  f,  can  be  divided  into  sectors  as  shown  in 
Figure  14a  with  the  corresponding  phase  protrait  for  0  >  0  given  in  Figure  IS.  The  situation 
is  similar  for  0  >  0.  There  is  a  saddle-node  bifurcation  on  each  of  the  solid  curves  in 
Figure  14a.  The  analysis  in  a  neighborhood  of  the  dotted  line  X  =  -0  requires  further  study 
and  will  now  be  given. 


Figure  14a. 


Figure  14b. 
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Region  1 


Region  2 


Region  2' 


Region  3 


e 

i 
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Region  4 
Figure  15. 
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For  a  =  0,  e  =  0  in  (8.12),  the  phase  portrait  is  determined  by  the  first  integral 
V(p,y)  in  (8.13)  and  is  shown  in  Figure  16.  There  is  a  heteroclinic  orbit  connecting  the 
two  saddle  points.  Let  (p0(t),  y0(t)),  p0(/)  >  0,  be  the  solution  of  (8.12)  for  a  =  e  =  0 
such  that  p0(r)  — *  0  as  t  — * ►  ±  <*»,  >>0(f)  •— *  -  V4  as  t  — >  -<*,  y0(r)  -1 *  Vi  as  t  — *  °».  Using 
the  remarks  after  Theorem  7.2,  one  can  obtain  a  function  G(a,  e )  for  |a|,  |e|  <  5  such  that 


DYNAMIC'  BIFURCATION  THKOK Y 


51 


1 


equation  (8.12)  has  a  heteroclinic  orbit  if  and  only  if  G(a,  e)  =  0.  Furthermore,  the  func¬ 
tion  G  satisfies 


G’(a,  e)  =  a  j +  e  f_~P%(0dt  +  0((\a\  +  |e|)2) 

as  a,  e  — *  0.  The  Implicit  Function  Theorem  implies  there  is  a  unique  solution  a  =  a  *(e) 
of  this  equation  for  e  sufficiently  small,  a*(0)  =  0  and 

V-  4£S>.  _  J-  p;(„a,//-  „•««.  =  -i 

This  gives  a  curve  X  =  X*(0)  in  the  original  parameter  space  (X,  (3)  where  there  is  a  hetero¬ 
clinic  orbit  for  equation  (8.10).  This  curve  is  given  approximately  by  X  *(/3)  =  -0  -  f>02. 

Near  the  line  X  =  -0  (or  a  =  0),  there  is  also  a  Hopf  bifurcation.  Near  a  =  0,  e  -  0. 
equation  (8.12)  has  an  equilibrium  solution  given  approximately  by  y0  =  -(a  +  e/4)/2. 
p0  =  1/2  +  y0.  Analyzing  the  stability  properties  of  this  solution,  we  see  that  it  has  eigen¬ 
values  on  the  imaginary  axis  along  a  curve  given  approximately  by  a  =  -3e/4  or  X  =  -0  - 
3f ]2/4.  One  can  actually  show  there  is  a  Hopf  bifurcation  along  this  curve.  Also,  one  can 
show  there  is  a  unique  periodic  orbit  between  this  curve  and  the  one  defining  the  horno- 
clinic  orbit.  This  implies  the  complete  bifurcation  diagram  near  the  line  X  =  -0  is  the  one 
shown  in  Figure  14b.  The  flow  in  region  2'  is  shown  in  Figure  15.  The  flow  on  the  homo¬ 
clinic  curve  is  shown  in  Figure  16.  For  the  details  of  the  above  computations  see  Chow  and 
Hale  (1). 

It  remains  to  relate  the  solutions  of  the  approximate  equations  (8.12)  to  the  complete 
equations  (8.8).  Since  the  perturbation  terms  are  periodic  in  t,  it  is  not  difficult  to  show 
that  the  equilibrium  points  of  (8.12)  become  periodic  solutions  of  (8.8)  and  periodic  orbits 
of  (8.12)  become  invariant  tori  for  (8.8).  Around  the  homoclinic  orbit,  one  expects  a  be¬ 
havior  similar  to  the  one  discussed  in  §7.  The  precise  behavior,  of  course,  depends  upon 
the  higher  order  terms  in  (8.8)  and  the  method  of  analysis  will  be  similar  to  the  one  in  §7. 

Equations  (8.8)  are  the  generic  situation  for  two  purely  imaginary  and  one  zero  eigen¬ 
value.  If  further  symmetries  occur  in  the  problem,  there  may  be  no  second  order  terms  in 
the  normal  form  for  the  vector  field.  In  this  case,  the  simplest  case  for  the  approximate 
equations  are 

p  =  p(  X  -  ar2  -  by2 ),  y  =  y(0  +  cr 2  +  <ly2 ) 

with  a,  b,  c,  d  fixed  nonzeto  constants  and  X,  0  small  bifurcation  parameters.  In  the  case 
of  a  fourth  order  equation  with  two  purely  imaginary  roots,  the  same  equations  occur  in  a 
natural  way  coupled  with  two  angle  variables.  The  complete  bifurcation  diagram  for  this 
equation  cannot  be  obtained  without  the  addition  of  terms  of  order  five.  It  is  much  more 
difficult  and  the  reader  is  referred  to  Holmes  (2) ,  Guckenheimer  [4] .  Chow  and  Hale  1 1 1 
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9.  A  framework  for  infinite  dimensions 

For  infinite  dimensional  systems,  the  generic  theory  analogous  to  the  one  discussed  in 
§3  is  in  its  infancy,  in  this  section,  we  outline  an  approach  to  the  development  of  such  a 
theory  for  a  special  class  of  semigroups  of  transformations.  This  class  is  general  enough  to 
include  some  types  of  parabolic  and  hyperbolic  partial  differential  equations  as  well  as  re¬ 
tarded  functional  differential  equations  and  some  neutral  functional  differential  equations. 

Let  X,  Y,  Z  be  Banach  spaces  and  let  Xr  =  Cr(Y,  Z),r>  1 ,  be  the  set  of  functions 
from  Y  to  Z  which  are  bounded  and  uniformly  continuous  together  with  their  derivatives  up 
through  order  r.  We  impose  the  usual  topology  on  Xr.  For  each  / 6  Xr,  let  T^t):  X  — <■  X, 
t  >  0  be  a  strongly  continuous  semigroup  of  transformations  on  X.  For  each  x  £  X,  we 
suppose  TAt)x  is  defined  for  each  t  >0  and  is  Cr  in  x. 

In  applications,  one  often  is  interested  in  open  subsets  of  X,  Y,  but  this  presents  mainly 
notational  difficulties  and  therefore  will  not  be  discussed. 

We  say  that  a  point  x0  has  a  backward  extension  (relative  to  7y(f))  if  there  is  a  func¬ 
tion  0:  (-<*>,  0]  — *  X  such  that  0(0)  =  x0  and  Tf(t)<p(r )  =  0(f  +  r)  for  all  0  <  t  <  -  t. 
t  G  (-<*>,  0) .  If  x0  has  a  backward  extension  0,  we  define  T^(t)x0,  t  <  0,  by  the  relation 
Tf(t)x0  =  0(r)  and  say  that  Tf(t)x0  is  defined  for  t  <  0. 

The  positive  orbit  0+(xQ)  through  x0  is  O*  (x0)  =  U,  >0  Tf(t)x0 .  The  negative  orbit 
O  ~(jc0)  through  x0  is  O  ~(x0)  =  U,<07y(0*0  if  Tf(t)x0  is  defined  for  /  0.  The  u> -limit 

set  w(x0)  of  x0  and  a-limit  set  a(x0)  of  x0  are  defined  as 

w(*o>=  D  cl  U  Tf(t)x0,  a(x0)=  D  Cl  U  Tf(t)x0. 

T>0  f>r  T<0  t 

A  set  M  C  X  is  invariant  ( relative  to  T^t))  if,  for  any  x0  €  M,  T^(t)x0  is  defined  for 
t  <0  and  Tf(t)x0  €  M  for  all  t  6  (-»,  °°).  An  equilibrium  point  of  Tf(t)  is  an  invariant 
set  consisting  of  a  single  point;  that  is,  an  x0  G  X  such  that  Tf(t)x0  =  x0  for  all  /  6  R.  A 
periodic  orbit  y  of  Tf(t)  is  an  invariant  set  which  is  a  closed  curve;  that  is.  there  are  an 
x0  G  X  and  an  w  >  0  such  that  Tf(t)x0  =  Tf(t  +  w)x0  for  t  G  R,  to  is  the  least  period, 
and  y  =  U ,Tf{t)x0. 

The  following  result  is  easy  to  prove  (see,  for  example.  Hale  [5] ). 

Proposition  9.1 .  If  0+(x0)  is  precompact,  then  «(x0)  is  a  nonempty  compact  con¬ 
nected  invariant  set  and  (list (Tf(t)x0,  w(x0))  — *  0  as  t  — 
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(9.1)  A{  =  {jc  £  X:  Tj\t)x  is  defined  and  bounded  for  r  <  0} . 

This  set  contains  all  equilibrium  points,  periodic  orbits  and  the  co-limit  sets  of  any  compact 
orbit.  The  set  A ,  is  invariant.  For  any  bounded  open  set  N  C  X.  let 

(9.2)  ■4|-(/V)  =  {.r  6  .V  Tf(t)x  is  defined  and  belongs  to  N  for  /  <  0} 

Pkoi*osition  9  .2 .  If  A  ^  is  compact,  it  is  maximal  compact  invariant  and  attracts  points 
of  X.  If,  in  addition.  Tf{t)  is  one-to-one  on  A then  T^(t)  is  a  continuous  group  on  Af 

The  proof  is  elementary. 

Di-tinition  9.3.  For  /,  g  £  Xr.  we  say  f  is  equivalent  to  g,  f  ~  g,  if  there  is  a  home- 
omorphism  h :  A f  — »  Ag  such  that  h  maps  the  orbits  of  Tf(t)  on  A f  onto  the  orbits  of  Tg{t) 
on  Ag  preserving  the  sense  of  direction  in  time.  An  /€  Kr  is  structurally  stable  if  there  is 
a  neighborhood  V  of  /  such  that  f  ~  g  for  every  g  £  V.  An/£  Xr  is  a  bifurcation  point 
if  it  is  not  structurally  stable. 

Du  iNiTtoN  9.4.  For  f  g  £  Xr.  we  say  /  is  locally  equivalent  to  £  in  a  neighborhood 
of  ,r0  if  there  are  neighborhoods  N.  Af  of  ,t0  and  a  homeomorphism  h  Af(N)  — *  Ag(M) 
such  that  h  takes  the  orbits  of  Tf(t)  on  Af(N)  onto  the  orbits  of  Tg(t)  on  Ag(M)  preserving 
the  sense  of  direction  in  time.  An  /£  Xr  is  locally  structurally  stable  if  there  is  a  neighbor¬ 
hood  V  of  /  such  that  /  is  locally  equivalent  to  g  at  x0  for  every  g  £  V 

It  is  important  to  notice  that  the  flows  defined  by  Tf(t).  Tg(t)  are  compared  only  on 
the  invariant  sets  Af.  Ag.  However,  in  order  for  there  to  be  a  homeomorphism  between 
and  Ag  for  each  g  in  a  neighborhood  of  f  the  set  Af  must  have  some  strong  type  of  stability 
as  an  invariant  set  of  Tf(t)  We  shall  see  several  illustrations  of  this  remark  in  infinite  di¬ 
mensional  problems,  but  the  idea  is  easily  understood  by  comparing  this  definition  with  the 
one  given  in  §3  for  ordinary  differential  equations  in  Rn. 

Consider  the  equation 

(9.3)  X  =  Bx  +  fix) 

where  .t  £  R",  /( 0)  =  0,  df(0)/3x  =  0.  If  W“  is  the  unstable  manifold  of  zero  (it  could 
contain  only  llie  point  zero),  suppose  there  is  a  neighborhood  N  of  x  ~  0  such  that  if 
x0  £  N\W“.  then  there  is  a  r  <0  such  that  Tf(t)x0  6  d/V.  It  follows  that  Af(N)  ~  W“  n  /V. 

If /is  locally  structurally  stable  at  zero  according  to  Definition  9.4,  then  one  can  show  that 
the  equilibrium  point  x  =  0  of  (9.3)  must  he  hyperbolic;  that  is.  Re  \B  =£  0.  Also,  there  is 
a  neighborhood  V  of  / such  that,  for  each  g  £  V,  there  are  a  neighborhood  Af  of  zero,  a 
unique  equilibrium  point  xg  of  g  £  Af  and  /^(Af)  =  W14^)  O  Af.  In  Definition  9.4,  the 
trajectories  in  a  full  neighborhood  of  zero  are  not  compared,  and  they  are  not  even  com¬ 
pared  on  the  stable  manifold.  It  frequently  happens  in  applications  in  infinite  dimensions 
that  the  stable  manifold  is  infinite  dimensional  and  the  unstable  manifold  is  finite  dimen¬ 
sional.  Comparison  of  trajectories  on  the  infinite  dimensional  parts  seems  to  be  impossible 
Thus,  the  definition  of  equivalence  is  chosen  as  above. 
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To  define  a  large  class  ol  semigroups  Tf(t)  which  will  guarantee  that  A}  is  compact, 
maximal  invariant  and  satisfies  some  stability  properties,  we  need  some  additional  notation. 

Let  q (B)  be  the  Kuratowskii  measure  of  noncompactness  of  a  bounded  set  B  in  a 
Banach  space  X\  that  is,  a(B )  =  inf{J  >  0:  B  has  a  finite  cover  with  each  element  of  diam¬ 
eter  <  J}.  II  T'.  X  ♦  X  is  continuous,  we  say  7"  is  a  conditional  a-contraction  if  there  is  a 
constant  A:  €  (0,  I)  such  that  a(7B)  <  ka(B)  for  ail  bounded  sets  B  in  X  for  which  TB  is 
bounded.  The  map  is  an  a-contraction  if  it  is  a  conditional  a-contraction  and  takes  bounded 
sets  into  bounded  sets.  For  general  properties  of  the  measure  of  noncompactness  and  a- 
contractions,  see  Sadovskii  [lj  or  Martin  Jlj .  A  map  T  is  conditionally  completely  con¬ 
tinuous.  if,  for  any  bounded  set  B  for  which  TB  is  bounded,  the  closure  of  TB  is  compact. 

It  7V),  /  >  0,  is  a  family  of  maps,  we  say  it  is  a  conditional  a-contraction  |or  conditionally 
completely  continuous)  if  for  each  t  and  each  bounded  set  B  for  which  [T(t)B,  0  <  r  t\ 
is  bounded.  a(7’(r)S)  <  k(t)a(B)  for  some  k(t)  6  )0,  1)  [or  the  closure  of  T(t)B  is  com¬ 
pact)  . 

If  re(T)  denotes  the  radius  of  the  essential  spectrum  of  a  linear  operator  T  and 
re(T)  <  1 ,  then  there  is  an  equivalent  norm  in  X  such  that  T  -  S  +  U  where  S.  U  are  linear 
operators  with  |S|  <  1  and  (/compact  (see  Leggett  [1] ,  [2) ,  Massatt  [3)).  Thus,  there  is 
an  equivalent  norm  in  X  such  that  T  is  an  a-contraction. 

The  semigroup  Tf(t)  is  a  conditional  a-contraction  if  it  is  a  conditional  a-contraction 
for  each  t  >  0.  If  Tf(t)  is  linear  with  re(Tf(t))  <  exp(-(3f)  for  some  p  >  0.  then  there  is  an 
equivalent  norm  in  X  such  that  Tf(t)  is  an  a-contraction  for  each  i  >  0. 

Proi'osition  9.5.  If  Tj-(t)  is  a  conditional  a-contraction,  then  each  bounded  orbit  is 
precompact. 

The  proof  is  trivial  since  O+(x0)  is  invariant  under  T^(t)  for  each  /  >  0. 

A  set  K  C  X  is  said  to  attract  a  set  M  C  X  (relative  to  Tf(t))  if  dist(7>(/)A/,  A  )  0 

as  /  — »  oo  A  set  A  in  X  is  stable  (relative  to  T^(t))  if,  for  any  neighborhood  V  of  A,  there 
is  a  neighborhood  V  of  A  such  that  Tf(t)V  C  U,  t  >  0.  A  set  A  in  X  is  asymptotically 
stable  ( uniformly  asymptotically  stable)  if  it  is  stable  and  there  is  a  neighborhood  W  of  A 
such  that  A  attracts  points  of  W  (A  attracts  IV). 

The  semigroup  T^(t)  is  point  ( local)  ( compact )  dissipative  if  there  is  a  bounded  set  II 
in  X  such  that  B  attracts  each  point  (some  neighborhood  of  each  point)  (each  compact  set) 
of  X. 

Proposition  9.6.  If  Tf(t)  is  a  conditional  a-contraction  and  compact  dissipative, 
then  Af  is  a  maxima!  compact  invariant  set  and  Aj  is  uniformly  asymptotically  stable.  If 
Tf(t)  is  a  conditional  a-contraction,  point  dissipative  and  there  is  some  neighborhood  Ox  of 
each  x  £  X  such  that  U,  >0  Tf.(t)Ox  is  bounded,  then  the  same  conclusion  holds.  If  Tr(t ) 
is  point  dissipative  and  conditionally  completely  continuous,  then  the  same  conclusion  holds. 
If,  in  addition,  the  orbits  of  bounded  sets  are  bounded,  then  A ^  attracts  bounded  sets  of  A'. 

For  a  proof  of  these  results,  see  Cooperman  [I),  Hale  [6|,  Massatt  |1J. 
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In  many  applications,  the  semigroup  Tf(t)  is  not  a  general  a -contraction  but  satisfies 
the  condition  T^{i)  =  S^(t)  +  U^(t)  where  Sy(r)  is  a  linear  contraction  for  each  /  >  0  and 
the  mapping  U^(t)  is  conditionally  completely  continuous.  Also,  the  family  of  mappings 
Tf(t)  is  a  semigroup  on  two  Banach  spaces  Xx ,  X2  with  Xx  compactly  imbedded  in  X2 . 
Suppose  also  that,  for  any  set  B  in  Xx  for  which  B  and  Uf(l  )B  are  bounded  in  X2,  it  fol¬ 
lows  that  Cfj-(l)  is  bounded  in  Xx.  In  this  situation,  Massatt  [5)  has  recently  proved  the  im¬ 
portant  result  that  point  dissipative  in  X2  is  equivalent  to  compact  dissipative  in  X2 .  This 
is  the  same  result  as  the  one  that  is  stated  in  Proposition  9.6  for  the  case  when  SM)  =  0  for 
all  r.  Other  properties  of  the  same  type  of  maps  have  been  discussed  by  Massatt  |2) . 

Our  primary  objective  is  to  study  how  the  set  A ^  varies  with  f.  The  easiest  result  to 
obtain  is  the  upper  semicontinuity  in  /  and  is  stated  precisely  in  the  following 

Proposition  9.7.  Suppose  there  is  a  neighborhood  V  of  f  and  a  bounded  set  B  C  X 
such  that,  for  each  g  6  V,  Tg(t)  is  an  a-contraction,  B  attracts  compact  sets  of  X  relative  to 
Tg(t)  and  { Tf(t)H,  t  >  0}  is  bounded  for  each  bounded  set  H  C  X.  Then  Af  is  upper  semi- 
continuous  at  f\  that  is,  for  any  neighborhood  U  of  A^,  there  is  a  neighborhood  W  of  f  such 
that  Ag  £.  U  for  g  €  W. 

The  idea  of  the  proof  is  the  following.  From  Proposition  9.6,  Ag  is  compact  and 
Ag  Q  B  for  every  g  6  V.  Also,  Proposition  9.6  implies  Af  is  uniformly  asymptotically  stable 
and  attracts  B.  This  is  enough  to  complete  the  proof.  Results  similar  to  Proposition  9.7  for 
locally  compact  spaces  have  been  proved  by  Marchetti  et  al.  f  1 ) . 

Before  proceeding  further,  let  us  give  some  examples  of  semigroups  which  are  o -con¬ 
tractions.  Suppose  A  is  a  linear  operator  (bounded  or  unbounded)  for  which  -  A  is  the 
infinitesimal  generator  of  a  strongly  continuous  semigroup  S(t)  on  Z.  Suppose  A"  is  a  Banach 
space  which  can  be  continuously  imbedded  in  Z  and  f.X  — ►  Z  is  a  given  function.  The  ex¬ 
amples  will  be  special  cases  of  an  evolutionary  equation 

(9.4)  u  +  Au  =  jXu ) 

which  generates  a  semigroup  7y(/)  for  which  the  variation  of  constants  formula  holds, 

(9.5)  Tf(t)x  =  S(f)x  +  JqSU  -  r)f(u(r))dT  d=f  5(/)jr  +  U(t)x 

where  the  radius  rf{S(t))  of  the  essential  spectrum  of  S(t)  satisfies 

(q  h)  re(S(t))  <  exp(- pt),  t>  0, 

for  some  0  >  0  and 

(9.7)  U{t)  is  conditionally  completely  continuous. 

From  (9.6),  (9.7),  it  follows  that  there  is  an  equivalent  norm  with  the  property  that 
Tf[t)  is  a  conditional  o -contraction. 

The  idea  for  the  proof  is  simple.  Since  rr{S(t))  <  cxp(-0 /)  for  any  r  >  0.  the  space 
Z  can  be  decomposed  as  Z  =  Z,  ©  Zz  where  Z,.  Zz  arc  independent  of  t,  invariant  under 
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S(t),  Z2  is  finite  dimensional  and  the  spectrum  of  S(t)  restricted  to  Z2  is  inside  the  disk 
centered  at  zero  with  radius  exp(-^r)  with  p  >  0.  Now  one  can  renorm  so  that  £(/)  re¬ 
stricted  to  Z2  has  norm  <  exp(-0f),  (3  >  0. 

Our  first  example  concerns  the  case  where  A  is  sectorial.  A  linear  operator  A  on  Z  is 
sectorial  if  it  is  closed,  densely  defined,  such  that .  for  some  <j>  in  (0,  rr/2)  and  M  >  l  .a  C  R. 
the  sector 

Sa.<t>  =  {X:  0  <  I arg(  A  -  a)\  <  rr,  X  *  a } 

is  in  the  resolvent  set  of  A  and  |(X  -  A)-1 1  <  M/  |X  -  a|  for  X  €  A  semigroup  of 
operators  S(t)  on  Z  is  analytic  if  the  map  t  t-»  T(t)z  is  real  analytic  on  0  <  /  <  °°  for  each 
z  e  Z.  If  A  is  sectorial,  then  -A  is  the  generator  of  an  analytic  semigroup  S(t)  and  con¬ 
versely  (see  Friedman  [I],  Henry  [  1  ] ,  Pazy  [1],  Martin  |1]).  In  addition,  if  A  lias  compact 
resolvent,  then  S(r)  is  compact  for  t  >  0,  re(S(t))  =  0  for  t  >  0  and  so  relation  (9.6)  in  an 
appropriate  norm  is  satisfied  for  any  p  >  0. 

If  A  is  sectorial,  then  one  can  define  the  fractional  powers  (A  +  al)a  of  A  +  al  and 
the  spaces  Za  =  D((A  +  alf  )  with  the  graph  norm.  |i/eX'  =  Cr(Za,  Z),  0  <a  <  1 . 
r  >  1,  then  equation  (9.4)  generates  a  semigroup  TXt)  on  Z°  satisfying  (9.5),  (9.7)  provided 
that  /  takes  bounded  sets  into  bounded  sets  and  the  resolvent  set  of  A  is  compact.  The  spaces 
X,  Y,  Z  are  X  =  Y  =  Z“,  Z  =  Z.  To  satisfy  (9.6),  additional  restrictions  must  be  imposed  on  A. 
A  specific  example  of  the  latter  situation  is 

(9.8)  du/dt  -  Au  =  /(*,  u,  ux)  in  12,  u  =  0  on  dfl 

where  S2  is  a  bounded  open  set  in  R"  with  smooth  boundary, /e  Xr  =  Cr(Sl  x  R  *  R.  R). 
In  W  =  L2(£l,  R),  the  operator  A  =  -  A  with  domain  H'0(£l)  n  H2(Sl)  =  n  h'z  :(i2) 

is  sectorial  with  compact  resolvent.  For  some  restrictions  on  the  rate  of  growth  of  /(•,«.  u) 
in  u,  v,  there  is  an  a  in  (0,  1)  such  that  equation  (9.8)  generates  a  strongly  continuous  semi¬ 
group  on  Wa .  In  this  case  X  =  If,  f  =  fi  *  R  x  R,  Z=R. 

Other  problems  that  could  be  considered  for  (9.5)  are  cases  when  fix,  u.  u)  is  indepen¬ 
dent  of  u;  that  is,  Y  =  fl  x  R,  or  f(x,  u,  v)  independent  of  x,  u;  that  is,  Y  =  R.  As  the 
function  /  is  restricted  to  a  smaller  class,  the  generic  theory  will  become  more  difficult. 

One  could  also  change  the  boundary  conditions  to  some  other  form  or  consider  sys¬ 
tems  of  parabolic  equations  and  obtain  semigroups  of  the  same  qualitative  type. 

With  O  as  in  (9.8),/ £  Xr  =  C7r(J2  x  R,  R),  p  >  0  a  constant,  consider  the  equation 

(9.9)  utt  -  pAut  -  Au  =  f(x,  u )  in  T2,  u  =  0  on  3J2. 

Webb  (I )  has  shown  that  (9.9)  generates  a  semigroup  Tf(t)  on  D(A)  x  l2{il,  R),  A  =  A, 
D(A)  =  //q (£1)  n  H2(Yl).  It  is  possible  to  show  that  T^(t)  also  satisfies  (9.5)— (9.7)  (see, 
for  example,  Massatt  [4]). 

In  an  appropriate  Sobolev  space,  there  is  a  semigroup  generated  by  the  beam  equation 


(9.10) 
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satisfying  (95),  (9.6).  Relation  (9.7)  does  not  appear  to  be  satisfied.  Ball  [1) ,  however, 
has  obtained  interesting  qualitative  results  on  this  equation  exploiting  (9.5),  (9.6)  and  con¬ 
vergence  in  the  weak  topology.  This  suggests  that  axiom  (9.7)  should  be  weakened  in  some 
way. 

Retarded  functional  differential  equations  generate  semigroups  satisfying  (9.5)— (9.7). 
Suppose  r  >  0,  C  =  C\\-r,  Oj ,  Rn),/e  Xk  =  Ck(C,  R"),  k  >  1  and  consider  the  equation 

(9.11)  x(t)  =  f(xt) 

where  xt(0)  ~  x(t  +  8),  -r  <  6  <  0.  For  any  0  £  C,  there  is  a  unique  solution  x(8)  of 

(9.1 1)  with  initial  value  0  at  r  =  0.  If  we  assume  that  x(0)(r)  is  defined  for  t  >  0  and  let 
7/(r)0  =  xf(0),  t>  0,  then  7y(/)0  is  a  strongly  continuous  semigroup  on  C  Let  S{t)  be 
the  semigroup  on  C  generated  by  x(t)  =  0;  that  is,  S(t)  =  T0(r), 

\  0(f  +  9),  t  +  8  <  0, 

S(f)<*0)=  ' 

( 0(0),  r  +  e  >  0. 

Then  re(S(r))  =  0  for  /  >  0.  Also, 

Tfu)m  =  mm  +  j‘0s(t  -  ^x^ewT^mdr 

where  X0(6)  -  0  for  6  <  0,  =  l  for  6  =  0,  S(t  -  r)X0(6)  =  0  for  t  -  r  +  0  <  0,  =  /  for 
t  -  t  +  8  >  0.  Thus, 

(9.12)  7>(r)0  =  5(r)0  +  JQfS(r  -  T)X0f(Tf(r)<p)dT  +  f/(r)0 

with  the  interpretation  as  above.  Relation  (9.12)  is  the  analogue  of  (9.5)  and  one  can  show 
that  £/(/)  is  conditionally  completely  continuous  if  /  takes  bounded  sets  into  bounded  sets 
(see  Hale  (3]).  In  this  case,  our  spaces  X,  Y,  Z  are  X  -  C  =  Y,  Z  =  R". 

For  the  difference-differential  equation 

(9.13)  x(f)  =  /(x(r),  x(r  -  r)) 

the  spaces  X,  Y.  Z  are  X  =  C,  Y  —  R”  *  R",  2  =  R"  For  the  equation 
(9  .14)  x(r)  = /(x(f  -  r)) 

the  spaces  are  X  =  C,  Y  =  R",  Z  =  R" . 

It  is  possible  to  discuss  these  retarded  equations  in  other  spaces  than  C,  for  example. 
R"  x  l-2([-r.  01 ,  R").  Relation  (9.12)  and  rr{S(t))  =  0  for  t  >  0  still  remain  true  (see 
Hale  [3]  for  references). 

For  the  case  of  infinite  delay  in  equation  (9.1 1),  it  is  not  difficult  to  obtain  spaces  X 
of  initial  data  for  which  one  obtains  a  strongly  continuous  semigroup  7^(0  on  X  satisfying 
(9.12)  with  S(t )  a  strongly  continuous  semigroup  and  U(t)  conditionally  completely  contin¬ 
uous.  For  the  number  r,(5(/))  to  satisfy  (9.6),  one  must  impose  some  additional  conditions 
on  the  space  X.  For  example,  if  X  is  a  fading  memory  space,  the  kernel  approaching  zero 
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exponentially  is  sufficient.  Many  other  spaces  assure  that  (9.6)  is  satisfied  (see  Hale  and 
Kato  [  1  ] ,  Schumacher  [  1  ] ).  The  survey  articles  of  Hale  [4] ,  Corduneanu  and  Laksmikan- 
tham  [1]  should  be  consulted  for  references  and  more  specific  properties  of  equations  with 
infinite  delays. 

Some  functional  differential  equations  of  neutral  type  also  will  generate  semigroups 
satisfying  all  of  the  properties  mentioned  above.  For  example,  this  will  be  true  for  the 
equation 


(9.15) 


Jt  •*(')-  £  Akx(t  -  rk)- h(f°rA(6)x(t  +  6)dd'j 


=  Six, ) 


with  h  e  C*( R,  R"),  g 6  Ck(C,  R"),  rk,  r  >  0,  A(&)  a  continuous  n  x  n  matrix,  each  Ak  is 
an  n  x  n  constant  matrix  and  the  zero  solution  of  the  difference  equation 


(9.16) 


AO  -  £  /ifcy(t-rt)  =  0 


k  =  l 


is  uniformly  asymptotically  stable.  The  semigroup  S(t)  is  the  one  generated  by  the  equation 


1 1*')  -  £i  AkAt  ~  r*)J  =  0. 

The  parameter  / in  the  semigroup  for  (9.15)  is  {Ak,  A(),  g,  h} .  These  equations  arise 
naturally  from  certain  linear  hyperbolic  partial  differential  equations  with  nonlinear  bound¬ 
ary  conditions  (see  Hale  [3]  for  references  and  details). 

For  g  =  0  in  (9.15),  one  also  obtains  semigroups  generated  by  the  functional  equation 


(9-17)  y(t)  -  £  Akx(t  -  rk)  -  h  °  A(0)x(r  +  0)c/<?)  =  0 

with  initial  data  <j>  restricted  to  C0,  where  Ca  for  a  e  R"  is  defined  by 

Ca  =  J0  €  C:  <K0)  -  ZAkA-rk)  -  h(f°_rA(0)mdo)  =  a 

In  this  case,  C0  is  a  closed  subset  of  C-  a  nonlinear  manifold  which  depends  on  the  coeffi¬ 
cients  Ak  and  the  functions  ,4,  h.  From  the  qualitative  point  of  view,  it  may  be  desirable 
not  to  consider  (9.17),  but  to  consider  (9.14)  with  g  =  0  realizing  that  Ca  for  each  a  e  R” 
is  invariant  under  the  semigroup. 

There  is  a  tremendous  literature  on  semigroups  generated  by  evolutionary  equations 
with  delays-a  combination  of  equations  (9.4)  and  (9.11).  Both  the  theory  and  applications 
are  fairly  well  developed  and  it  is  impossible  here  to  go  into  this  interesting  subject.  A  repre¬ 
sentative  selection  of  the  literature  can  be  obtained  by  consulting  the  references  in  Slemrod 
(l),  (2],  Infante  and  Walker  (1),  Dafermos  and  Nohel  (1),  Fitzgibbon  [  1 } . 
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The  above  examples  certainly  are  sufficient  motivation  to  study  in  more  detail  semi¬ 
groups  Tf(t)  which  satisfy  (9.5)-(9.7).  All  questions,  remarks,  conjectures,  etc.  in  the  fol¬ 
lowing  pages  are  made  for  semigroups  which  satisfy  at  least  these  conditions.  Proposition  9.2 
motivates  the  following  query. 

Question  9.8.  If  A  f  is  compact,  when  is  there  a  generic  set  of  f  such  that  Tf(t)  is 
one-to-one  on  Ay? 

Note  that  the  question  is  for  Ay  and  not  all  of  X.  It  does  not  seem  to  be  reasonable 
to  ask  the  same  question  for  all  of  X. 

For  parabolic  partial  differential  equations,  general  conditions  are  known  which  imply 
that  Ty(f)  is  one-to-one  on  all  of  X  (see  Henry  [1] ,  K.  Miller  [1] ;  for  further  references,  see 
Manselli  and  Miller  [1]).  In  particular,  it  is  true  if  the  elliptic  part  of  the  operator  is  the 
Laplacian.  For  retarded  functional  differential  equations  (9.1 1)  with  /  analytic,  T^(t)  is  one- 
to-one  on  Af  even  though  it  may  not  be  one-to-one  on  C.  This  follows  because  the  func¬ 
tions  defining  must  be  analytic  (see  Nussbaum  [1  ] ).  The  same  result  is  true  for  neutral 
equations  with  (9.16)  uniformly  asymptotically  stable  (see  Hale  [3]).  For  other  references 
on  retarded  equations,  see  Hale  [3] .  Mallet-Paret  has  given  an  example  (unpublished)  of  a 
retarded  functional  differential  equation  for  which  Tj-(t)  is  not  one-to-one  on  Ay. 

Question  9.9.  If  f  is  structurally  stable,  when  is  Tf(t)  one-to-one  on  Ay? 

Question  9. 1 0.  When  is  Af  generically  a  manifold  or  a  finite  union  of  manifolds' ? 

For  the  case  of  retarded  equations  (9.1 1),  defined  on  a  compact  manifold  M  without 
boundary,  which  are  in  some  sense  close  to  an  ordinary  differential  equation,  Kurzwcil  [1] 
has  shown  that  Ay  is  diffeomorphic  to  Af.  Oliva  [l  j  has  generalized  these  results  giving 
other  conditions  which  imply  Ay  is  diffeomorphic  to  Af.  Henry  ( 1  ]  has  discussed  this  ques¬ 
tion  for  certain  gradient  systems  of  parabolic  equations  and  shown  that  Af  is  the  union  of  a 
finite  number  of  manifolds.  We  mention  later  special  retarded  functional  differential  equa¬ 
tions  for  which  Ay  is  the  union  of  a  finite  number  of  manifolds. 

Of  course,  an  affirmative  answer  to  cither  of  the  above  questions  requires  conditions 
on  the  manner  in  which  Tj(t)  depends  on  f  and  a  specification  of  a  class  Xr  of  f  s  which 
gives  enough  flexibility  to  move  the  orbits  in  any  desired  direction  by  variations  in  f. 

Theorem  9.1 1 .  //Ay  is  compact  and,  for  each  t  >  0,  x  €  Ay.  Dx  T^(t)x  is  the  sum 
of  a  contraction  operator  and  a  completely  continuous  operator,  then  Ay  has  finite  Haus- 
dorff  dimension. 

For  X  a  separable  Hilbert  space,  this  result  was  proved  by  Mallet-Paret  1 1 J .  The  gen¬ 
eral  case  in  Theorem  9.1 1  for  X  an  arbitrary  Banach  space  was  proved  by  MafTe  [2] .  Lady- 
zenskaya  (!)  has  obtained  related  results  for  the  Navier-Stokes  equation.  Mand  [2]  proved 
even  more-that  the  limit  capacity  of  Af  is  finite.  The  limit  capacity  of  Ay  is  defined  as 
the  exponential  growth  rate  of  the  number  of  balls  of  radius  exp(-/)  that  are  required  to 
cover  A y  when  /  — *  °°.  The  limit  capacity  is  always  at  least  as  large  as  the  Hausdorff  di¬ 
mension.  Using  results  of  Cartwright  [2],  [3],  one  obtains  the  following  interesting  conse¬ 
quence  of  Theorem  9.1 1 . 
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Corollary  9.12.  If  the  conditions  of  Theorem  9.1 1  are  satisfied,  then  there  is  an 
integer  N  such  that,  if  there  is  an  x  €  X  such  that  Tf(t)x  is  almost  periodic,  then  Tf{t)x  is 
quasiperiodic  with  the  number  of  basic  frequencies  being  <  N. 

Question  9.13.  For  x  €  A^,  when  is  T^(t)x  continuously  differentiable  in  t  for 

r  e  (-«,  °°)? 

Note  again  that  this  question  pertains  to  A ^  and  not  all  of  X. 

For  retarded  equations  (9.1 1),  this  is  obviously  true  because  Tf(t)x  is  defined  for 
t  <  0.  For  neutral  equations  for  which  equation  (9.16)  is  uniformly  asymptotically  stable 
and  A(6)  has  a  continuous  derivative,  it  is  also  true.  The  proof  uses  the  fact  that  the  semi¬ 
group  S(t)  in  (9.S)  has  essential  spectrum  inside  the  unit  circle  for  t  >  0  (see  Hale  [3]).  It 
can  be  generalized  to  a  neutral  equation  with  the  essential  spectrum  of  S(f)  bounded  away 
from  the  unit  circle  for  t  >  0  (see  deOliveira  [1]).  Some  special  cases  also  have  been  con¬ 
sidered  for  functional  equations  by  Hale  and  deOliveira  ( 1  ] .  We  remark  that  TM) x  is  dif¬ 
ferentiable  in  t  for  x  £  Af,  but  is  not  differentiable  for  every  x  €  C. 

For  abstract  evolutionary  equations,  no  attempt  has  been  made  to  show  that  7y(t)x 
is  differentiable  in  t  for  x  only  in  Aj.  One  approach  has  been  to  show  that'7y(/)x  is  con¬ 
tinuously  differentiable  jointly  in  t,  x  for  t  >  0  and  x  E  Y,  a  Banach  space  continuously  im¬ 
bedded  in  X.  Typical  hypotheses  require  that  7y(f)  is  a  semigroup  also  on  a  Banach  space  V 
which  can  be  continuously  imbedded  in  X,  Y  belongs  to  the  domain  of  the  generator  A  of 
Tf(t)  and  is  a  bounded  operator  from  Y  to  X  (see  Dorroh  and  Marsden  [1] ,  Marsden  and 
McCracken  [1]).  For  the  case  where  A  is  sectorial,  Henry  [1]  gives  more  detailed  results  on 
the  differentiability  in  f. 

It  should  be  easier  to  prove  differentiability  in  t  on  Aj  and  there  should  be  an  abstract 
generalization  of  the  known  results  for  functional  differential  equations.  It  will  be  impossi¬ 
ble  to  obtain  any  type  of  generic  theory  without  having  this  type  of  smoothness  since  one 
cannot  study  the  local  behavior  of  orbits  on  Aj  by  taking  a  linear  approximation. 

Our  next  objective  is  to  discuss  more  detailed  properties  of  Aj  and,  in  particular,  hy- 
perbolicity. 

Definition  9.14.  An  equilibrium  point  x  is  hyperbolic  if  there  are  a  decomposition 
of  X  as  X  =  Xs  ©  Xu  and  positive  constants  K,  a  such  that  Xs  n  Xu  =  { 0},  Xs,  Xu  are 
invariant  under  DTf(t)x  and 

\(DTf(t)x)v\  <  Ke -**,  t  >  0,  u  6  X1 ,  \(DTf(t)x)v \  «  t  <  0,  u  6  Xu . 

It  is  implicitly  assumed  that  DTf(t)x  is  defined  on  Xu  for  t  <  0.  If  Tf(t)  is  an  a- 
contraction,  then  DT^(t)x  is  an  a -contraction,  the  space  Xu  is  finite  dimensional  and 
DTf(t)x  is  defined  on  Xu  for  t  <  0. 

Suppose  T  =  { Tf{t)x,  0  <  t  <  w}  is  an  w-periodic  orbit  for  which  T^{t)x  is  continu¬ 
ously  differentiable  in  t.  (This  is  one  reason  we  need  an  answer  to  Question  9.13.)  Let  // 
be  a  hyperplane  in  X\  that  is,  H  =  [y  €  X:  h(y)  =  0  where  h  is  a  nontrivial  continuous 
linear  functional  on  Jf).  Suppose  there  are  a  neighborhood  W  of  zero  in  H  and  a  function 
0:  W  — *  X  such  that  0(0)  =  x  and  <tf  W)  is  a  C1 -submanifold  of  codimension  one.  let  M 
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be  the  tangent  space  to  <t>(W)  at  x.  The  set  V  =  is  said  to  be  a  transversal  to  F  at  x 
if  dTf(t)x/dt  does  not  belong  to  M.  Let  rr  be  the  Poincare'  map  defined  on  a  neighborhood 
U  of  x  in  a  transversal  V  to  T  at  .x.  Then  ny  is  continuously  differentiable  for  y  £  U. 

Definition  9.15.  Suppose  the  notation  is  as  in  the  previous  paragraph.  A  periodic 
orbit  T  is  hyperbolic  if  there  are  a  decomposition  of  M  as  M  =  Ms  ©  Mu  and  positive  con¬ 
stants  K,  a  such  that  M3  O  Mu  =  {0},  M3,  Mu  are  invariant  under  Dir(x)  and 

|(Dn(x))"u|  <  Ke~an,  O0,u6W',  \(Dir(x))nvl  <  Kean ,  n  <  0,  v  e  Mu . 

One  can  define  the  stable  and  unstable  manifolds  for  hyperbolic  equilibrium  points  and 
periodic  orbits  in  the  usual  way.  For  functional  differential  equations  (see  Hale  [3])  and 
nonlinear  evolutionary  equations  (9.4)  with  A  sectorial  (see  Henry  [1]).  the  unstable  and 
stable  manifold  intersect  transversally  at  the  point  or  orbit,  the  unstable  manifold  is  an  im¬ 
mersed  submanifold  of  finite  dimension  and  the  local  stable  manifold  is  diffeomorphic  locally 
to  the  stable  manifold  of  the  linear  approximation.  In  general,  it  is  not  possible  to  define  a 
global  stable  manifold  because  the  solution  operator  7y(r)  is  not  one-to-one  for  negative  t 
when  it  is  defined.  An  abstract  generalization  of  these  results  is  needed. 

Definition  9.16.  Suppose  TAt)x  is  continuously  differentiable  in  t  if  x  belongs  to  a 
periodic  orbit.  Then  Tf(t)  is  said  to  be  Kupka-Sntale  if: 

(i)  All  equilibrium  points  are  hyperbolic. 

(ii)  All  periodic  orbits  are  hyperbolic. 

(iii)  The  local  stable  and  global  unstable  manifolds  of  equilibrium  points  and  periodic 
orbits  intersect  transversally. 

Question  9.17.  Is  the  set  of  f  6  Xr  corresponding  to  Kupka-Smale  semigroups  generic 
in  Xr 1 

Some  results  have  recently  been  obtained.  For  the  general  retarded  functional  differen¬ 
tial  equations  (9.1 1)  with  / e  Xr,  the  set  of  Cr-functions  from  C([-/i,  0) .  R")  to  R".  Mallct- 
Paret  [2] ,  (3 ]  has  shown  the  answer  to  the  question  is  affirmative.  II  /  E  Cr(R”  x  R".  R") 
as  in  (9.13),  he  has  also  shown  the  same  result  is  true.  If/€  C^R”,  R”)  as  in  (9.14),  the 
answer  to  the  question  is  not  known.  The  fundamental  part  of  the  proofs  always  involves 
(as  was  the  case  in  finite  dimension)  showing  that  the  set  of  /  for  which  the  periodic  orbits 
with  one  not  in  the  spectrum  of  Dtt(x)  is  dense. 

For  the  neutral  equations  discussed  above,  deOlivcira  [1]  has  shown  that  the/ such 
that  the  equilibrium  points  are  hyperbolic  is  open  and  dense.  The  general  KS  property  has 
not  been  discussed. 

For  the  parabolic  equation 

(9.18)  u,  =  uxx  +  fix,  u).  0  <  jp  <  1 . 

(9.19)  au(t,  0)  +  &ux(t.  0)  =  0,  yu(t.  1)  +  8ux(r,  1)  =  0, 

the  co-limit  set  of  every  hounded  orbit  is  a  single  equilibrium  point;  that  is,  a  solution  of 

(9.20)  nxx  +  fix,  u)  =  0 
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satisfying  the  boundary  conditions  at  x  =  0,  x  =  1  (see  Matano  [2]).  It  is  not  difficult  to 
show  that  these  equilibrium  points  are  hyperbolic  genetically  in  /£  C’fR  *  R,  R).  The 
same  result  has  recently  been  proved  by  Brunovsky  and  Chow  [1]  for  the  equation 

(9.21)  ut  =  uxx  +  /(u),  0  <  x  <  1 , 

with  Neumann  or  Dirichlet  boundary  conditions,  where  the  function  / €  C2(R,  R)  with  the 
Whitney  topology  is  independent  of  the  space  variable  x.  The  restriction  to  this  smaller  class 
of  functions  makes  the  result  nontrivial.  The  methods  appear  to  apply  to  the  general  bound¬ 
ary  conditions  (9.19). 

To  see  why  this  is  nontrivial,  consider  the  Neumann  boundary  condition  ux( 0)  = 
1^(1)  =  0  for  the  equilibrium  equation  (9.20).  If  ui(b,  /)  is  the  period  of  the  solution  of 
(9.20)  through  the  point  (u(0),  ux(0))  =  ( b ,  0),  then  one  must  show  that  the  solutions  b  of 
the  equation  u> (b,  /)  -  2=0  satisfy  3cj(3,  /)/3A  0  for  a  residual  set  of/e  C2( R,  R). 

There  are  classical  examples  of  nonlinear  functions  /  where  ui(b,  /)  is  a  constant  function  of 
b  (see,  for  example,  Urabe  [1]).  This  implies  that  the  derivatives  of  u>(b,  /)  with  icspect  to 
b  are  very  complicated  functions  of  the  derivatives  of /. 

For  the  equation 

(9.22)  ut  =  Au  +  f(u)  in  ft,  au  +  (idu/dn  =  0  on  3ft 

where  £2  is  a  bounded  open  connected  set  in  R"  with  smooth  boundary,  the  corresponding 
result  is  not  known.  In  this  case,  it  is  reasonable  to  consider  the  parameters  as  /  and  ft  and 
to  prove  that  the  equilibrium  points  are  hyperbolic  generically  in  (/  ft).  An  important  role 
should  be  played  by  results  similar  to  the  ones  of  Uhlenbeck  [I]  on  the  simplicity  of  the 
eigenvalues  of  the  Laplacian  operator  generically  in  ft. 

To  the  author’s  knowledge,  property  (iii)  in  the  definition  of  Kupka-Smalc  systems 
has  not  been  discussed  for  the  general  equation  (9.22).  Some  partial  results  of  Henry  |1] . 
[2)  are  mentioned  below. 

Definition  9.18.  Let  ft(/)  be  the  nonwandering  set  of  Tj-(r).  A  semigroup  TAt)  is 
said  to  be  Morse-Smale  if  ft(/)  is  the  union  of  a  finite  number  of  equilibrium  points  and 
periodic  orbits,  each  hyperbolic  with  stable  and  unstable  manifolds  intersecting  transvcrsally. 
Question  9.19.  Are  Morse-Smale  systems  open  and  structurally  stable ? 

One  of  the  first  important  problems  in  any  attempt  to  discuss  these  last  general  ques¬ 
tions  is  to  construct  several  examples.  The  first  step  would  be  to  consider  systems  which 
have  properties  analogous  to  gradient  systems  for  ordinary  differential  equations. 

Gradient  flows  defined  by  parabolic  equations  have  been  analyzed  by  Henry  (Ij.  (2|. 
In  particular,  it  is  shown  that  the  equation 

u,  =  uxx  +  Mu  -  u3),  0  <  x  <  ir,  t  >  0. 

(9.23) 

u  =  0  at  x  =  0,  n 

is  structurally  stable  in  //q(0,  1 )  for  X  €  (0,  9),  X  =£  1 ,  X  #  4,  is  a  point  for  0  <  \  <  I . 
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Ak  is  a  smooth  closed  arc  with  three  equilibrium  points  for  1  <  X  <  4,  A^  is  2-dimensional 
and  contains  5  equilibrium  points  for  4  <  X  <  9. 

It  is  not  difficult  to  see  why  this  should  be  true  if  we  use  our  detailed  knowledge  of 
the  maximal  compact  invariant  set.  In  fact,  for  any  X  >  0,  one  can  use  the  invariance  orm- 
ciple  and  the  function 


to  show  that  every  solution  is  bounded  and  approaches  a  solution  of 

(9.24)  uxx  +  \(u  -  u3)  =  0 

with  i/(0)  =  u(ir)  =  0.  For  each  X  =  0,  this  equation  with  these  boundary  conditions  has  a  finite 
number  of  solutions  {t/,,/  =  1 , 2, . . . ,  »(X)}  with  unstable  manifolds  Wu{Uj)  of  finite  dimension. 
One  can  show  that  this  implies  Ax  is  compact  and  attracts  bounded  sets,>4x  =  U/iV ^u(w/)> 
and  is  upper  semicontinuous  in  X.  If  X€(0,  1),  n(X)  =  l,u,  =  0,A^={0}.  ForX  =  4,the 
half-period  for  the  solutions  of  the  linear  part  of  (9.24)  is  rr.  For  X  =  4  and  close  to  4,  there  are 
two  additional  solutions  tz2,  i/3  near  zero,  the  solution  u  =  0  is  unstable  with  dim  ^’“(0)  =  1 . 
Thus,  n(X)  =  3  andAx  =  Wu(0)  U  {«2 }  U  (i/3}is  a  smooth  closed  arc.  At  X  =  9,  there  is  an¬ 
other  bifurcation  at  zero  in  a  direction  independent  of  the  previous  bifurcation.  This  yields 
two  more  solutions  each  unstable-and  the  two  dimensional  set  Ax  mentioned  before. 

For  retarded  functional  differential  equations,  it  does  not  seem  possible  to  have  a  sys¬ 
tem  which  is  a  gradient  system.  However,  there  are  some  equations  which  have  some  of  the 
same  qualitative  properties.  We  discuss  such  an  example  in  some  detail  since  it  brings  out 
many  of  the  difficulties  encountered  in  infinite  dimensions  and  it  also  leads  to  many  inter¬ 
esting  and  specific  unsolved  problems. 

Consider  the  equation 

(9.25)  x(f)  =  -  J°|  a(-0)g{x(t  +  6))dd 

where  g  €  C2(R,  R),  a  £  C2((0,  I] ,  R),  f£g  — *  00  as  |x|  — *  <*,  cr(l)  =  0,  a(s)  >  0.  a(s)  <  0. 
a(s)  >  0. 

Under  these  hypotheses,  every  solution  of  (9.25)  is  defined  and  bounded  for  t  >  0.  If 
C  =  C([-  1,  0) ,  R),  then  the  semigroup  Ta  g(t):  C  — >  C,  /  >  0,  is  well  defined  by  the  rela¬ 
tion  Ta  g(t)<p(0)  =  +  6),  -  1  <  8  <  0,  where  x(0)  is  the  solution  of  (9.25)  through  0. 

If  there  is  an  s  €  [0,  1  ]  such  that  a(s)  >  0,  then  the  w-limit  set  of  any  solution  is  an  equi¬ 
librium  point  of  (9.25);  that  is,  a  zero  of g.  If  a{s)  =  0  for  all  s  G  (0,  1]  (that  is,  a  is  linear), 
then  the  cu-limit  set  of  any  solution  is  either  an  equilibrium  point  or  a  closed  orbit  given  by 
(pt ,  0  <  t  <  1},  where  p(t)  is  a  one-periodic  solution  of  the  ordinary  differential  equation 

(9.26)  y  +  e(0)g(y)  =  0. 


For  a  proof  of  these  results,  see  Hale  (3,  p  122). 
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Let  us  consider  first  the  case  where  d(s)  >  0  for  all  s  £  [0,  1  j  and  g  has  a  finite  num¬ 
ber  of  simple  zeros  {o/(  j  —  1,2,...,*}.  Let 

JP*(a;-)  =  {<p  6  C:  Ta  — *  the  constant  function  Oy  as  t  — ►  ooj ^ 

^“(ay)  =  {0  €  C:  Ta  g(t)<P  — *  the  constant  function  cty  as  t  — ►-<*>}. 

It  is  not  difficult  to  show  that  each  a;.  is  hyperbolic  and  dim  Wu(<xj)  <  1.  Then 

Aa,g  =  U  *“<<>,) 

/=  1 

is  a  union  of  a  finite  number  of  manifolds  and  Ta  g(t)  is  one-to-one  on  Aa  g. 

It  is  possible  to  construct  open  sets  A  C  C2((0,  1  ] ,  R),  G  in  the  Whitney  topology 
on  C2(R,  R)  such  that  Aa  is  upper  semicontinuous  at  each  (a,  g)  £  A  x  G.  The  way  to 
do  this  can  be  discovered  from  the  proof  of  the  above  results  in  Hale  [3]  and  the  method 
of  proof  of  Proposition  9.7.  We  can  now  prove 

Proposition  9.20.  // d(s)  >  0  for  all  s  e  (0,  1  ] ,  the  zeros  {ay,  j  =  1 ,  2 . k]  of 

g  are  simple,  (a,  g)  €-  A  x  G,  and  the  oj-Iimit  set  of  each  Wu{a.j)\  {ay}  is  a  stable  equilibrium, 
then  (a,  g)  is  structurally  stable. 

Proof.  Since  Aa  g  is  upper  semicontinuous  and  uniformly  asymptotically  stable,  one 
can  choose  a  neighborhood  U  of  Aa  g  (as  small  as  desired)  and  a  neighborhood  V  of  (a,  g) 
in  A  x  G  such  that  Tb  h(t)U  C  U,  t  >  0,  ( b ,  h)  G  V.  Since  the  zeros  of  g  are  simple,  we 
have  dim  W“(ay)  <  1 . 

The  assumption  that  the  co-limit  set  of  each  Wu(a-)\{a/]  is  a  stable  equilibrium  im¬ 
plies  Aa  g  is  one  dimensional  as  shown  in  Figure  17.  The  theory  of  the  neighborhood  of  a 
saddle  point  in  Hale  [3]  implies  that  the  flow  in  U  is  given  as  shown  in  Figure  18.  where 
the  vertical  manifolds  have  codimension  1  and  Aa  has  dimension  one.  This  implies  that 
Ab  h  has  the  same  topological  structure  as  Aa  as  well  as  the  same  type  of  flow. 


•  > - • — <  «  >• - • 


Figure  17. 


One  can  now  ask  the  following  interesting  question. 

Question  9.21 .  Suppose  a  is  a  fixed  function  satisfying  d(s)  >  0  for  s6(0,  1 ).  Let 
Gk  =  {g  which  have  exactly  2k  +  1  zeros  which  are  all  simple).  How  many  different  con¬ 
nected  components  of  structurally  stable  systems  are  in  Gkn 

If  *  =  0,  that  is,  g  has  only  one  zero  a,  then  Aa  =  {a} ,  the  constant  function  a  and 
all#  €  G0  are  structurally  stable  and  there  is  only  one  component. 
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If  k  ~  I ,  then  g  €  Gt  has  three  simple  zeros,  o,  <  a2  <  oc3,  with  a, ,  a3  asymptoti¬ 
cally  stable  and  a2  a  saddle  point.  Since  the  unstable  manifold  at  a2  is  smooth,  one  dimen¬ 
sional  and  Aag  is  uniformly  asymptotically  stable,  it  follows  that  Aag  is  a  one  dimensional 
manifold  with  boundary  points  a, ,  a3.  Again,  all  elements  of  G,  are  structurally  stable. 

The  topological  structure  of  Aa  g  is  not  understood  for  the  case  when  g  has  given 
zeros  a,  <  a2  <  a3  <  a4  <a5.  If  we  consider  g  depending  on  some  parameters  p  with 
#(*,  n)  ~  xs  at  p  =  0,  then  Aa  #(p)  is  {0}  for  p  =  0.  There  is  a  center  manifold  in  a  neigh¬ 
borhood  of  x  =  0  which  is  one  dimensional  and  smooth  in  p.  As  p  varies  near  p  =  0,  g(x,  p) 
can  have  as  many  as  S  zeros  which  must  lie  on  this  center  manifold.  Thus,  Aa  #(p)  is  either 
a  point  or  a  one  dimensional  manifold  with  boundary  for  p  small.  As  p  increases,  it  is  con¬ 
ceivable  that  the  topological  structure  of  Aa  g  changes.  Let  us  give  some  intuitive  reasons 
for  why  this  is  possible.  The  author  wishes  to  acknowledge  conversations  with  John  Mallet- 
Paret  and  Shui-Nee  Chow  which  were  of  great  assistance  in  the  remaining  discussion  of  this 
section. 


g 


r  ■> 

t 

\  4  A 

-  -W_  A 

\  V-  A 

X  -/ 

r  ?  6  p 

* 

,  4 

'  >'  ' 
i  4 

h  '  1 

b  J 

r 

Figure  20.  Figure  21. 


Suppsoe  g  has  five  zeros  and  the  general  shape  shown  in  Figure  19.  If  a  is  strictly 
convex  and  very  close  to  the  6-function  at  zero,  then  Aa  g  is  shown  in  Figure  20.  The  com¬ 
plete  flow  near  Aa  g  is  shown  in  Figure  21  with  the  exponential  decay  toward  Aa  g  being 
very  rapid  and  much  greater  than  the  convergence  of  the  flow  on  Aa  g  toward  0. 

If  g'(P)  -  «.  then  the  flow  near  0  is  determined  by  the  roots  of  the  characteristic  equa¬ 
tion 

(9.27)  X  =  -o  J®  a(-0)eKO  <J0. 
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One  can  show  there  are  an  a0  and  a  strictly  convex  function  a0  with  the  property  that  the 
roots  of  this  equation  with  maximum  real  part  corresponds  to  a  double  root.  Change  g 
slightly  so  that  the  zero  5  disappears  and  g  has  only  three  zeros.  Because  of  the  nature  of  a 
stable  node  with  both  eigenvalues  equal,  there  are  two  distinguished  directions  along  which 
(3  can  be  approached  along  the  unstable  manifold  of  y  as  shown  in  Figure  22.  Choose  one 


Figure  22. 

of  these  directions.  Now,  one  should  be  able  to  move  g  back  so  that  the  double  zero  6 
appears  but  the  unstable  manifold  of  5  approaches  /3  from  the  other  distinguished  direction. 
If  this  is  the  case,  then  the  new  Aa  g  is  shown  in  Figure  23.  Now  one  can  further  change  g 


to  make  6  become  two  distinct  zeros  6,,  62  and  the  resulting  Aa  g  is  structurally  stable. 

Thus,  there  are  at  least  two  distinct  connected  components  if  g  has  five  zeros.  To  go  from  one 
of  these  classes  to  the  other  keeping  the  zeros  of  g  simple,  it  is  not  difficult  to  see  that  one 
must  have  a  saddle  connection  at  some  point.  Since  we  have  not  verified  that  it  is  possible 
to  move  g  as  specified  above,  we  state  this  as 

Conjecture  9.22.  There  is  a  saddle  connection  for  some  (a,  g)  with  a  strictly  con¬ 
vex  and  g  having  five  simple  zeros.  The  set  of  (a,  g)  for  which  this  is  true  is  not  generic. 

Another  interesting  conjecture  about  equation  (9.24)  is  the  following 

Conjecture  9.23.  lfd[s)  >  0  for  s  €  [0,  1 J  and  the  zeros  of  g  are  simple,  then  the 
set  W  =  {(V*(a;):  a;.  is  asymptotically  stable}  is  dense  in  C. 

One  could  try  to  prove  this  using  the  same  idea  that  Henry  [!J  employed  to  prove 
the  same  result  for  equation  (9.23).  Let  a  be  a  saddle  point  and  for  any  r>  0,r>  0,  let 
H^  r(a)  =  (0  €  C:  Ta  g(T)<p  e  B*(a),  I  Tag(r)<t>  ~a\<r}.  The  set  W’  r(a)  is  closed.  If  it 
contains  a  ball  B  =  (0,  +  ph:  |h|  =  1, 0  <p  <  r0 },  then  Ta  g{r)B  has  been  flattened  so 
that  it  is  contained  in  a  submanifold  of  codimension  one.  In  particular,  dTag(t)(<t>l  +  ph)/Sp 
at  p  =  0  is  in  the  tangent  space  of  ^(a).  This  function  D^^Tn  g(t)h  is  the  solution  of  the 
linear  variational  equation 
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m  =  a(-6)g\TaJtHi)y(t) 

with  y0  =  h.  If  the  solution  operator  of  the  adjoint  of  this  equation  is  one-to-one,  then 
the  assertion  that  W}  r(a)  contains  a  ball  is  false.  In  fact,  take  a  vector  tj  #  0  such  that 
<77,  <p)  =  0  for  all  0  in  the  tangent  space  to  K^a)  at  Ta  ^(r)0, .  Integrating  the  adjoint  equa¬ 
tion  for  a  function  z,  on  [0,  r]  with  zr  -  rj,  one  obtains,  for  all  h,  (z0,  h)  =  (zT,  h)  = 

(17,  ( Ta  g(f)h)  =  0.  Thus,  z0  =  0.  Since  the  solution  operator  for  the  adjoint  equation 

is  one-to-one,  it  follows  that  77  =  0.  Consequently,  the  set  W’  r(a)  contains  no  ball,  which 
implies  that  W’  r  is  closed  and  nowhere  dense.  Thus,  ^(a)  has  no  interior.  This  will  prove 
the  conjecture. 
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10.  Bifurcation  in  infinite  dimensions 

For  systems  of  parabolic  partial  differential  equations,  retarded  functional  differential 
equations  and  certain  types  of  neutral  equations,  several  problems  of  local  bifurcation  near 
an  equilibrium  point  can  be  discussed  in  essentially  the  same  way  as  for  ordinary  differen¬ 
tial  equations.  In  particular,  the  results  in  §4  on  bifurcation  at  a  zero  eigenvalue  and  in  §5 
on  bifurcation  at  a  focus  extend  verbatum  (see,  for  example,  Marsden  and  McCracken  [1] , 
deOIiveira  and  Hale  [1],  Kielhofer  (lj).  The  center  manifold  theorem  plays  an  important 
role  in  the  theoretical  justification  of  the  results.  For  illustrations  of  the  results  in  §8  for 
the  case  of  a  double  eigenvalue  zero,  see  Carr  (1  ] ,  Howard  [1] . 

Very  little  attention  has  been  devoted  to  the  manner  in  which  local  bifurcation  in¬ 
fluences  the  global  flow  defined  by  the  equations.  In  particular,  how  does  the  maximal 
compact  invariant  set  change  at  bifurcation  points? 

Marchetti  et  al.  (1)  have  some  results  for  locally  compact  spaces.  More  specifically, 
they  relate  the  change  in  stability  of  A{  to  the  appearance  of  new  invariant  sets.  It  would 
be  interesting  to  extend  this  to  Banach  spaces  for  the  type  of  semigroups  mentioned  in  the 
previous  section.  For  the  gradient  equation  (9.23),  we  have  discussed  in  the  previous  section 
the  results  of  Henry  (1),  [2)  on  the  manner  in  which  the  maximal  compact  invariant  set 
changes  with  a  parameter. 

For  the  more  general  equation  (9.21)  with  /  depending  only  on  u  and  boundary  con¬ 
ditions  (9.19),  Brunovsky  and  Chow  [1]  have  shown  the  following  interesting  fact:  generi- 
cally  in  /€  C2( R,  R),  the  bifurcations  always  are  saddle-node  type;  that  is,  a  saddle  and 
node  coalesce  and  disappear.  In  particular,  this  has  the  following  interesting  implication 
for  second  order  ordinary  differential  equations.  The  period  u>(b,  /)  of  the  solution  of 
uxx  +  /(M)  =  9  through  the  point  u( 0)  =  a,  ux(0)  =  0  is  a  Morse  function  generically 
in  /  Smoller  and  Wasserman  (I)  have  also  discussed  this  function  in  detail  for  specific 
functions  / 

For  the  equation  x(f)  =  ~/(x(r  -  1)),  where  /  is  piecewise  linear  with  a  finite  number 
of  jump  discontinuities,  Walther  (1  j  and  Chow  and  Walther  [1]  have  discussed  the  set  Af 
and  the  dynamics  on  Af  in  some  detail.  The  hypotheses  on  /serve  as  a  model  for  the  equa¬ 
tion  of  Wright  (see  Nussbaum  (21). 

For  equation  (9. 25),  we  discussed  how  the  set  Aa  could  bifurcate  by  the  creation  of 
an  arc  connecting  two  saddle  points.  This  is  always  a  bifurcation  caused  by  the  global  prop¬ 
erties  of  the  flow. 
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The  analogue  of  the  Generic  Hopf  Bifurcation  Theorem  (Theorem  5.6)  for  infinite 
dimensional  systems  is  prevalent  in  the  literature.  It  seems  reasonable  to  say  that  it  is  the 
generic  situation.  On  the  other  hand,  we  must  be  careful  because  there  are  so  many  differ¬ 
ent  ways  to  model  a  problem  in  infinite  dimensions.  For  example,  in  equation  (9.25),  the 
vector  field  should  be  considered  determined  by  the  function  a(s)g(x),  s  €  [0,  1] ,  x  €  R; 
that  is,  as  a  product  of  a  function  of  s  and  a  function  of  x.  This  is  certainly  an  easier  model 
to  discuss  than  one  which  uses  a  general  function  h(s,  x ).  However,  it  is  still  feasible  that 
all  interesting  qualitative  behavior  of  solutions  can  be  obtained  by  considering  the  integrand 
as  a  product  of  two  functions.  For  equation  (9.25),  the  function  a0,  g,  satisfying  xgfx)  >  0 
for  x  ^  0,  g'(0)  =  1,  a0(s)  =  4jt2(  1  -  s),  is  a  bifurcation  point  for  every  g.  The  linear  vari¬ 
ational  equation  at  zero  has  all  eigenvalues  with  negative  real  parts  except  the  purely  imagi¬ 
nary  ones  ±2 id.  Following  the  procedure  in  §5,  one  can  compute  the  bifurcation  function 
G(b,  a,  g)  for  the  periodic  orbits  near  zero  with  the  amplitude  of  the  periodic  orbits  being 
given  approximately  by  b.  For  the  Generic  Hopf  Bifurcation  Theorem  to  be  applicable,  one 
should  have  G(b,  a0,  g)  =  P„ogb3  +  oflil3)  as  131  — *■  0,  0aQg  a*  0.  However,  one  can  show 
that  0aQ  g  =  0  for  all  g  in  the  above  class!  Thus,  the  bifurcation  from  the  focus  is  never  a 
bifurcation  of  order  one.  This  could  not  happen  if  <j(s)g(x)  were  replaced  by  a  general  func¬ 
tion  h(s,  x).  For  a  detailed  discussion  of  this  point  as  well  as  the  nature  of  the  bifurcation 
point,  see  Hale  [7] . 

In  the  previous  section,  we  noted  that  it  was  not  known  if  the  Kupka-Smale  systems 
were  generic  for  x(r)  =  f(x(t  -  1)).  This  is  a  problem  of  the  same  type  as  in  the  previous 
paragraph  where  severe  restrictions  are  imposed  on  the  vector  field.  We  have  seen  before 
that  the  same  difficulties  occur  in  parabolic  equations  when  the  nonlinearities  depend  only 
on  the  dependent  variable. 

Bifurcation  of  a  periodic  orbit  to  a  torus  for  infinite  dimensional  systems  has  been 
discussed  in  some  detail  (see  Iooss  and  Joseph  {1]  for  references).  Holmes  and  Marsden  [1] 
show  that  equation  (9.10)  with  a  periodic  forcing  can  have  homoclinic  points. 

There  are  some  bifurcation  problems  that  are  unique  to  infinite  dimensions.  For  ex¬ 
ample,  consider  the  scalar  parabolic  equation 

(10.1)  ut  =  Au  +  f(u)  in  S2,  du/dn  =  0  on  dfi 

where  £1  is  a  bounded  open  connected  set  in  R"  with  smooth  boundary  3 fl.  If  W  =  L2(£l,  R), 
suppose  this  equation  generates  a  strongly  continuous  semigroup  on  Wa ,  a  €  (0,  1).  Then 
every  bounded  orbit  is  precompact  and  the  co-limit  set  of  any  bounded  orbit  can  be  shown 
to  belong  to  the  set  E  of  equilibrium  solutions  of  (10.1);  that  is,  the  set  of  solutions  of  the 
equation 

(10.2)  Au+/(u)  =  0  inn,  3u/3n  =  0  on  3n 

(see  Matano  (!]).  For  n  —  1,  Matano  (2]  has  shown  that  the  co-limit  of  any  bounded  orbit 
is  a  single  point  in  E.  The  same  result  for  arbitrary  n  is  contained  in  Hale  and  Massatt  [1] 
with  the  additional  restriction  that  every  point  in  a  connected  component  of  E  (which  is  not 


DYNAMIC  BIFURCATION  THI  OKY 


71 


a  single  point)  has  zero  as  a  simple  eigenvalue  of  the  corresponding  linear  variational  equa¬ 
tion.  In  particular,  if  all  solutions  are  bounded  and  the  set  of  equilibrium  points  is  finite, 
then  Af  is  compact,  Af  =  U,  Wu(atj)  where  ti'u(a/)  is  the  unstable  manifold  of  the  equilib¬ 
rium  point  a;. 

For  12  convex,  Casten  and  Holland  [1 1 ,  Matano  (1|  have  shown  that  the  only  stable 
equilibrium  points  of  (10.1)  are  spatially  homogeneous;  that  is,  they  correspond  to  constant 
functions-the  zeros  of  f.  Chafee  [2 ]  previously  has  proved  this  result  for  n  =  1.  Bardos, 
Matano  and  Smoller  [1]  have  proved  a  similar  result  for  one  equation  of  the  form  (1 0.1) 
coupled  with  some  ordinary  differential  equations.  Casten  and  Holland  |1],  Matano  |1| 
also  have  proved  this  type  of  instability  for  some  other  special  domains  £2. 

If  n  >  2 ,f(u)  =  u  -  u3  (or  any  function  with  similar  properties),  Matano  [1|  has 
shown  that  there  are  stable  equilibrium  points  of  (10.1)  which  arc  not  spatially  homogene¬ 
ous  if  the  region  £2  has  certain  properties.  For  example,  it  is  sufficient  to  have  £2  consist 
of  disjoint  open  convex  connected  sets  £2,,  £22  joined  by  a  channel  £23  which  is  not  too 
wide  compared  to  the  size  of  £2,,  £22  (see  Figure  24). 
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This  leads  to  the  following  interesting  bifurcation  problem.  Let  /(«)  -  \u  -  u  \ 

X  e  R,  and  let  £2f,  e  €  R,  be  a  simply  connected  open  set  in  R"  with  smooth  boundary,  £20 
convex.  Let  TK  e(t)  be  the  semigroup  generated  by  (10. 1)  for  this/and  £2( .  Then  T0  0(f) 
has  a  unique  equilibrium  point  u  =  0  which  is  uniformly  asymptotically  stable  and  attracts 
bounded  sets  of  Wa .  The  maximal  compact,  invariant  set  AK  c  for  (X,  e)  =  (0,  0)  is  shown 
in  Figure  25a.  For  X  >  0  and  small,  there  are  only  three  equilibrium  points  0,  ±X'/!  and 
■4^  0  as  shown  in  Figure  25b.  Now  suppose  that  the  region  £2,  with  increasing  e  becomes 
nonconvcx  so  that  the  second  eigenvalue  p2(e)  for  the  Neumann  problem  lor  the  Laplacian 
satisfies  /’( 0)  =  X  <  p2(e)  for  0  <  e  <  t0.J\ 0)  =  p2(f0)  and  J  '(0)  >  p2(e)  for  e  >  e0. 

The  origin  becomes  a  bifurcation  point  at  e  -  e0,  and  AK  t  for  e  >  ca  in  a  neighborhood  of 
e0  is  shown  in  Figure  25c.  The  solutions  bifurcating  from  0  are  spatially  nonhomogeneous. 
Suppose  tii at  £2e  as  e  — *  00  approaches  the  set  £2„  =  £2,  U  £22  U  /.  where  £2, ,  £22,  /.  arc 
disjoint,  £2, ,  £22  are  open  convex  connected  sets  and  L  is  an  (rt  -  I )  dimensional  closed 
manifold  of  codimension  I  joining  £2,  to  £22  (in  R2,  an  arc).  There  will  be  some  point  r, 
such  that  the  conditions  for  Matano’s  theorem  arc  satisfied  for  e  >  e,.  There  will  be  a  sta¬ 
ble  spatially  nonhomogeneous  equilibrium  point.  Thus,  the  set  AK  f  has  had  to  undergo 
another  bifurcation.  This  could  not  have  occurred  at  the  equilibrium  points  0,  t  X  /;  since 
±  X/;  arc  uniformly  asymptotically  stable  and  0  has  a  two  dimensional  unstable  manifold. 
Thus,  it  is  reasonable  to  conjecture  that  the  unstable  spatially  nonhomogeneous  solutions 
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undergo  a  bifurcation  as  shown  in  Figure  25d  at  least  if  the  family  of  regions  fif  are  chosen 
appropriately.  The  bifurcation  diagram  is  conjectured  to  be  the  one  shown  in  Figure  25e  for 
certain  fle.  All  the  solutions  are  spatially  nonhomogeneous  except  0,  ±  and  the  stable 
ones  occur  as  a  secondary  bifurcation.  It  would  be  interesting  to  obtain  these  curves 
numerically.  Hale  and  Vegas  [1]  showed  there  is  an  e0  >  0  such  that  the  bifurcation 
diagram  for  e  >  e0  is  the  one  shown  in  Figure  25e  provided  that  the  region  satisfies 
some  conditions  in  e.  The  most  interesting  condition  is  that  the  third  eigenvalue  of  the 
laplacian  on  f le  is  bounded  away  from  zero  for  all  e.  This  restriction  arises  because  the 
analysis  consists  in  treating  the  problem  as  a  bifurcation  from  a  double  eigenvalue  at  e  =  0. 
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If  the  bifurcation  diagram  is  the  one  depicted  in  Figure  25e  and  one  introduces  an¬ 
other  parameter  in  the  vector  field  /,  then  one  should  be  able  to  break  all  of  the  symmetries 
in  the  problem  and  obtain  a  jumping  from  one  state  to  another  as  in  buckling  problems. 

Chafee  and  Infante  ( 1  ]  have  investigated  the  Dirichlet  problem 

(10.3)  ut  =  Au  +  \f(u)  in  fi,  u  =  0  on  bSl 

for  Q.  =  [0,  it)  ,  X  e  R.  Under  suitable  conditions  on  /  including  uf(u )  >  0  for  u  #  0,  they 
show  that  the  only  stable  equilibrium  solution  is  u  =  0  if  0  <  X  <  X, ,  and  either  of  two 
functions  u±(X)  for  X  >  X,  where  u+(X)  is  positive  on  (0,  tt)  and  u.  (X)  is  negative  on  (0,  it). 
For  fl  £.  R",  one  should  encounter  a  much  more  interesting  type  of  behavior  by  changing 
the  shape  of  the  region  SI  as  we  did  with  the  Neumann  boundary  conditions.  To  the  author’s 
knowledge,  this  problem  has  not  been  investigated. 

It  would  also  be  interesting  to  discuss  the  same  equation  with  mixed  boundary  condi¬ 
tions  varying  both  the  boundary  conditions  and  the  region. 

For  systems  of  reaction  diffusion  equations  with  one  space  variable  and  Neumann 
boundary  conditions,  one  can  also  obtain  stable  equilibrium  solutions  which  are  not  spatially 
homogeneous.  This  is  a  bifurcation  problem  where  one  creates  an  instability  of  the  zero 
solution  by  making  the  linear  approximation  have  a  zero  eigenvalue  by  varying  the  diffusion 
coefficients  and  the  linear  coupling  terms  (for  references  and  a  detailed  discussion,  see  Fife 
[1]).  This  corresponds  to  primary  bifurcation  from  the  trivial  solution  whereas  the  stable 
solution  in  the  previous  discussion  was  created  through  a  secondary  bifurcation. 

For  reaction  diffusion  equations  in  an  unbounded  domain,  there  are  some  very  inter¬ 
esting  bifurcation  problems  associated  with  traveling  wave  solutions.  Due  to  limitations  in 
space,  we  can  only  refer  to  Fife  [1]  for  references. 

Another  problem  which  becomes  important  in  infinite  dimensions  is  the  discussion  of 
semigroups  Tx(t)  which  are  continuous  in  a  parameter  X,  but  are  not  continuously  differentia¬ 
ble  in  X.  In  particular,  it  is  possible  to  discuss  how  invariant  sets  change  with  X,  how  the  stable 
manifolds  change  with  X,  etc.?  In  the  qualitative  theory,  these  arc  fundamental  questions.  If 
7\(r)  has  a  maximal  compact  invariant  set  then  all  of  the  interesting  properties  of  the  orbits 
of  rA(f)  are  determined  by  the  behavior  of  the  orbits  on  .  The  orbits  on  Ax  are  bounded 
and  defined  on  (-°°,  «*>).  As  we  indicated  earlier,  for  some  important  types  of  semigroups,  the 
fact  that  orbits  on  Ax  are  defined  on  (-<*>, “>)  implies  that  Tx(tpc  forxCA  is  continuously  dif¬ 
ferentiable  in  t;  that  is,  the  elements  in  <4A  arc  “smoother"  than  a  generic  clement  of  the  under¬ 
lying  space  X.  Can  this  additional  smoothness  on  A  A  be  exploited  to  discuss  how  A  A  varies  with 
X?  For  differential  difference  equations  with  the  parameter  X  being  the  delays,  we  indicate  why 
this  is  feasible.  No  corresponding  results  have  been  obtained  for  partial  differential  equations. 

The  following  lemma  is  crucial  to  the  analysis.  It  asserts  that  the  fixed  points  of  a 
map  can  be  shown  to  be  continuously  differentiable  in  a  parameter  X  by  requiring  only  the 
derivative  of  the  map  be  differentiable  with  respect  to  X  on  the  fixed  point  set.  This  lemma 
was  stated  in  Hale  [3|  with  the  omission  of  the  obviously  necessary  fact  that  the  derivative 
in  (iv)  below  is  continuous. 


74 


JACK  K.  HALF 


Lkmma  10.1 .  l£t  F  be  a  closed  subset  of  a  Banach  space  X,  Int  F  #  0,  A  be  an 
open  subset  of  a  Banach  space  Y.  where  Int  F  denotes  the  interior  of  F.  Assume  that  T: 

F  x  A  — *  F,  (x,  X)  — - ►  T[x,  X),  satisfies  the  following  set  of  hypotheses: 

(i)  7\x,  •):  A  —*  F  is  continuous; 

(ii)  T\- ,  X):  F  — ►  F  is  continuous  and  has,  for  each  \  a  unique  fixed  point  x(X) 
which  is  continuous  in  X; 

(iii)  if  x(A)  =  F{ ,  then  1\x,  X)  is  continuously  differentiable  in  X  for  (x,  X)  G  F,  x  A ; 

(iv)  there  are  an  open  set  F{  C  X,  F  C  F2,  and  a  6  €  [0,  1)  such  that  the  derivative 
of  7\x,  X)  with  respect  to  x  is  continuous  and  has  norm  <  S  for  all  (x,  X)  6  F2  x  A. 

Then  the  fixed  point  x(  X),  X  6  A ,  of  7\x,  X)  is  continuously  differentiable  in  X. 

Proof.  The  proof  is  an  adaptation  of  the  usual  proof  of  differentiability  with  re¬ 
spect  to  parameters.  We  merely  give  an  outline  (see  P.  Lima  [  1 )  for  details).  Let  TK(x,  X) 
=  dT(x,  X)/3X,  Tx(x,  X)  =  dT(x,  X)/3X  and  consider  the  equation 

z  -  r^WX),  X)z  =  rx(x(X),  \)h 

for  given  h  €  Y.  This  equation  has  a  unique  solution  z(h,  X)  which  is  linear  in  h.  If 
z(h,  X)  =  z(\)h,  then  one  easily  shows  that  z(X)  is  continuous  in  X. 

To  complete  the  proof,  one  must  show  that  z(X)  satisfies  the  definition  of  the  deriva¬ 
tive  of  x(X);  that  is, 

cu  ds  x(X  +  h)  -  *(X)  -  z( \)h  =  oil/tl) 

as  |/i|  — ♦  0.  An  easy  computation  shows  that  (/  -  Tx(x(X),  X))]u>  -  R( u>,  X,  h)  =  o(|/t|) 
as  |/i|  — *  0  where  R(y,  X.  h)  =  l\x  +  y,  X)  -  T\x,  X)  -  Tx(x,  X).v  is  o(|^i)  as  |v|  — *  0. 

One  now  proves  that  this  implies  cj  =  o(|/r|)  as  |A|  — *  0. 

We  illustrate  the  application  of  this  theorem  to  prove  the  Hopf  Bifurcation  Theorem 
with  respect  to  the  delays  for  a  differential  difference  equation.  We  impose  more  restrictions 
than  in  Hale  [8|  because  the  proof  here  will  be  based  on  the  center  manifold  theorem.  By 
showing  there  is  an  asymptotically  stable  local  center  manifold  of  dimension  two  which  is 
continuously  differentiable  in  X,  the  problem  is  reduced  to  the  usual  one  in  ordinary  differ¬ 
ential  equations. 

Suppose  ft  C  R*  is  an  open  set  (the  parameter  space),  C  =  CQ-r,  0],  R").  /:  ft  x 
C  — ►  R",  ft  x  C  — ►  R"  are  continuous,  L(a)<l>  is  linear  in  <t>,f(a ,  <f>)  has  continuous  first 

and  second  derivatives  in  <t>,  /(a ,  0)  =  0,  3/(a ,  0)/3a  =  0.  IX  he  notation  xt(6)  =  x(t  +  3), 
-r  <0  <  0,  consider  the  equation 

00.4)  x(t)  =  L(a)x,  +  /(a,  „,). 

The  first  hypothesis  is: 

(H, )  The  matrix  A(a,  X)  =  X/  -  i.(a)cK'i,  where  /  is  the  identity  matrix,  is  continu¬ 
ously  differentiable  in  a,  there  is  a  pair  of  simple  purely  imaginary  roots  t  iV0,  i>0  >  0.  for 
a  =  or0  and  all  other  roots  of  the  characteristic  equation 
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(>0.5)  det  Ma,  X)  =  0 

for  a  -  a0  have  negative  real  parts. 

Under  hypothesis  (H, ),  there  is  a  6  >  0  and  a  simple  characteristic  root  \(a)  which  is 
C1  in  or  for  |or  -  a0|  <  S  and  X(a0)  =  iv0.  Then  Re  X(a)  =  (a  -  o0)  ■  f(cr)  +  o(|a  -  ct0|) 
as  or  *■  a0,  where  f(a)  £  R*  is  Cl  in  or.  Our  next  hypotheses  are: 

(H2)  f(a0)  *  0. 

(H3)  Considering  L(a)0,  /(a,  <j>)  restricted  to  R*  x  C'([-r,  0] ,  R"),  they  have  a 
derivative  in  or  which  is  continuous  in  or,  <t>  in  the  topology  of  R*  x  C'([-r,  0] ,  R"). 

For  a  a  scalar,  hypothesis  (H2)  implies  there  are  two  eigenvalues  crossing  the  imaginary 
axis  at  a  =  a0 .  Hypothesis  (H3)  does  not  imply  that  L(a)<t>,  /(or,  0)  have  a  derivative  in  a 
for  (or,  0)  £  C.  For  example,  if  a  £  (0,  r),  the  function  L(a)0  =  0(-a)  satisfies  hypothesis 
(Hj).  This  function  clearly  has  no  continuous  derivative  in  a  for  (a,  0)  £  R  x  C  since  the 
linear  operator  L(a)  is  not  continuous  in  a  in  the  operator  topology.  In  Hale  [7] ,  (H3) 
was  stated  incorrectly.  The  hypothesis  did  not  state  that  the  continuity  should  be  in  a 
and  0. 

Theorem  10.2.  If  (H,),  (H2),  (H3)  are  satisfied,  then  there  is  an  e  >  0  such  that  for 
<r  £  R,  |a|  <  e,  there  is  a  Cx -manifold  ra  €  R*  of  codimension  one,  ra  is  Cl  in  a,  ro  = 

{a  £  R*:  Re  X(a)  =  0,  |a  -  a0|  <  e)  and,  for  every  a  €  T^,  there  are  a  function  co(a,  a), 
an  co(a,  afperiodic  function  x*(a,  a)  which  is  Cl  in  a  and  a,  w(a0,  0)  =  w0  =  2ir/v0, 
x*(a,  0)  =  0  and 

0)(f)  =  ay  cos  v0t  +  o(|a|)  as  )o|  — *•  0 
where  y  cos  v0t  is  a  solution  of  x(t)  =  L(a0)xt  with  y  £  R",  |y|  =  ] . 

To  indicate  the  proof,  decompose  the  space  C  by  the  eigenvalues  (X(a),  X(a))  as 
c  ~  Qa  where  Pa  is  two  dimensional  and  spanned  by  the  solutions  ofx(t)  =  L(a)xt 
corresponding  to  the  eigenvalue  X(a),  X(a).  Let  0a  be  an  n  x  2  matrix  whose  columns  are 
a  basis  for  Pa.  If  Ta(t)  is  the  semigroup  generated  by  this  linear  equation  and  xt  =  <t>n  y(t) 

+  zt  with  y  €  R2,  z,  €  Q,  then  equation  (10.4)  is  equivalent  to 

y  -  Bay  +  caf(a,  *ay  +  zt), 

zt  -  Ta(0<l'  +  fg  Ta(t  -  s)*f(«,  *„y(s)  +  Z,)ds, 

where  Ba,  Ca  are  two-by-two  matrices  C1  in  a,  with  the  eigenvalues  of  Ba  being  X(a), 

X(a),  0  =  *0,  ♦  is  the  projection  of  the  n  x  n  matrix  X0  onto  Qa,  X0(9)  =  0  for  6  <  0, 
*0(°)  =  /  (see  Hale  [3J  for  details).  There  is  an  estimate  || Ta(t  -  s)ICal  II  <  Ke~at,  t>  0, 
for  some  positive  constants  K,  a. 

To  construct  a  local  center  manifold  by  means  of  simple  looking  formulas,  let  us  sup¬ 
pose  /  and  its  derivative  are  bounded  on  fi  x  C.  If  T  =  C'(R2,  Q),  (y,  h)  £  R2  x  r,  let 
i(t.  y.  h),  f(0,  y,  h)  =  y,  be  the  solution  of  the  equation  f  =  +  Caf(a,  <J>Qf  +  h({)) 
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and  define  the  operator 

(10.6)  A'(a,  y)h  =  /_°_  TJ-s)*A*.  <*a{(s,  y,  h))  +  Mf(*.  y.  h))ds. 

A  fixed  point  of  K( a,  h )  is  a  center  manifold  of  (10.4).  An  application  of  the  Implicit 
Function  Theorem  yields  a  center  manifold  Ma  =  {(>,  4>)'.  &  =  h( a,  >)}  with  the  function 
h(a.  y)  continuous  in  a,  y  and  having  as  many  continuous  derivatives  in  y  as  the  function 
/(a,  <t> )  has  in  <t>  (finite,  of  course),  h(a,  0)  =  0.  The  flow  on  the  center  manifold  is  deter¬ 
mined  by  the  ordinary  equation 

00.7)  y  =  Bay  +  CaAa,  <t>a y(l)  +  h(a, y(r))) 

and  is  given  by  (X0.  M®.  Xf)))  All  solutions  with  initial  values  on  Ma  are  defined  for 
t  6  (-<*»,  »).  If  we  were  considering  only  a  local  center  manifold,  we  need  here  t  6  [-r,  (J), 
0  >  0,  which  would  be  no  loss  in  generality.  If  y(t)  is  a  solution  of  (10.7),  this  implies  that 
xt  -  *ay(t)  +  h(a,  X0)  is  defined  for  /€(-«>,  °°).  Since  xt  satisfies  (10.4)  and  4>a(0)  is 
C1  in  0  this  implies  h(a,  XO))(0)  >s  continuously  differentiable  in  6  for  6  6  [-r,  0] .  Thus, 
the  fixed  point  set  F,  for  K(a,  y)  must  consist  of  continuously  differentiable  functions.  Hy¬ 
pothesis  (H3)  implies  that  condition  (iv)  of  Lemma  10.1  is  satisfied.  Thus  one  can  conclude 
that  h(a,  y)  is  C1  in  (a,  >>)  (or  Ck  if  all  hypotheses  are  satisfied  for  derivatives  up  through 
order  k )  and  that  the  vector  field  in  (10.7)  is  C1  in  a,  y.  The  proof  of  Theorem  10.2  is 
completed  by  applying  the  corresponding  result  for  ordinary  equations. 

Some  very  interesting  bifurcations  occur  in  problems  with  several  delays  (see  Hale  [8] , 
Nussbaum  [1]). 

Marsden  and  McCracken  [1,  p.  255]  have  a  version  of  the  Hopf  Bifurcation  Theorem 
for  general  semigroups  of  transformations  satisfying  some  smoothness  properties.  Using 
Lemma  10.1,  it  should  be  possible  to  improve  those  results. 
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